Monte-Carlo dosimetry on a realistic cell monolayer geometry exposed to alpha particles

The energy and specific energy absorbed in the main cell compartments (nucleus and cytoplasm) in typical radiobiology experiments are usually estimated by calculations as they are not accessible for a direct measurement. In most of the work, the cell geometry is modelled using the combination of simple mathematical volumes. We propose a method based on high resolution confocal imaging and ion beam analysis (IBA) in order to import realistic cell nuclei geometries in Monte-Carlo simulations and thus take into account the variety of different geometries encountered in a typical cell population. Seventy-six cell nuclei have been imaged using confocal microscopy and their chemical composition has been measured using IBA. A cellular phantom was created from these data using the ImageJ image analysis software and imported in the Geant4 Monte-Carlo simulation toolkit. Total energy and specific energy distributions in the 76 cell nuclei have been calculated for two types of irradiation protocols: a 3 MeV alpha particle microbeam used for targeted irradiation and a 239Pu alpha source used for large angle random irradiation. Qualitative images of the energy deposited along the particle tracks have been produced and show good agreement with images of DNA double strand break signalling proteins obtained experimentally. The methodology presented in this paper provides microdosimetric quantities calculated from realistic cellular volumes. It is based on open-source oriented software that is publicly available.

Colocalization of α-actinin and Synaptopodin in the Pyramidal Cell Axon Initial Segment

The cisternal organelle that resides in the axon initial segment (AIS) of neocortical and hippocampal pyramidal cells is thought to be involved in regulating the Ca(2+) available to maintain AIS scaffolding proteins, thereby preserving normal AIS structure and function. Through immunocytochemistry and correlative light and electron microscopy, we show here that the actin-binding protein ?-actinin is present in the typical cistenal organelle of rodent pyramidal neurons as well as in a large structure in the AIS of a subpopulation of layer V pyramidal cells that we have called the "giant saccular organelle." Indeed, this localization of ?-actinin in the AIS is dependent on the integrity of the actin cytoskeleton. Moreover, in the cisternal organelle of cultured hippocampal neurons, ?-actinin colocalizes extensively with synaptopodin, a protein that interacts with both actin and ?-actinin, and they appear concomitantly during the development of these neurons. Together, these results indicate that ?-actinin and the actin cytoskeleton are important components of the cisternal organelle that are probably required to stabilize the AIS.

Analysis of large-scale digital optical neural networks by Feynman diagrams

A new method to study large scale neural networks is presented in this paper. The basis is the use of Feynman- like diagrams. These diagrams allow the analysis of collective and cooperative phenomena with a similar methodology to the employed in the Many Body Problem. The proposed method is applied to a very simple structure composed by an string of neurons with interaction among them. It is shown that a new behavior appears at the end of the row. This behavior is different to the initial dynamics of a single cell. When a feedback is present, as in the case of the hippocampus, this situation becomes more complex with a whole set of new frequencies, different from the proper frequencies of the individual neurons. Application to an optical neural network is reported.

Large scale multiphysics modeling of neurites

Situado en el límite entre Ingeniería, Informática y Biología, la mecánica computacional de las neuronas aparece como un nuevo campo interdisciplinar que potencialmente puede ser capaz de abordar problemas clínicos desde una perspectiva diferente. Este campo es multiescala por naturaleza, yendo desde la nanoescala (como, por ejemplo, los dímeros de tubulina) a la macroescala (como, por ejemplo, el tejido cerebral), y tiene como objetivo abordar problemas que son complejos, y algunas veces imposibles, de estudiar con medios experimentales. La modelización computacional ha sido ampliamente empleada en aplicaciones Neurocientíficas tan diversas como el crecimiento neuronal o la propagación de los potenciales de acción compuestos. Sin embargo, en la mayoría de los enfoques de modelización hechos hasta ahora, la interacción entre la célula y el medio/estímulo que la rodea ha sido muy poco explorada. A pesar de la tremenda importancia de esa relación en algunos desafíos médicos—como, por ejemplo, lesiones traumáticas en el cerebro, cáncer, la enfermedad del Alzheimer—un puente que relacione las propiedades electrofisiológicas-químicas y mecánicas desde la escala molecular al nivel celular todavía no existe. Con ese objetivo, esta investigación propone un marco computacional multiescala particularizado para dos escenarios respresentativos: el crecimiento del axón y el acomplamiento electrofisiológicomecánico de las neuritas. En el primer caso, se explora la relación entre los constituyentes moleculares del axón durante su crecimiento y sus propiedades mecánicas resultantes, mientras que en el último, un estímulo mecánico provoca deficiencias funcionales a nivel celular como consecuencia de sus alteraciones electrofisiológicas-químicas. La modelización computacional empleada en este trabajo es el método de las diferencias finitas, y es implementada en un nuevo programa llamado Neurite. Aunque el método de los elementos finitos es también explorado en parte de esta investigación, el método de las diferencias finitas tiene la flexibilidad y versatilidad necesaria para implementar mode los biológicos, así como la simplicidad matemática para extenderlos a simulaciones a gran escala con un coste computacional bajo. Centrándose primero en el efecto de las propiedades electrofisiológicas-químicas sobre las propiedades mecánicas, una versión adaptada de Neurite es desarrollada para simular la polimerización de los microtúbulos en el crecimiento del axón y proporcionar las propiedades mecánicas como función de la ocupación de los microtúbulos. Después de calibrar el modelo de crecimiento del axón frente a resultados experimentales disponibles en la literatura, las características mecánicas pueden ser evaluadas durante la simulación. Las propiedades mecánicas del axón muestran variaciones dramáticas en la punta de éste, donde el cono de crecimiento soporta las señales químicas y mecánicas. Bansándose en el conocimiento ganado con el modelo de diferencias finitas, y con el objetivo de ir de 1D a 3D, este esquema preliminar pero de una naturaleza innovadora allana el camino a futuros estudios con el método de los elementos finitos. Centrándose finalmente en el efecto de las propiedades mecánicas sobre las propiedades electrofisiológicas- químicas, Neurite es empleado para relacionar las cargas mecánicas macroscópicas con las deformaciones y velocidades de deformación a escala microscópica, y simular la propagación de la señal eléctrica en las neuritas bajo carga mecánica. Las simulaciones fueron calibradas con resultados experimentales publicados en la literatura, proporcionando, por tanto, un modelo capaz de predecir las alteraciones de las funciones electrofisiológicas neuronales bajo cargas externas dañinas, y uniendo lesiones mecánicas con las correspondientes deficiencias funcionales. Para abordar simulaciones a gran escala, aunque otras arquitecturas avanzadas basadas en muchos núcleos integrados (MICs) fueron consideradas, los solvers explícito e implícito se implementaron en unidades de procesamiento central (CPU) y unidades de procesamiento gráfico (GPUs). Estudios de escalabilidad fueron llevados acabo para ambas implementaciones mostrando resultados prometedores para casos de simulaciones extremadamente grandes con GPUs. Esta tesis abre la vía para futuros modelos mecánicos con el objetivo de unir las propiedades electrofisiológicas-químicas con las propiedades mecánicas. El objetivo general es mejorar el conocimiento de las comunidades médicas y de bioingeniería sobre la mecánica de las neuronas y las deficiencias funcionales que aparecen de los daños producidos por traumatismos mecánicos, como lesiones traumáticas en el cerebro, o enfermedades neurodegenerativas como la enfermedad del Alzheimer. ABSTRACT Sitting at the interface between Engineering, Computer Science and Biology, Computational Neuron Mechanics appears as a new interdisciplinary field potentially able to tackle clinical problems from a new perspective. This field is multiscale by nature, ranging from the nanoscale (e.g., tubulin dimers) to the macroscale (e.g., brain tissue), and aims at tackling problems that are complex, and sometime impossible, to study through experimental means. Computational modeling has been widely used in different Neuroscience applications as diverse as neuronal growth or compound action potential propagation. However, in the majority of the modeling approaches done in this field to date, the interactions between the cell and its surrounding media/stimulus have been rarely explored. Despite of the tremendous importance of such relationship in several medical challenges—e.g., traumatic brain injury (TBI), cancer, Alzheimer’s disease (AD)—a bridge between electrophysiological-chemical and mechanical properties of neurons from the molecular scale to the cell level is still lacking. To this end, this research proposes a multiscale computational framework particularized for two representative scenarios: axon growth and electrophysiological-mechanical coupling of neurites. In the former case, the relation between the molecular constituents of the axon during its growth and its resulting mechanical properties is explored, whereas in the latter, a mechanical stimulus provokes functional deficits at cell level as a consequence of its electrophysiological-chemical alterations. The computational modeling approach chosen in this work is the finite difference method (FDM), and was implemented in a new program called Neurite. Although the finite element method (FEM) is also explored as part of this research, the FDM provides the necessary flexibility and versatility to implement biological models, as well as the mathematical simplicity to extend them to large scale simulations with a low computational cost. Focusing first on the effect of electrophysiological-chemical properties on the mechanical proper ties, an adaptation of Neurite was developed to simulate microtubule polymerization in axonal growth and provide the axon mechanical properties as a function of microtubule occupancy. After calibrating the axon growth model against experimental results available in the literature, the mechanical characteristics can be tracked during the simulation. The axon mechanical properties show dramatic variations at the tip of the axon, where the growth cone supports the chemical and mechanical signaling. Based on the knowledge gained from the FDM scheme, and in order to go from 1D to 3D, this preliminary yet novel scheme paves the road for future studies with FEM. Focusing then on the effect of mechanical properties on the electrophysiological-chemical properties, Neurite was used to relate macroscopic mechanical loading to microscopic strains and strain rates, and simulate the electrical signal propagation along neurites under mechanical loading. The simulations were calibrated against experimental results published in the literature, thus providing a model able to predict the alteration of neuronal electrophysiological function under external damaging load, and linking mechanical injuries to subsequent acute functional deficits. To undertake large scale simulations, although other state-of-the-art architectures based on many integrated cores (MICs) were considered, the explicit and implicit solvers were implemented for central processing units (CPUs) and graphics processing units (GPUs). Scalability studies were done for both implementations showing promising results for extremely large scale simulations with GPUs. This thesis opens the avenue for future mechanical modeling approaches aimed at linking electrophysiological- chemical properties to mechanical properties. Its overarching goal is to enhance the bioengineering and medical communities knowledge on neuronal mechanics and functional deficits arising from damages produced by direct mechanical insults, such as TBI, or neurodegenerative evolving illness, such as AD.

Neurite, a finite difference large scale parallel program for the simulation of the electrical signal propagation in neurites under mechanical loading

With the growing body of research on traumatic brain injury and spinal cord injury, computational neuroscience has recently focused its modeling efforts on neuronal functional deficits following mechanical loading. However, in most of these efforts, cell damage is generally only characterized by purely mechanistic criteria, function of quantities such as stress, strain or their corresponding rates. The modeling of functional deficits in neurites as a consequence of macroscopic mechanical insults has been rarely explored. In particular, a quantitative mechanically based model of electrophysiological impairment in neuronal cells has only very recently been proposed (Jerusalem et al., 2013). In this paper, we present the implementation details of Neurite: the finite difference parallel program used in this reference. Following the application of a macroscopic strain at a given strain rate produced by a mechanical insult, Neurite is able to simulate the resulting neuronal electrical signal propagation, and thus the corresponding functional deficits. The simulation of the coupled mechanical and electrophysiological behaviors requires computational expensive calculations that increase in complexity as the network of the simulated cells grows. The solvers implemented in Neurite-explicit and implicit-were therefore parallelized using graphics processing units in order to reduce the burden of the simulation costs of large scale scenarios. Cable Theory and Hodgkin-Huxley models were implemented to account for the electrophysiological passive and active regions of a neurite, respectively, whereas a coupled mechanical model accounting for the neurite mechanical behavior within its surrounding medium was adopted as a link between lectrophysiology and mechanics (Jerusalem et al., 2013). This paper provides the details of the parallel implementation of Neurite, along with three different application examples: a long myelinated axon, a segmented dendritic tree, and a damaged axon. The capabilities of the program to deal with large scale scenarios, segmented neuronal structures, and functional deficits under mechanical loading are specifically highlighted.

Study of the influence of actin-binding proteins using linear analyses of cell deformability

The actin cytoskeleton plays a key role in the deformability of the cell and in mechanosensing. Here we analyze the contributions of three major actin cross-linking proteins, myosin II, a-actinin and filamin, to cell deformability, by using micropipette aspiration of Dictyostelium cells. We examine the applicability of three simple mechanical models: for small deformation, linear viscoelasticity and drop of liquid with a tense cortex; and for large deformation, a Newtonian viscous fluid. For these models, we have derived linearized equations and we provide a novel, straightforward methodology to analyze the experiments. This methodology allowed us to differentiate the effects of the cross-linking proteins in the different regimes of deformation. Our results confirm some previous observations and suggest important relations between the molecular characteristics of the actin-binding proteins and the cell behavior: the effect of myosin is explained in terms of the relation between the lifetime of the bond to actin and the resistive force; the presence of a-actinin obstructs the deformation of the cytoskeleton, presumably mainly due to the higher molecular stiffness and to the lower dissociation rate constants; and filamin contributes critically to the global connectivity of the network, possibly by rapidly turning over crosslinks during the remodeling of the cytoskeletal network, thanks to the higher rate constants, flexibility and larger size. The results suggest a sophisticated relationship between the expression levels of actinbinding proteins, deformability and mechanosensing.

A logical analysis of T cell activation and anergy

Interaction of the antigen-specific receptor of T lymphocytes with its antigenic ligand can lead either to cell activation or to a state of profound unresponsiveness (anergy). Although subtle changes in the nature of the ligand or of the antigen-presenting cell have been shown to affect the outcome of T cell receptor ligation, the mechanism by which the same receptor can induce alternative cellular responses is not completely understood. A model for explaining both positive (cell proliferation and cytokine production) and negative (anergy induction) signaling of T lymphocytes is described herein. This model relies on the autophosphorylative properties of the tyrosine kinases associated with the T cell receptor. One of its basic assumptions is that the kinase activity of these receptor-associated enzymes remains above background level after ligand removal and is responsible for cellular unresponsiveness. Using a simple Boolean formalism, we show how the timing of the binding and intracellular signal-transduction events can affect the properties of receptor signaling and determine the type of cellular response. The present approach integrates into a common framework a large body of experimental observations and allows specification of conditions leading to cellular activation or to anergy.

Essential role of germinal vesicle material in the meiotic cell cycle of Xenopus oocytes

In almost all animal species, immature oocytes are arrested naturally in the first meiotic prophase, with a large nucleus called the germinal vesicle. A number of previous studies showed that both activation of maturation/M phase-promoting factor (MPF) (assayed by semiquantitative cytological methods) and some other maturational events occur essentially normally in enucleated oocytes from many amphibian species and mice. Hence, for nearly three decades, it has generally been believed that nuclear material is dispensable for MPF activation and the meiotic cell cycle in vertebrate oocytes. Here, we have challenged this view by examining the histone H1 kinase activities and the molecular forms of MPF in experimentally manipulated Xenopus oocytes. We show that oocytes injected with nuclear material undergo much more rapid MPF activation and maturation than uninjected control oocytes. Conversely, enucleated oocytes, unlike nucleated counterparts, undergo only weak MPF activation in meiosis I and no detectable MPF reactivation in meiosis II, the latter accompanying inhibitory tyrosine phosphorylation of cdc2 kinase, the catalytic subunit of MPF. These results argue strongly that nuclear material is indispensable for the meiotic cell cycle, particularly MPF reactivation (or cdc2 tyrosine dephosphorylation) on entry into meiosis II, in Xenopus oocytes. The classical and general view may thus need reconsideration.

Impaired fetal T cell development and perinatal lethality in mice lacking the cAMP response element binding protein

CREB, the cAMP response element binding protein, is a key transcriptional regulator of a large number of genes containing a CRE consensus sequence in their upstream regulatory regions. Mice with a hypomorphic allele of CREB that leads to a loss of the CREBα and Δ isoforms and to an overexpression of the CREBβ isoform are viable. Herein we report the generation of CREB null mice, which have all functional isoforms (CREBα, β, and Δ) inactivated. In contrast to the CREBαΔ mice, CREB null mice are smaller than their littermates and die immediately after birth from respiratory distress. In brain, a strong reduction in the corpus callosum and the anterior commissures is observed. Furthermore, CREB null mice have an impaired fetal T cell development of the αβ lineage, which is not affected in CREBαΔ mice on embryonic day 18.5. Overall thymic cellularity in CREB null mice is severely reduced affecting all developmental stages of the αβ T cell lineage. In contrast γδ T cell differentiation is normal in CREB mutant mice.

T cell positive selection by a high density, low affinity ligand

Interaction of the αβ T cell receptor (TCR) with major histocompatibility (MHC) molecules occupied with any of a large collection of peptides derived from self proteins is a critical step in driving T cell “positive” selection in the thymus. Interaction with this same pool of self-peptide/MHC ligands deletes T cells with potential self-reactivity. To examine how T cells survive both of these processes to form a self-tolerant mature repertoire, mice were constructed whose entire class II MHC IEk specific repertoire was positively selected on a single peptide covalently attached to the IEk molecule. In these mice T cells were identified that could respond to a variant of the positively selecting peptide bound to IEk. The affinities of the TCRs from these T cells for the positively selecting ligand were extremely low and at least 10-fold less than those for the activating ligand. These results support the theory that positive selection is driven by TCR affinities lower than those involved in T cell deletion or activation and that, if present at high concentration, even very low affinity ligands can positively select.

Protein kinase mutants of human ATR increase sensitivity to UV and ionizing radiation and abrogate cell cycle checkpoint control

In fission yeast, the rad3 gene product plays a critical role in sensing DNA structure defects and activating damage response pathways. A structural homologue of rad3 in humans (ATR) has been identified based on sequence similarity in the protein kinase domain. General information regarding ATR expression, protein kinase activity, and cellular localization is known, but its function in human cells remains undetermined. In the current study, the ATR protein was examined by gel filtration of protein extracts and was found to exist predominantly as part of a large protein complex. A kinase-inactivated form of the ATR gene was prepared by site-directed mutagenesis and was used in transfection experiments to probe the function of this complex. Introduction of this kinase-dead ATR into a normal fibroblast cell line, an ATM-deficient fibroblast line derived from a patient with ataxia–telangiectasia, or a p53 mutant cell line all resulted in significant losses in cell viability. Clones expressing the kinase-dead ATR displayed increased sensitivity to x-rays and UV and a loss of checkpoint control. We conclude that ATR functions as a critical part of a protein complex that mediates responses to ionizing and UV radiation in human cells. These responses include effects on cell viability and cell cycle checkpoint control.

Down-regulation of transmembrane carbonic anhydrases in renal cell carcinoma cell lines by wild-type von Hippel-Lindau transgenes

To discover genes involved in von Hippel-Lindau (VHL)-mediated carcinogenesis, we used renal cell carcinoma cell lines stably transfected with wild-type VHL-expressing transgenes. Large-scale RNA differential display technology applied to these cell lines identified several differentially expressed genes, including an alpha carbonic anhydrase gene, termed CA12. The deduced protein sequence was classified as a one-pass transmembrane CA possessing an apparently intact catalytic domain in the extracellular CA module. Reintroduced wild-type VHL strongly inhibited the overexpression of the CA12 gene in the parental renal cell carcinoma cell lines. Similar results were obtained with CA9, encoding another transmembrane CA with an intact catalytic domain. Although both domains of the VHL protein contribute to regulation of CA12 expression, the elongin binding domain alone could effectively regulate CA9 expression. We mapped CA12 and CA9 loci to chromosome bands 15q22 and 17q21.2 respectively, regions prone to amplification in some human cancers. Additional experiments are needed to define the role of CA IX and CA XII enzymes in the regulation of pH in the extracellular microenvironment and its potential impact on cancer cell growth.

Conditional switching of vascular endothelial growth factor (VEGF) expression in tumors: Induction of endothelial cell shedding and regression of hemangioblastoma-like vessels by VEGF withdrawal

We have recently shown that VEGF functions as a survival factor for newly formed vessels during developmental neovascularization, but is not required for maintenance of mature vessels. Reasoning that expanding tumors contain a significant fraction of newly formed and remodeling vessels, we examined whether abrupt withdrawal of VEGF will result in regression of preformed tumor vessels. Using a tetracycline-regulated VEGF expression system in xenografted C6 glioma cells, we showed that shutting off VEGF production leads to detachment of endothelial cells from the walls of preformed vessels and their subsequent death by apoptosis. Vascular collapse then leads to hemorrhages and extensive tumor necrosis. These results suggest that enforced withdrawal of vascular survival factors can be applied to target preformed tumor vasculature in established tumors. The system was also used to examine phenotypes resulting from over-expression of VEGF. When expression of the transfected VEGF cDNA was continuously “on,” tumors became hyper-vascularized with abnormally large vessels, presumably arising from excessive fusions. Tumors were significantly less necrotic, suggesting that necrosis in these tumors is the result of insufficient angiogenesis.

Intrastriatal injection of an adenoviral vector expressing glial-cell-line-derived neurotrophic factor prevents dopaminergic neuron degeneration and behavioral impairment in a rat model of Parkinson disease

Glial-cell-line-derived neurotrophic factor (GDNF) is a potent neurotrophic factor for adult nigral dopamine neurons in vivo. GDNF has both protective and restorative effects on the nigro-striatal dopaminergic (DA) system in animal models of Parkinson disease. Appropriate administration of this factor is essential for the success of its clinical application. Since it cannot cross the blood–brain barrier, a gene transfer method may be appropriate for delivery of the trophic factor to DA cells. We have constructed a recombinant adenovirus (Ad) encoding GDNF and injected it into rat striatum to make use of its ability to infect neurons and to be retrogradely transported by DA neurons. Ad-GDNF was found to drive production of large amounts of GDNF, as quantified by ELISA. The GDNF produced after gene transfer was biologically active: it increased the survival and differentiation of DA neurons in vitro. To test the efficacy of the Ad-mediated GDNF gene transfer in vivo, we used a progressive lesion model of Parkinson disease. Rats received injections unilaterally into their striatum first of Ad and then 6 days later of 6-hydroxydopamine. We found that mesencephalic nigral dopamine neurons of animals treated with the Ad-GDNF were protected, whereas those of animals treated with the Ad-β-galactosidase were not. This protection was associated with a difference in motor function: amphetamine-induced turning was much lower in animals that received the Ad-GDNF than in the animals that received Ad-β-galactosidase. This finding may have implications for the development of a treatment for Parkinson disease based on the use of neurotrophic factors.

The identification of CD4+ T cell epitopes with dedicated synthetic peptide libraries

For a large number of T cell-mediated immunopathologies, the disease-related antigens are not yet identified. Identification of T cell epitopes is of crucial importance for the development of immune-intervention strategies. We show that CD4+ T cell epitopes can be defined by using a new system for synthesis and screening of synthetic peptide libraries. These libraries are designed to bind to the HLA class II restriction molecule of the CD4+ T cell clone of interest. The screening is based on three selection rounds using partial release of 14-mer peptides from synthesis beads and subsequent sequencing of the remaining peptide attached to the bead. With this approach, two peptides were identified that stimulate the β cell-reactive CD4+ T cell clone 1c10, which was isolated from a newly diagnosed insulin-dependent diabetes mellitus patient. After performing amino acid-substitution studies and protein database searches, a Haemophilus influenzae TonB-derived peptide was identified that stimulates clone 1c10. The relevance of this finding for the pathogenesis of insulin-dependent diabetes mellitus is currently under investigation. We conclude that this system is capable of determining epitopes for (autoreactive) CD4+ T cell clones with previously unknown peptide specificity. This offers the possibility to define (auto)antigens by searching protein databases and/or to induce tolerance by using the peptide sequences identified. In addition the peptides might be used as leads to develop T cell receptor antagonists or anergy-inducing compounds.