A hierarchical approach to model parameter optimization for developmental systems.

In the context of Systems Biology, computer simulations of gene regulatory networks provide a powerful tool to validate hypotheses and to explore possible system behaviors. Nevertheless, modeling a system poses some challenges of its own: especially the step of model calibration is often difficult due to insufficient data. For example when considering developmental systems, mostly qualitative data describing the developmental trajectory is available while common calibration techniques rely on high-resolution quantitative data. Focusing on the calibration of differential equation models for developmental systems, this study investigates different approaches to utilize the available data to overcome these difficulties. More specifically, the fact that developmental processes are hierarchically organized is exploited to increase convergence rates of the calibration process as well as to save computation time. Using a gene regulatory network model for stem cell homeostasis in Arabidopsis thaliana the performance of the different investigated approaches is evaluated, documenting considerable gains provided by the proposed hierarchical approach.

Reduction des risques au cabinet du généraliste chez les consommateurs de cocaine: une perspective motivationnelle [Risk reduction among cocaine users: a motivational approach in the GP's office]

Today, cocaine use is a public health issue. Cocaine is a powerfully addictive stimulant drug which use is increasing among some part of the population. After a brief description of the physical and psychological effects of cocaine use, the article presents a motivational way for general practitioners to deal with risk-reduction issues. Based on the Transtheoretical Model of human behavior change and providing clinical examples, the article focuses particularly on the two earliest stages of change: "pre-contemplation" and "contemplation".

A model of cellular dosimetry for macroscopic tumors in radiopharmaceutical therapy.

PURPOSE: In the radiopharmaceutical therapy approach to the fight against cancer, in particular when it comes to translating laboratory results to the clinical setting, modeling has served as an invaluable tool for guidance and for understanding the processes operating at the cellular level and how these relate to macroscopic observables. Tumor control probability (TCP) is the dosimetric end point quantity of choice which relates to experimental and clinical data: it requires knowledge of individual cellular absorbed doses since it depends on the assessment of the treatment's ability to kill each and every cell. Macroscopic tumors, seen in both clinical and experimental studies, contain too many cells to be modeled individually in Monte Carlo simulation; yet, in particular for low ratios of decays to cells, a cell-based model that does not smooth away statistical considerations associated with low activity is a necessity. The authors present here an adaptation of the simple sphere-based model from which cellular level dosimetry for macroscopic tumors and their end point quantities, such as TCP, may be extrapolated more reliably. METHODS: Ten homogenous spheres representing tumors of different sizes were constructed in GEANT4. The radionuclide 131I was randomly allowed to decay for each model size and for seven different ratios of number of decays to number of cells, N(r): 1000, 500, 200, 100, 50, 20, and 10 decays per cell. The deposited energy was collected in radial bins and divided by the bin mass to obtain the average bin absorbed dose. To simulate a cellular model, the number of cells present in each bin was calculated and an absorbed dose attributed to each cell equal to the bin average absorbed dose with a randomly determined adjustment based on a Gaussian probability distribution with a width equal to the statistical uncertainty consistent with the ratio of decays to cells, i.e., equal to Nr-1/2. From dose volume histograms the surviving fraction of cells, equivalent uniform dose (EUD), and TCP for the different scenarios were calculated. Comparably sized spherical models containing individual spherical cells (15 microm diameter) in hexagonal lattices were constructed, and Monte Carlo simulations were executed for all the same previous scenarios. The dosimetric quantities were calculated and compared to the adjusted simple sphere model results. The model was then applied to the Bortezomib-induced enzyme-targeted radiotherapy (BETR) strategy of targeting Epstein-Barr virus (EBV)-expressing cancers. RESULTS: The TCP values were comparable to within 2% between the adjusted simple sphere and full cellular models. Additionally, models were generated for a nonuniform distribution of activity, and results were compared between the adjusted spherical and cellular models with similar comparability. The TCP values from the experimental macroscopic tumor results were consistent with the experimental observations for BETR-treated 1 g EBV-expressing lymphoma tumors in mice. CONCLUSIONS: The adjusted spherical model presented here provides more accurate TCP values than simple spheres, on par with full cellular Monte Carlo simulations while maintaining the simplicity of the simple sphere model. This model provides a basis for complementing and understanding laboratory and clinical results pertaining to radiopharmaceutical therapy.

Dispensable role of myeloid differentiation primary response gene 88 (MyD88) and MyD88-dependent toll-like receptors (TLRs) in a murine model of osteoarthritis.

OBJECTIVES: The aim of our study was to evaluate the role of cell-membrane expressed TLRs and the signaling molecule MyD88 in a murine model of OA induced by knee menisectomy (surgical partial removal of the medial meniscus [MNX]). METHODS: OA was induced in 8-10weeks old C57Bl/6 wild-type (WT) female (n=7) mice and in knockout (KO) TLR-1 (n=7), -2 (n=8), -4 (n=9) -6 (n=5), MyD88 (n=8) mice by medial menisectomy, using the sham-operated contralateral knee as a control. Cartilage destruction and synovial inflammation were evaluated by knee joint histology using the OARSI scoring method. Apoptotic chondrocytes and cartilage metabolism (collagen II synthesis and MMP-mediated aggrecan degradation) were analyzed using immunohistochemistry. RESULTS: Operated knees exhibited OA features at 8weeks post-surgery compared to sham-operated ones. In menisectomized TLR-1, -2, -4, and -6 deficient mice, cartilage lesions, synovial inflammation and cartilage metabolism were similar to that in operated WT mice. Accordingly, using the same approach, we found no significant protection in MyD88-deficient mice in terms of OA progression as compared to WT littermates. CONCLUSIONS: Deficiency of TLRs or their signalling molecule MyD88 did not impact on the severity of experimental OA. Our results demonstrate that MyD88-dependent TLRs are not involved in this murine OA model. Moreover, the dispensable role of MyD88, which is also an adaptor for IL-1 receptor signaling, suggests that IL-1 is not a key mediator in the development of OA. This latter hypothesis is strengthened by the lack of efficiency of IL-1β antagonist in the treatment of OA.

Visual functions in the healthy and schizophrenic brain : from stimulus control to illusory perceptions using electrical neuroimaging

ABSTRACT (FRENCH)Ce travail de thèse basé sur le système visuel chez les sujets sains et chez les patients schizophrènes, s'articule autour de trois articles scientifiques publiés ou en cours de publication. Ces articles traitent des sujets suivants : le premier article présente une nouvelle méthode de traitement des composantes physiques des stimuli (luminance et fréquence spatiale). Le second article montre, à l'aide d'analyses de données EEG, un déficit de la voie magnocellulaire dans le traitement visuel des illusions chez les patients schizophrènes. Ceci est démontré par l'absence de modulation de la composante PI chez les patients schizophrènes contrairement aux sujets sains. Cette absence est induite par des stimuli de type illusion Kanizsa de différentes excentricités. Finalement, le troisième article, également à l'aide de méthodes de neuroimagerie électrique (EEG), montre que le traitement des contours illusoires se trouve dans le complexe latéro-occipital (LOC), à l'aide d'illusion « misaligned gratings ». De plus il révèle que les activités démontrées précédemment dans les aires visuelles primaires sont dues à des inférences « top- down ».Afin de permettre la compréhension de ces trois articles, l'introduction de ce manuscrit présente les concepts essentiels. De plus des méthodes d'analyses de temps-fréquence sont présentées. L'introduction est divisée en quatre parties : la première présente le système visuel depuis les cellules retino-corticales aux deux voix du traitement de l'information en passant par les régions composant le système visuel. La deuxième partie présente la schizophrénie par son diagnostic, ces déficits de bas niveau de traitement des stimuli visuel et ces déficits cognitifs. La troisième partie présente le traitement des contours illusoires et les trois modèles utilisés dans le dernier article. Finalement, les méthodes de traitement des données EEG seront explicitées, y compris les méthodes de temps-fréquences.Les résultats des trois articles sont présentés dans le chapitre éponyme (du même nom). De plus ce chapitre comprendra les résultats obtenus à l'aide des méthodes de temps-fréquenceFinalement, la discussion sera orientée selon trois axes : les méthodes de temps-fréquence ainsi qu'une proposition de traitement de ces données par une méthode statistique indépendante de la référence. La discussion du premier article en montrera la qualité du traitement de ces stimuli. La discussion des deux articles neurophysiologiques, proposera de nouvelles d'expériences afin d'affiner les résultats actuels sur les déficits des schizophrènes. Ceci pourrait permettre d'établir un marqueur biologique fiable de la schizophrénie.ABSTRACT (ENGLISH)This thesis focuses on the visual system in healthy subjects and schizophrenic patients. To address this research, advanced methods of analysis of electroencephalographic (EEG) data were used and developed. This manuscript is comprised of three scientific articles. The first article showed a novel method to control the physical features of visual stimuli (luminance and spatial frequencies). The second article showed, using electrical neuroimaging of EEG, a deficit in spatial processing associated with the dorsal pathway in chronic schizophrenic patients. This deficit was elicited by an absent modulation of the PI component in terms of response strength and topography as well as source estimations. This deficit was orthogonal to the preserved ability to process Kanizsa-type illusory contours. Finally, the third article resolved ongoing debates concerning the neural mechanism mediating illusory contour sensitivity by using electrical neuroimaging to show that the first differentiation of illusory contour presence vs. absence is localized within the lateral occipital complex. This effect was subsequent to modulations due to the orientation of misaligned grating stimuli. Collectively, these results support a model where effects in V1/V2 are mediated by "top-down" modulation from the LOC.To understand these three articles, the Introduction of this thesis presents the major concepts used in these articles. Additionally, a section is devoted to time-frequency analysis methods not presented in the articles themselves. The introduction is divided in four parts. The first part presents three aspects of the visual system: cellular, regional, and its functional interactions. The second part presents an overview of schizophrenia and its sensoiy-cognitive deficits. The third part presents an overview of illusory contour processing and the three models examined in the third article. Finally, advanced analysis methods for EEG are presented, including time- frequency methodology.The Introduction is followed by a synopsis of the main results in the articles as well as those obtained from the time-frequency analyses.Finally, the Discussion chapter is divided along three axes. The first axis discusses the time frequency analysis and proposes a novel statistical approach that is independent of the reference. The second axis contextualizes the first article and discusses the quality of the stimulus control and direction for further improvements. Finally, both neurophysiologic articles are contextualized by proposing future experiments and hypotheses that may serve to improve our understanding of schizophrenia on the one hand and visual functions more generally.

iLOGP: A Simple, Robust, and Efficient Description of n-Octanol/Water Partition Coefficient for Drug Design Using the GB/SA Approach.

The n-octanol/water partition coefficient (log Po/w) is a key physicochemical parameter for drug discovery, design, and development. Here, we present a physics-based approach that shows a strong linear correlation between the computed solvation free energy in implicit solvents and the experimental log Po/w on a cleansed data set of more than 17,500 molecules. After internal validation by five-fold cross-validation and data randomization, the predictive power of the most interesting multiple linear model, based on two GB/SA parameters solely, was tested on two different external sets of molecules. On the Martel druglike test set, the predictive power of the best model (N = 706, r = 0.64, MAE = 1.18, and RMSE = 1.40) is similar to six well-established empirical methods. On the 17-drug test set, our model outperformed all compared empirical methodologies (N = 17, r = 0.94, MAE = 0.38, and RMSE = 0.52). The physical basis of our original GB/SA approach together with its predictive capacity, computational efficiency (1 to 2 s per molecule), and tridimensional molecular graphics capability lay the foundations for a promising predictor, the implicit log P method (iLOGP), to complement the portfolio of drug design tools developed and provided by the SIB Swiss Institute of Bioinformatics.

Prediction of 18-month survival in patients with primary myelodysplastic syndrome. A regression model and scoring system based on the combination of chromosome findings and the Bournemouth score.

The predictive potential of six selected factors was assessed in 72 patients with primary myelodysplastic syndrome using univariate and multivariate logistic regression analysis of survival at 18 months. Factors were age (above median of 69 years), dysplastic features in the three myeloid bone marrow cell lineages, presence of chromosome defects, all metaphases abnormal, double or complex chromosome defects (C23), and a Bournemouth score of 2, 3, or 4 (B234). In the multivariate approach, B234 and C23 proved to be significantly associated with a reduction in the survival probability. The similarity of the regression coefficients associated with these two factors means that they have about the same weight. Consequently, the model was simplified by counting the number of factors (0, 1, or 2) present in each patient, thus generating a scoring system called the Lausanne-Bournemouth score (LB score). The LB score combines the well-recognized and easy-to-use Bournemouth score (B score) with the chromosome defect complexity, C23 constituting an additional indicator of patient outcome. The predicted risk of death within 18 months calculated from the model is as follows: 7.1% (confidence interval: 1.7-24.8) for patients with an LB score of 0, 60.1% (44.7-73.8) for an LB score of 1, and 96.8% (84.5-99.4) for an LB score of 2. The scoring system presented here has several interesting features. The LB score may improve the predictive value of the B score, as it is able to recognize two prognostic groups in the intermediate risk category of patients with B scores of 2 or 3. It has also the ability to identify two distinct prognostic subclasses among RAEB and possibly CMML patients. In addition to its above-described usefulness in the prognostic evaluation, the LB score may bring new insights into the understanding of evolution patterns in MDS. We used the combination of the B score and chromosome complexity to define four classes which may be considered four possible states of myelodysplasia and which describe two distinct evolutional pathways.

Verification at the protein level of the PIF4-mediated external coincidence model for the temperature-adaptive photoperiodic control of plant growth in Arabidopsis thaliana.

Plant circadian clock controls a wide variety of physiological and developmental events, which include the short-days (SDs)-specific promotion of the elongation of hypocotyls during de-etiolation and also the elongation of petioles during vegetative growth. In A. thaliana, the PIF4 gene encoding a phytochrome-interacting basic helix-loop-helix (bHLH) transcription factor plays crucial roles in this photoperiodic control of plant growth. According to the proposed external coincidence model, the PIF4 gene is transcribed precociously at the end of night specifically in SDs, under which conditions the protein product is stably accumulated, while PIF4 is expressed exclusively during the daytime in long days (LDs), under which conditions the protein product is degraded by the light-activated phyB and also the residual proteins are inactivated by the DELLA family of proteins. A number of previous reports provided solid evidence to support this coincidence model mainly at the transcriptional level of the PIF 4 and PIF4-traget genes. Nevertheless, the diurnal oscillation profiles of PIF4 proteins, which were postulated to be dependent on photoperiod and ambient temperature, have not yet been demonstrated. Here we present such crucial evidence on PIF4 protein level to further support the external coincidence model underlying the temperature-adaptive photoperiodic control of plant growth in A. thaliana.

The functions of DNA methylation by CcrM in Caulobacter crescentus: a global approach.

DNA methylation is involved in a diversity of processes in bacteria, including maintenance of genome integrity and regulation of gene expression. Here, using Caulobacter crescentus as a model, we exploit genome-wide experimental methods to uncover the functions of CcrM, a DNA methyltransferase conserved in most Alphaproteobacteria. Using single molecule sequencing, we provide evidence that most CcrM target motifs (GANTC) switch from a fully methylated to a hemi-methylated state when they are replicated, and back to a fully methylated state at the onset of cell division. We show that DNA methylation by CcrM is not required for the control of the initiation of chromosome replication or for DNA mismatch repair. By contrast, our transcriptome analysis shows that >10% of the genes are misexpressed in cells lacking or constitutively over-expressing CcrM. Strikingly, GANTC methylation is needed for the efficient transcription of dozens of genes that are essential for cell cycle progression, in particular for DNA metabolism and cell division. Many of them are controlled by promoters methylated by CcrM and co-regulated by other global cell cycle regulators, demonstrating an extensive cross talk between DNA methylation and the complex regulatory network that controls the cell cycle of C. crescentus and, presumably, of many other Alphaproteobacteria.

hMMP9 as predictive factor for response and progression free survival in breast cancer patients treated with bevacizumab and pegylated liposomal doxorubicin (PLD)

Background: The anti-angiogenic drug, bevacizumab (Bv), is currently used in the treatment of different malignancies including breast cancer. Many angiogenesis-associated molecules are found in the circulation of cancer patients. Until now, there are no prognostic or predictive factors identified in breast cancer patients treated with Bv. We present here the first results of the prospective monitoring of 6 angiogenesis-related molecules in the peripheral blood of breast cancer patients treated with a combination of Bv and PLD in the phase II trial, SAKK 24/06. Methods: Patients were treated with PLD (20 mg/m2) and Bv (10 mg/kg) on days 1 and 15 of each 4-week cycle for a maximum of 6 cycles, followed by Bv monotherapy maintenance (10 mg/m2 q2 weeks) until progression or severe toxicity. Plasma and serum samples were collected at baseline, after 2 months of therapy, then every 3 months and at treatment discontinuation. Enzyme-linked immunosorbent assays (Quantikine, R&D Systems and Reliatech) were used to measure the expression levels of human vascular endothelial growth factor (hVEGF), placental growth factor (hPlGF), matrix metalloproteinase 9 (hMMP9) and soluble VEGF receptors hsVEGFR-1, hsVEGFR-2 and hsVEGFR-3. The log-transformed data (to reduce the skewness) for each marker was analyzed using an analysis of variance (ANOVA) model to determine if there was a difference between the mean of the subgroups of interest (where α = 0.05). The untransformed data was also analyzed in the same manner as a "sensitivity" check. Results: 132 blood samples were collected in 41 out of 43 enrolled patients. Baseline levels of the molecules were compared to disease status according to RECIST. There was a statistically significant difference in the mean of the log-transformed levels of hMMP9 between responders [CR+PR] versus the mean in patients with PD (p-value=0.0004, log fold change=0.7536), and between patients with disease control [CR+PR+SD] and those with PD (p-value=<0.0001, log fold change=0.81559), with the log-transformed level of hMMP9 being higher for the responder group. The mean of the log-transformed levels of hsVEGFR-1 was statistically significantly different between patients with disease control [CR+PR+SD] and those with PD (p-value=0.0068, log fold change=-0.6089), where the log-transformed level of hsVEGFR-1 was lower for the responder group. The log-transformed level of hMMP9 at baseline was identified as a significant prognostic factor in terms of progression free survival (PFS): p-value=0.0417, hazard ratio (HR)=0.574 with a corresponding 95% confidence interval (0.336 - 0.979)). No strong correlation was shown either between the log-transformed levels of hsVEGF, hPlGF, hsVEGFR-2 or hsVEGFR-3 and clinical response or the occurrence of severe toxicity, or between the levels of the different molecules. Conclusions: Our results suggest that baseline plasma level of the matrix metalloproteinase, hMMP9, could predict tumor response and PFS in patients treated with a combination of Bv and PLD. These data justify further investigation in breast cancer patients treated with anti-angiogenic therapy.

Pathological complete response after neoadjuvant chemotherapy is an independent predictive factor irrespective of simplified breast cancer intrinsic subtypes: a landmark and two-step approach analyses from the EORTC 10994/BIG 1-00 phase III trial.

BACKGROUND: Pathological complete response (pCR) following chemotherapy is strongly associated with both breast cancer subtype and long-term survival. Within a phase III neoadjuvant chemotherapy trial, we sought to determine whether the prognostic implications of pCR, TP53 status and treatment arm (taxane versus non-taxane) differed between intrinsic subtypes. PATIENTS AND METHODS: Patients were randomized to receive either six cycles of anthracycline-based chemotherapy or three cycles of docetaxel then three cycles of eprirubicin/docetaxel (T-ET). pCR was defined as no evidence of residual invasive cancer (or very few scattered tumour cells) in primary tumour and lymph nodes. We used a simplified intrinsic subtypes classification, as suggested by the 2011 St Gallen consensus. Interactions between pCR, TP53 status, treatment arm and intrinsic subtype on event-free survival (EFS), distant metastasis-free survival (DMFS) and overall survival (OS) were studied using a landmark and a two-step approach multivariate analyses. RESULTS: Sufficient data for pCR analyses were available in 1212 (65%) of 1856 patients randomized. pCR occurred in 222 of 1212 (18%) patients: 37 of 496 (7.5%) luminal A, 22 of 147 (15%) luminal B/HER2 negative, 51 of 230 (22%) luminal B/HER2 positive, 43 of 118 (36%) HER2 positive/non-luminal, 69 of 221(31%) triple negative (TN). The prognostic effect of pCR on EFS did not differ between subtypes and was an independent predictor for better EFS [hazard ratio (HR) = 0.40, P < 0.001 in favour of pCR], DMFS (HR = 0.32, P < 0.001) and OS (HR = 0.32, P < 0.001). Chemotherapy arm was an independent predictor only for EFS (HR = 0.73, P = 0.004 in favour of T-ET). The interaction between TP53, intrinsic subtypes and survival outcomes only approached statistical significance for EFS (P = 0.1). CONCLUSIONS: pCR is an independent predictor of favourable clinical outcomes in all molecular subtypes in a two-step multivariate analysis. CLINICALTRIALSGOV: EORTC 10994/BIG 1-00 Trial registration number NCT00017095.

Amélioration de l'hyperglycémie hospitalière: création d'un nouveau modèle de gouvernance [Improvement of hospital hyperglycemia: creation of a new governance model?].

The report of significant decrease of the inpatient hospital mortality and morbidity with an efficient insulin therapy has demonstrated the need of a good glycaemic control for patients hospitalised in acute care. However, one is faced with numerous difficulties in the hospital management of patients with hyperglycaemia, errors often occur when prescribing insulin, and the management skills are insufficient. Our goal is to change the medical and nursing practices to evolve towards an efficient and safe management of the hospitalised patient. The model we lay out in this article is based upon observation of the therapeutic support of patients with a chronic condition, whilst using a systemic management approach.

Defective tumor necrosis factor-alpha-dependent control of astrocyte glutamate release in a transgenic mouse model of Alzheimer disease.

The cytokine tumor necrosis factor-alpha (TNFalpha) induces Ca2+-dependent glutamate release from astrocytes via the downstream action of prostaglandin (PG) E2. By this process, astrocytes may participate in intercellular communication and neuromodulation. Acute inflammation in vitro, induced by adding reactive microglia to astrocyte cultures, enhances TNFalpha production and amplifies glutamate release, switching the pathway into a neurodamaging cascade (Bezzi, P., Domercq, M., Brambilla, L., Galli, R., Schols, D., De Clercq, E., Vescovi, A., Bagetta, G., Kollias, G., Meldolesi, J., and Volterra, A. (2001) Nat. Neurosci. 4, 702-710). Because glial inflammation is a component of Alzheimer disease (AD) and TNFalpha is overexpressed in AD brains, we investigated possible alterations of the cytokine-dependent pathway in PDAPP mice, a transgenic model of AD. Glutamate release was measured in acute hippocampal and cerebellar slices from mice at early (4-month-old) and late (12-month-old) disease stages in comparison with age-matched controls. Surprisingly, TNFalpha-evoked glutamate release, normal in 4-month-old PDAPP mice, was dramatically reduced in the hippocampus of 12-month-old animals. This defect correlated with the presence of numerous beta-amyloid deposits and hypertrophic astrocytes. In contrast, release was normal in cerebellum, a region devoid of beta-amyloid deposition and astrocytosis. The Ca2+-dependent process by which TNFalpha evokes glutamate release in acute slices is distinct from synaptic release and displays properties identical to those observed in cultured astrocytes, notably PG dependence. However, prostaglandin E2 induced normal glutamate release responses in 12-month-old PDAPP mice, suggesting that the pathology-associated defect involves the TNFalpha-dependent control of secretion rather than the secretory process itself. Reduced expression of DENN/MADD, a mediator of TNFalpha-PG coupling, might account for the defect. Alteration of this neuromodulatory astrocytic pathway is described here for the first time in relation to Alzheimer disease.

A model of cellular dosimetry for macroscopic tumors in radiopharmaceutical therapy.

PURPOSE: In the radiopharmaceutical therapy approach to the fight against cancer, in particular when it comes to translating laboratory results to the clinical setting, modeling has served as an invaluable tool for guidance and for understanding the processes operating at the cellular level and how these relate to macroscopic observables. Tumor control probability (TCP) is the dosimetric end point quantity of choice which relates to experimental and clinical data: it requires knowledge of individual cellular absorbed doses since it depends on the assessment of the treatment's ability to kill each and every cell. Macroscopic tumors, seen in both clinical and experimental studies, contain too many cells to be modeled individually in Monte Carlo simulation; yet, in particular for low ratios of decays to cells, a cell-based model that does not smooth away statistical considerations associated with low activity is a necessity. The authors present here an adaptation of the simple sphere-based model from which cellular level dosimetry for macroscopic tumors and their end point quantities, such as TCP, may be extrapolated more reliably. METHODS: Ten homogenous spheres representing tumors of different sizes were constructed in GEANT4. The radionuclide 131I was randomly allowed to decay for each model size and for seven different ratios of number of decays to number of cells, N(r): 1000, 500, 200, 100, 50, 20, and 10 decays per cell. The deposited energy was collected in radial bins and divided by the bin mass to obtain the average bin absorbed dose. To simulate a cellular model, the number of cells present in each bin was calculated and an absorbed dose attributed to each cell equal to the bin average absorbed dose with a randomly determined adjustment based on a Gaussian probability distribution with a width equal to the statistical uncertainty consistent with the ratio of decays to cells, i.e., equal to Nr-1/2. From dose volume histograms the surviving fraction of cells, equivalent uniform dose (EUD), and TCP for the different scenarios were calculated. Comparably sized spherical models containing individual spherical cells (15 microm diameter) in hexagonal lattices were constructed, and Monte Carlo simulations were executed for all the same previous scenarios. The dosimetric quantities were calculated and compared to the adjusted simple sphere model results. The model was then applied to the Bortezomib-induced enzyme-targeted radiotherapy (BETR) strategy of targeting Epstein-Barr virus (EBV)-expressing cancers. RESULTS: The TCP values were comparable to within 2% between the adjusted simple sphere and full cellular models. Additionally, models were generated for a nonuniform distribution of activity, and results were compared between the adjusted spherical and cellular models with similar comparability. The TCP values from the experimental macroscopic tumor results were consistent with the experimental observations for BETR-treated 1 g EBV-expressing lymphoma tumors in mice. CONCLUSIONS: The adjusted spherical model presented here provides more accurate TCP values than simple spheres, on par with full cellular Monte Carlo simulations while maintaining the simplicity of the simple sphere model. This model provides a basis for complementing and understanding laboratory and clinical results pertaining to radiopharmaceutical therapy.

Selective in vivo visualization of immune-cell infiltration in a mouse model of autoimmune myocarditis by fluorine-19 cardiac magnetic resonance.

BACKGROUND: The goal of this study was to characterize the performance of fluorine-19 ((19)F) cardiac magnetic resonance (CMR) for the specific detection of inflammatory cells in a mouse model of myocarditis. Intravenously administered perfluorocarbons are taken up by infiltrating inflammatory cells and can be detected by (19)F-CMR. (19)F-labeled cells should, therefore, generate an exclusive signal at the inflamed regions within the myocardium. METHODS AND RESULTS: Experimental autoimmune myocarditis was induced in BALB/c mice. After intravenous injection of 2×200 µL of a perfluorocarbon on day 19 and 20 (n=9) after immunization, in vivo (19)F-CMR was performed at the peak of myocardial inflammation (day 21). In 5 additional animals, perfluorocarbon combined with FITC (fluorescein isothiocyanate) was administered for postmortem immunofluorescence and flow-cytometry analyses. Control experiments were performed in 9 animals. In vivo (19)F-CMR detected myocardial inflammation in all experimental autoimmune myocarditis-positive animals. Its resolution was sufficient to identify even small inflammatory foci, that is, at the surface of the right ventricle. Postmortem immunohistochemistry and flow cytometry confirmed the presence of perfluorocarbon in macrophages, dendritic cells, and granulocytes, but not in lymphocytes. The myocardial volume of elevated (19)F signal (rs=0.96; P<0.001), the (19)F signal-to-noise ratio (rs=0.92; P<0.001), and the (19)F signal integral (rs=0.96; P<0.001) at day 21 correlated with the histological myocarditis severity score. CONCLUSIONS: In vivo (19)F-CMR was successfully used to visualize the inflammation specifically and robustly in experimental autoimmune myocarditis, and thus allowed for an unprecedented insight into the involvement of inflammatory cells in the disease process.