Phase I and phase II xenobiotic reactions and metabolism of the food-borne carcinogen 2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline in aggregating liver cell cultures.

Aggregating fetal liver cell cultures were tested for their ability to metabolize xenobiotics using ethoxycoumarin-O-deethylase (ECOD), as marker of phase I metabolism, and glutathione S-transferase (GST), as marker for phase II reactions. Significant basal activities, stable over 14 days in culture were measured for both ECOD and GST activities. The prototype cytochrome P450 inducers, 3-methylcholanthrene (3-MC) and phenobarbital (PB), increased ECOD and GST activities reaching an optimum 7 days after culturing, followed by a decline in activity. This decline was partially prevented by 1% dimethyl sulfoxide (DMSO) added chronically to the culture medium. DMSO was also found to induce ECOD activity and to a lesser extent GST activity. Furthermore, it potentiated in a dose-dependent manner the induction of ECOD by PB. The food-borne carcinogen 2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline (MeIQx) is metabolically transformed through a number of pathways in vivo. It was therefore used to examine the metabolic capacity in fetal and adult liver cell aggregates. Metabolism of MeIQx was mainly through N2-conjugation, resulting in formation of the N2-glucuronide and sulfamate conjugates for non-induced fetal liver cells. These metabolites were also found in large amounts in non-induced adult liver cells. Low levels of cytochrome P450-mediated ring-hydroxylated metabolites were detected in both non-induced fetal and adult liver cells. After induction with arochlor (PCB) or 3-MC, the major pathway was ring-hydroxylation (cytochrome P450 dependent), followed by conjugation to beta-glucuronic or sulfuric acid. The presence of the glucuronide conjugate of N-hydroxy-MeIQx, a mutagenic metabolite, suggested an induction of P450 CYP1A2. The metabolism of MeIQx by liver cell aggregates is very similar to that observed in vivo and suggests that aggregating liver cell cultures are a useful model for in vitro metabolic studies in toxicology.

NEUROPATHOPHYSIOLOGICAL MECHANISMS IN GLUTARIC ACIDURIA TYPE I: 3-HYDROXYGLUTARIC ACID LEADS TO AMMONIA INCREASE AND NON-APOPTOTIC CELL DEATH

A 3D in vitro model of rat organotypic brain cell cultures in aggregates was used to investigate neurotoxicity mechanisms in glutaric aciduria type I (GA-I). 1 mM glutarate (GA) or 3-hydroxyglutarate (3OHGA) were repeatedly added to the culture media at two different time points. In cultures treated with 3OHGA, we observed an increase in lactate in the medium, pointing to a possible inhibition of Krebs cycle and respiratory chain. We further observed that 3OHGA and to a lesser extend GA induced an increase in ammonia production with concomitant decrease of glutamine concentrations, which may suggest an inhibition of the astrocytic enzyme glutamine synthetase. These previously unreported findings may uncover a pathogenic mechanism in this disease which has deleterious effects on early stages of brain development. By immunohistochemistry we showed that 3OHGA increased non-apoptotic cell death. On the cellular level, 3OHGA and to a lesser extend GA led to cell swelling and loss of astrocytic fibers whereas a loss of oligodendrocytes was only observed for 3OHGA. We conclude that 3OHGAwas the most toxic metabolite in our model for GA-I. 3OHGA induced deleterious effects on glial cells, an increase of ammonia production, and resulted in accentuated cell death of non-apoptotic origin.

FastEpistasis: a high performance computing solution for quantitative trait epistasis.

Motivation: Genome-wide association studies have become widely used tools to study effects of genetic variants on complex diseases. While it is of great interest to extend existing analysis methods by considering interaction effects between pairs of loci, the large number of possible tests presents a significant computational challenge. The number of computations is further multiplied in the study of gene expression quantitative trait mapping, in which tests are performed for thousands of gene phenotypes simultaneously. Results: We present FastEpistasis, an efficient parallel solution extending the PLINK epistasis module, designed to test for epistasis effects when analyzing continuous phenotypes. Our results show that the algorithm scales with the number of processors and offers a reduction in computation time when several phenotypes are analyzed simultaneously. FastEpistasis is capable of testing the association of a continuous trait with all single nucleotide polymorphism ( SNP) pairs from 500 000 SNPs, totaling 125 billion tests, in a population of 5000 individuals in 29, 4 or 0.5 days using 8, 64 or 512 processors.

Ammonium accumulation and cell death in a rat 3D brain cell model of glutaric aciduria type I.

Glutaric aciduria type I (glutaryl-CoA dehydrogenase deficiency) is an inborn error of metabolism that usually manifests in infancy by an acute encephalopathic crisis and often results in permanent motor handicap. Biochemical hallmarks of this disease are elevated levels of glutarate and 3-hydroxyglutarate in blood and urine. The neuropathology of this disease is still poorly understood, as low lysine diet and carnitine supplementation do not always prevent brain damage, even in early-treated patients. We used a 3D in vitro model of rat organotypic brain cell cultures in aggregates to mimic glutaric aciduria type I by repeated administration of 1 mM glutarate or 3-hydroxyglutarate at two time points representing different developmental stages. Both metabolites were deleterious for the developing brain cells, with 3-hydroxyglutarate being the most toxic metabolite in our model. Astrocytes were the cells most strongly affected by metabolite exposure. In culture medium, we observed an up to 11-fold increase of ammonium in the culture medium with a concomitant decrease of glutamine. We further observed an increase in lactate and a concomitant decrease in glucose. Exposure to 3-hydroxyglutarate led to a significantly increased cell death rate. Thus, we propose a three step model for brain damage in glutaric aciduria type I: (i) 3-OHGA causes the death of astrocytes, (ii) deficiency of the astrocytic enzyme glutamine synthetase leads to intracerebral ammonium accumulation, and (iii) high ammonium triggers secondary death of other brain cells. These unexpected findings need to be further investigated and verified in vivo. They suggest that intracerebral ammonium accumulation might be an important target for the development of more effective treatment strategies to prevent brain damage in patients with glutaric aciduria type I.

Dynamic simulation of regulatory networks using SQUAD.

BACKGROUND: The ambition of most molecular biologists is the understanding of the intricate network of molecular interactions that control biological systems. As scientists uncover the components and the connectivity of these networks, it becomes possible to study their dynamical behavior as a whole and discover what is the specific role of each of their components. Since the behavior of a network is by no means intuitive, it becomes necessary to use computational models to understand its behavior and to be able to make predictions about it. Unfortunately, most current computational models describe small networks due to the scarcity of kinetic data available. To overcome this problem, we previously published a methodology to convert a signaling network into a dynamical system, even in the total absence of kinetic information. In this paper we present a software implementation of such methodology. RESULTS: We developed SQUAD, a software for the dynamic simulation of signaling networks using the standardized qualitative dynamical systems approach. SQUAD converts the network into a discrete dynamical system, and it uses a binary decision diagram algorithm to identify all the steady states of the system. Then, the software creates a continuous dynamical system and localizes its steady states which are located near the steady states of the discrete system. The software permits to make simulations on the continuous system, allowing for the modification of several parameters. Importantly, SQUAD includes a framework for perturbing networks in a manner similar to what is performed in experimental laboratory protocols, for example by activating receptors or knocking out molecular components. Using this software we have been able to successfully reproduce the behavior of the regulatory network implicated in T-helper cell differentiation. CONCLUSION: The simulation of regulatory networks aims at predicting the behavior of a whole system when subject to stimuli, such as drugs, or determine the role of specific components within the network. The predictions can then be used to interpret and/or drive laboratory experiments. SQUAD provides a user-friendly graphical interface, accessible to both computational and experimental biologists for the fast qualitative simulation of large regulatory networks for which kinetic data is not necessarily available.

The InterPro Database, 2003 brings increased coverage and new features.

InterPro, an integrated documentation resource of protein families, domains and functional sites, was created in 1999 as a means of amalgamating the major protein signature databases into one comprehensive resource. PROSITE, Pfam, PRINTS, ProDom, SMART and TIGRFAMs have been manually integrated and curated and are available in InterPro for text- and sequence-based searching. The results are provided in a single format that rationalises the results that would be obtained by searching the member databases individually. The latest release of InterPro contains 5629 entries describing 4280 families, 1239 domains, 95 repeats and 15 post-translational modifications. Currently, the combined signatures in InterPro cover more than 74% of all proteins in SWISS-PROT and TrEMBL, an increase of nearly 15% since the inception of InterPro. New features of the database include improved searching capabilities and enhanced graphical user interfaces for visualisation of the data. The database is available via a webserver (http://www.ebi.ac.uk/interpro) and anonymous FTP (ftp://ftp.ebi.ac.uk/pub/databases/interpro).

Extent of hypoattenuation on CT angiography source images in basilar artery occlusion: prognostic value in the Basilar Artery International Cooperation Study.

BACKGROUND AND PURPOSE: The posterior circulation Acute Stroke Prognosis Early CT Score (pc-ASPECTS) quantifies the extent of early ischemic changes in the posterior circulation with a 10-point grading system. We hypothesized that pc-ASPECTS applied to CT angiography source images predicts functional outcome of patients in the Basilar Artery International Cooperation Study (BASICS). METHODS: BASICS was a prospective, observational registry of consecutive patients with acute symptomatic basilar artery occlusion. Functional outcome was assessed at 1 month. We applied pc-ASPECTS to CT angiography source images of patients with CT angiography for confirmation of basilar artery occlusion. We calculated unadjusted and adjusted risk ratios (RRs) of pc-ASPECTS dichotomized at ≥8 versus <8. Primary outcome measure was favorable outcome (modified Rankin Scale scores 0-3). Secondary outcome measures were mortality and functional independence (modified Rankin Scale scores 0-2). RESULTS: Of 158 patients included, 78 patients had a CT angiography source images pc-ASPECTS≥8. Patients with a pc-ASPECTS≥8 more often had a favorable outcome than patients with a pc-ASPECTS<8 (crude RR, 1.7; 95% CI, 0.98-3.0). After adjustment for age, baseline National Institutes of Health Stroke Scale score, and thrombolysis, pc-ASPECTS≥8 was not related to favorable outcome (RR, 1.3; 95% CI, 0.8-2.2), but it was related to reduced mortality (RR, 0.7; 95% CI, 0.5-0.98) and functional independence (RR, 2.0; 95% CI, 1.1-3.8). In post hoc analysis, pc-ASPECTS dichotomized at ≥6 versus <6 predicted a favorable outcome (adjusted RR, 3.1; 95% CI, 1.2-7.5). CONCLUSIONS: pc-ASPECTS on CT angiography source images independently predicted death and functional independence at 1 month in the CT angiography subgroup of patients in the BASICS registry.

Computer-assisted cup placement techniques in total hip arthroplasty improve accuracy of placement.

Malposition of the acetabular component during hip arthroplasty increases the occurrence of impingement, reduces range of motion, and increases the risk of dislocation and long-term wear. To prevent malpositioned hip implants, an increasing number of computer-assisted orthopaedic systems have been described, but their accuracy is not well established. The purpose of this study was to determine the reproducibility and accuracy of conventional versus computer-assisted techniques for positioning the acetabular component in total hip arthroplasty. Using a lateral approach, 150 cups were placed by 10 surgeons in 10 identical plastic pelvis models (freehand, with a mechanical guide, using computer assistance). Conditions for cup implantations were made to mimic the operating room situation. Preoperative planning was done from a computed tomography scan. The accuracy of cup abduction and anteversion was assessed with an electromagnetic system. Freehand placement revealed a mean accuracy of cup anteversion and abduction of 10 degrees and 3.5 degrees, respectively (maximum error, 35 degrees). With the cup positioner, these angles measured 8 degrees and 4 degrees (maximum error, 29.8 degrees), respectively, and using computer assistance, 1.5 degrees and 2.5 degrees degrees (maximum error, 8 degrees), respectively. Computer-assisted cup placement was an accurate and reproducible technique for total hip arthroplasty. It was more accurate than traditional methods of cup positioning.

3D Automated Lung Nodule Segmentation in HRCT

A fully-automated 3D image analysis method is proposed to segment lung nodules in HRCT. A specific gray-level mathematical morphology operator, the SMDC-connection cost, acting in the 3D space of the thorax volume is defined in order to discriminate lung nodules from other dense (vascular) structures. Applied to clinical data concerning patients with pulmonary carcinoma, the proposed method detects isolated, juxtavascular and peripheral nodules with sizes ranging from 2 to 20 mm diameter. The segmentation accuracy was objectively evaluated on real and simulated nodules. The method showed a sensitivity and a specificity ranging from 85% to 97% and from 90% to 98%, respectively.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

Should Koos grade I vestibular schwannomas be treated early with gamma knife surgery? A subgroup analysis in a series of 190 consecutive patients.

INTRODUCTION: Gamma knife surgery (GKS) for vestibular schwannomas (VS) has a long-term clinical and scientific track record. After a period of de-escalation of dose prescription, results show a high rate of tumor control with improvement of clinical outcome (less than 1% facial palsy, 50-70% hearing preservation). Régis et al. (J Neurosurg 2013;119 Suppl.:105-11) suggested recently that proactive GKS management in intracanalicular tumors is better than a « wait and see » strategy when hearing is still useful at the time of diagnosis. MATERIALS AND METHODS: Based on these previous findings, we prospectively analyzed 190 vestibular schwannomas (VS), treated with GKS as first intention over a period of 4 years (2010-2014). We concentrated on patient, tumor and dosimetric characteristics. Special attention was given on the dose to the cochlea and its impact in maintaining serviceable hearing. RESULTS: The mean follow-up period was 1.3years (range 0.6-3.6). Preoperative serviceable hearing was present in 63.11% patients. The mean maximal diameter was 15.1mm (range 5-29.5). The size and volume of the tumor corresponded to Koos grade I, II, III and IV in 15.9%, 34.8%, 45.4% and 3.8% of the cases, respectively. The mean target volume was 1.24cm(3) (0.017-7.8). The mean prescription isodose volume was 1.6 cc (0.032-8.5). The mean marginal dose was 12Gy (11-12). The mean maximal dose received by the cochlea in patients with GR class 1 and 2 was 4.1Gy (1.5-7.6). Our preliminary neuroradiological follow-up shows 97% tumor control, with 45% shrinkage. Patients presenting with GR class 1 and class 2 at baseline retained serviceable hearing in 85% of cases. Among the patients with a follow-up of at least one year, those with Koos I tumors had the highest probability to maintain identical level of hearing after GKS. CONCLUSION: Our preliminary data suggest that Koos I patients should be treated early with GKS, before tumor growth and/or hearing deterioration, as they have the highest probability of hearing preservation.

Inferring ultraviolet anatomical exposure patterns while distinguishing the relative contribution of radiation components

Exposure to solar ultraviolet (UV) radiation is the main causative factor for skin cancer. UV exposure depends on environmental and individual factors, but individual exposure data remain scarce. UV irradiance is monitored via different techniques including ground measurements and satellite observations. However it is difficult to translate such observations into human UV exposure or dose because of confounding factors (shape of the exposed surface, shading, behavior, etc.) A collaboration between public health institutions, a meteorological office and an institute specialized in computing techniques developed a model predicting the dose and distribution of UV exposure on the basis of ground irradiation and morphological data. Standard 3D computer graphics techniques were adapted to develop this tool, which estimates solar exposure of a virtual manikin depicted as a triangle mesh surface. The amount of solar energy received by various body locations is computed for direct, diffuse and reflected radiation separately. The radiation components are deduced from corresponding measurements of UV irradiance, and the related UV dose received by each triangle of the virtual manikin is computed accounting for shading by other body parts and eventual protection measures. The model was verified with dosimetric measurements (n=54) in field conditions using a foam manikin as surrogate for an exposed individual. Dosimetric results were compared to the model predictions. The model predicted exposure to solar UV adequately. The symmetric mean absolute percentage error was 13%. Half of the predictions were within 17% range of the measurements. This model allows assessing outdoor occupational and recreational UV exposures, without necessitating time-consuming individual dosimetry, with numerous potential uses in skin cancer prevention and research. Using this tool, we investigated solar UV exposure patterns with respect to the relative contribution of the direct, diffuse and reflected radiation. We assessed exposure doses for various body parts and exposure scenarios of a standing individual (static and dynamic postures). As input, the model used erythemally-weighted ground irradiance data measured in 2009 at Payerne, Switzerland. A year-round daily exposure (8 am to 5 pm) without protection was assumed. For most anatomical sites, mean daily doses were high (typically 6.2-14.6 SED) and exceeded recommended exposure values. Direct exposure was important during specific periods (e.g. midday during summer), but contributed moderately to the annual dose, ranging from 15 to 24% for vertical and horizontal body parts, respectively. Diffuse irradiation explained about 80% of the cumulative annual exposure dose. Acute diffuse exposures were also obtained for cloudy summer days. The importance of diffuse UV radiation should not be underestimated when advocating preventive measures. Messages focused on avoiding acute direct exposures may be of limited efficiency to prevent skin cancers associated with chronic exposure (e.g., squamous cell carcinomas).

Multi-slice three-dimensional myocardial strain tensor quantification using zHARP.

In this article we propose a novel method for calculating cardiac 3-D strain. The method requires the acquisition of myocardial short-axis (SA) slices only and produces the 3-D strain tensor at every point within every pair of slices. Three-dimensional displacement is calculated from SA slices using zHARP which is then used for calculating the local displacement gradient and thus the local strain tensor. There are three main advantages of this method. First, the 3-D strain tensor is calculated for every pixel without interpolation; this is unprecedented in cardiac MR imaging. Second, this method is fast, in part because there is no need to acquire long-axis (LA) slices. Third, the method is accurate because the 3-D displacement components are acquired simultaneously and therefore reduces motion artifacts without the need for registration. This article presents the theory of computing 3-D strain from two slices using zHARP, the imaging protocol, and both phantom and in-vivo validation.

The joint influence of gender and amount of smoking on weight gain one year after smoking cessation.

Weight gain is often associated with smoking cessation and may discourage smokers from quitting. This study estimated the weight gained one year after smoking cessation and examined the risk factors associated with weight gain in order to identify socio-demographic groups at higher risk of increased weight after quitting. We analyzed data from 750 adults in two randomized controlled studies that included smokers motivated to quit and found a gradient in weight gain according to the actual duration of abstinence during follow-up. Subjects who were abstinent for at least 40 weeks gained 4.6 kg (SD = 3.8) on average, compared to 1.2 kg (SD = 2.6) for those who were abstinent less than 20 weeks during the 1-year follow-up. Considering the duration of abstinence as an exposure variable, we found an age effect and a significant interaction between sex and the amount of smoking before quitting: younger subjects gained more weight than older subjects; among light smokers, men gained more weight on average than women one year after quitting, while the opposite was observed among heavy smokers. Young women smoking heavily at baseline had the highest risk of weight gain after quitting.