Stochastic models for the in silico simulation of synaptic processes

Background: Research in life sciences is benefiting from a large availability of formal description techniques and analysis methodologies. These allow both the phenomena investigated to be precisely modeled and virtual experiments to be performed in silico. Such experiments may result in easier, faster, and satisfying approximations of their in vitro/vivo counterparts. A promising approach is represented by the study of biological phenomena as a collection of interactive entities through process calculi equipped with stochastic semantics. These exploit formal grounds developed in the theory of concurrency in computer science, account for the not continuous, nor discrete, nature of many phenomena, enjoy nice compositional properties and allow for simulations that have been demonstrated to be coherent with data in literature. Results: Motivated by the need to address some aspects of the functioning of neural synapses, we have developed one such model for synaptic processes in the calyx of Held, which is a glutamatergic synapse in the auditory pathway of the mammalia. We have developed such a stochastic model starting from existing kinetic models based on ODEs of some sub-components of the synapse, integrating other data from literature and making some assumptions about non-fully understood processes. Experiments have confirmed the coherence of our model with known biological data, also validating the assumptions made. Our model overcomes some limitations of the kinetic ones and, to our knowledge, represents the first model of synaptic processes based on process calculi. The compositionality of the approach has permitted us to independently focus on tuning the models of the pre- and post- synaptic traits, and then to naturally connect them, by dealing with “interface” issues. Furthermore, we have improved the expressiveness of the model, e.g. by embedding easy control of element concentration time courses. Sensitivity analysis over several parameters of the model has provided results that may help clarify the dynamics of synaptic transmission, while experiments with the model of the complete synapse seem worth explaining short-term plasticity mechanisms. Conclusions: Specific presynaptic and postsynaptic mechanisms can be further analysed under various conditions, for instance by studying the presynaptic behaviour under repeated activations. The level of details of the description can be refined, for instance by further specifying the neurotransmitter generation and release steps. Taking advantage of the compositionality of the approach, an enhanced model could then be composed with other neural models, designed within the same framework, in order to obtain a more detailed and comprehensive model. In the long term, we are interested, in particular, in addressing models of synaptic plasticity, i.e. activity dependent mechanisms, which are the bases of memory and learning processes. More on the computer science side, we plan to follow some directions to improve the underlying computational model and the linguistic primitives it provides as suggested by the experiments carried out, e.g. by introducing a suitable notion of (spatial) locality

STRASSE: A Silicon Tracker for Quasi-free Scattering Measurements at the RIBF

STRASSE (Silicon Tracker for RAdioactive nuclei Studies at SAMURAI Experiments) is a new detection system under construction for quasi-free scattering (QFS) measurements at 200-250 MeV/nucleon at the RIBF facility of the RIKEN Nishina Center. It consists of a charged-particle silicon tracker coupled with a dedicated thick liquid hydrogen target (up to 150-mm long) in a compact geometry to fit inside large scintillator or germanium arrays. Its design was optimized for two types of studies using QFS: missing-mass measurements and in-flight prompt $\gamma$-ray spectroscopy. This article describes (i) the resolution requirements needed to go beyond the sensitivity of existing systems for these two types of measurements, (ii) the conceptual design of the system using detailed simulations of the setup and (iii) its complete technical implementation and challenges. The final tracker aims at a sub-mm reaction vertex resolution and is expected to reach a missing-mass resolution below 2 MeV in $\sigma$ for $(p,2p)$ reactions when combined with the CsI(Na) CATANA array.Comment: 25 pages, 29 figure

Hydroxybenzothiazoles as New Nonsteroidal Inhibitors of 17β-Hydroxysteroid Dehydrogenase Type 1 (17β-HSD1)

17β-estradiol (E2), the most potent estrogen in humans, known to be involved in the development and progession of estrogen-dependent diseases (EDD) like breast cancer and endometriosis. 17β-HSD1, which catalyses the reduction of the weak estrogen estrone (E1) to E2, is often overexpressed in breast cancer and endometriotic tissues. An inhibition of 17β-HSD1 could selectively reduce the local E2-level thus allowing for a novel, targeted approach in the treatment of EDD. Continuing our search for new nonsteroidal 17β-HSD1 inhibitors, a novel pharmacophore model was derived from crystallographic data and used for the virtual screening of a small library of compounds. Subsequent experimental verification of the virtual hits led to the identification of the moderately active compound 5. Rigidification and further structure modifications resulted in the discovery of a novel class of 17β-HSD1 inhibitors bearing a benzothiazole-scaffold linked to a phenyl ring via keto- or amide-bridge. Their putative binding modes were investigated by correlating their biological data with features of the pharmacophore model. The most active keto-derivative 6 shows IC50-values in the nanomolar range for the transformation of E1 to E2 by 17β-HSD1, reasonable selectivity against 17β-HSD2 but pronounced affinity to the estrogen receptors (ERs). On the other hand, the best amide-derivative 21 shows only medium 17β-HSD1 inhibitory activity at the target enzyme as well as fair selectivity against 17β-HSD2 and ERs. The compounds 6 and 21 can be regarded as first benzothiazole-type 17β-HSD1 inhibitors for the development of potential therapeutics

Reelin Controls Progenitor Cell Migration in the Healthy and Pathological Adult Mouse Brain

Understanding the signals that control migration of neural progenitor cells in the adult brain may provide new therapeutic opportunities. Reelin is best known for its role in regulating cell migration during brain development, but we now demonstrate a novel function for reelin in the injured adult brain. First, we show that Reelin is upregulated around lesions. Second, experimentally increasing Reelin expression levels in healthy mouse brain leads to a change in the migratory behavior of subventricular zone-derived progenitors, triggering them to leave the rostral migratory stream (RMS) to which they are normally restricted during their migration to the olfactory bulb. Third, we reveal that Reelin increases endogenous progenitor cell dispersal in periventricular structures independently of any chemoattraction but via cell detachment and chemokinetic action, and thereby potentiates spontaneous cell recruitment to demyelination lesions in the corpus callosum. Conversely, animals lacking Reelin signaling exhibit reduced endogenous progenitor recruitment at the lesion site. Altogether, these results demonstrate that beyond its known role during brain development, Reelin is a key player in post-lesional cell migration in the adult brain. Finally our findings provide proof of concept that allowing progenitors to escape from the RMS is a potential therapeutic approach to promote myelin repair

Human embryonic stem cell-derived neurons establish region-specific, long-range projections in the adult brain

While the availability of pluripotent stem cells has opened new prospects for generating neural donor cells for nervous system repair, their capability to integrate with adult brain tissue in a structurally relevant way is still largely unresolved. We addressed the potential of human embryonic stem cell-derived long-term self-renewing neuroepithelial stem cells (lt-NES cells) to establish axonal projections after transplantation into the adult rodent brain. Transgenic and species-specific markers were used to trace the innervation pattern established by transplants in the hippocampus and motor cortex. In vitro, lt-NES cells formed a complex axonal network within several weeks after the initiation of differentiation and expressed a composition of surface receptors known to be instrumental in axonal growth and pathfinding. In vivo, these donor cells adopted projection patterns closely mimicking endogenous projections in two different regions of the adult rodent brain. Hippocampal grafts placed in the dentate gyrus projected to both the ipsilateral and contralateral pyramidal cell layers, while axons of donor neurons placed in the motor cortex extended via the external and internal capsule into the cervical spinal cord and via the corpus callosum into the contralateral cortex. Interestingly, acquisition of these region-specific projection profiles was not correlated with the adoption of a regional phenotype. Upon reaching their destination, human axons established ultrastructural correlates of synaptic connections with host neurons. Together, these data indicate that neurons derived from human pluripotent stem cells are endowed with a remarkable potential to establish orthotopic long-range projections in the adult mammalian brain

Microsecond Isomer at the N=20 Island of Shape Inversion Observed at FRIB

Excited-state spectroscopy from the first Facility for Rare Isotope Beams (FRIB) experiment is reported. A 24(2)-$\mu$s isomer was observed with the FRIB Decay Station initiator (FDSi) through a cascade of 224- and 401-keV $\gamma$ rays in coincidence with $^{32}\textrm{Na}$ nuclei. This is the only known microsecond isomer ($1{\text{ }\mu\text{s}}\leq T_{1/2} < 1\text{ ms}$) in the region. This nucleus is at the heart of the $N=20$ island of shape inversion and is at the crossroads of spherical shell-model, deformed shell-model, and ab initio theories. It can be represented as the coupling of a proton hole and neutron particle to $^{32}\textrm{Mg}$, $^{32}\textrm{Mg}+\pi^{-1} + \nu^{+1}$. This odd-odd coupling and isomer formation provides a sensitive measure of the underlying shape degrees of freedom of $^{32}\textrm{Mg}$, where the onset of spherical-to-deformed shape inversion begins with a low-lying deformed $2^+$ state at 885 keV and a low-lying shape-coexisting $0_2^+$ state at 1058 keV. We suggest two possible explanations for the 625-keV isomer in $^{32}$Na: a $6^-$ spherical shape isomer that decays by $E2$ or a $0^+$ deformed spin isomer that decays by $M2$. The present results and calculations are most consistent with the latter, indicating that the low-lying states are dominated by deformation.Comment: 7 pages, 5 figures, accepted by Physical Review Letter

Estrus cyclicity of spinogenesis: underlying mechanisms

Hippocampal spine density varies with the estrus cycle. The cyclic change in estradiol levels in serum was hypothesized to underlie this phenomenon, since treatment of ovariectomized animals with estradiol induced an increase in spine density in hippocampal dendrites of rats, as compared to ovariectomized controls. In contrast, application of estradiol to hippocampal slice cultures did not promote spinogenesis. In addressing this discrepancy, we found that hippocampal neurons themselves are capable of synthesizing estradiol de novo. Estradiol synthesis can be suppressed by aromatase inhibitors and by knock-down of Steroid Acute Regulatory Protein (StAR) and enhanced by substrates of steroidogenesis. Expression of estrogen receptors (ERs) and synaptic proteins, synaptogenesis, and long-term potentiation (LTP) correlated positively with aromatase activity in hippocampal cultures without any difference between genders. All effects due to inhibition of aromatase activity were rescued by application of estradiol to the cultures. Most importantly, gonadotropin-releasing hormone (GnRH) increased estradiol synthesis dose-dependently via an aromatase-mediated mechanism and consistently increased spine synapse density and spinophilin expression. As a consequence, our data suggest that cyclic fluctuations in spine synapse density result from pulsative release of GnRH from the hypothalamus and its effect on hippocampal estradiol synthesis, rather than from varying levels of serum estradiol. This hypothesis is further supported by higher GnRH receptor (GnRH-R) density in the hippocampus than in the cortex and hypothalamus and the specificity of estrus cyclicity of spinogenesis in the hippocampus, as compared to the cortex

Insights in 17β-HSD1 Enzyme Kinetics and Ligand Binding by Dynamic Motion Investigation

BACKGROUND: Bisubstrate enzymes, such as 17beta-hydroxysteroid dehydrogenase type 1 (17beta-HSD1), exist in solution as an ensemble of conformations. 17beta-HSD1 catalyzes the last step of the biosynthesis of estradiol and, thus, it is a potentially attractive target for breast cancer treatment. METHODOLOGY/PRINCIPAL FINDINGS: To elucidate the conformational transitions of its catalytic cycle, a structural analysis of all available crystal structures was performed and representative conformations were assigned to each step of the putative kinetic mechanism. To cover most of the conformational space, all-atom molecular dynamic simulations were performed using the four crystallographic structures best describing apoform, opened, occluded and closed state of 17beta-HSD1 as starting structures. With three of them, binary and ternary complexes were built with NADPH and NADPH-estrone, respectively, while two were investigated as apoform. Free energy calculations were performed in order to judge more accurately which of the MD complexes describes a specific kinetic step. CONCLUSIONS/SIGNIFICANCE: Remarkably, the analysis of the eight long range trajectories resulting from this multi-trajectory/-complex approach revealed an essential role played by the backbone and side chain motions, especially of the betaF alphaG'-loop, in cofactor and substrate binding. Thus, a selected-fit mechanism is suggested for 17beta-HSD1, where ligand-binding induced concerted motions of the FG-segment and the C-terminal part guide the enzyme along its preferred catalytic pathway. Overall, we could assign different enzyme conformations to the five steps of the random bi-bi kinetic cycle of 17beta-HSD1 and we could postulate a preferred pathway for it. This study lays the basis for more-targeted biochemical studies on 17beta-HSD1, as well as for the design of specific inhibitors of this enzyme. Moreover, it provides a useful guideline for other enzymes, also characterized by a rigid core and a flexible region directing their catalysis

Species Used for Drug Testing Reveal Different Inhibition Susceptibility for 17beta-Hydroxysteroid Dehydrogenase Type 1

Steroid-related cancers can be treated by inhibitors of steroid metabolism. In searching for new inhibitors of human 17beta-hydroxysteroid dehydrogenase type 1 (17β-HSD 1) for the treatment of breast cancer or endometriosis, novel substances based on 15-substituted estrone were validated. We checked the specificity for different 17β-HSD types and species. Compounds were tested for specificity in vitro not only towards recombinant human 17β-HSD types 1, 2, 4, 5 and 7 but also against 17β-HSD 1 of several other species including marmoset, pig, mouse, and rat. The latter are used in the processes of pharmacophore screening. We present the quantification of inhibitor preferences between human and animal models. Profound differences in the susceptibility to inhibition of steroid conversion among all 17β-HSDs analyzed were observed. Especially, the rodent 17β-HSDs 1 were significantly less sensitive to inhibition compared to the human ortholog, while the most similar inhibition pattern to the human 17β-HSD 1 was obtained with the marmoset enzyme. Molecular docking experiments predicted estrone as the most potent inhibitor. The best performing compound in enzymatic assays was also highly ranked by docking scoring for the human enzyme. However, species-specific prediction of inhibitor performance by molecular docking was not possible. We show that experiments with good candidate compounds would out-select them in the rodent model during preclinical optimization steps. Potentially active human-relevant drugs, therefore, would no longer be further developed. Activity and efficacy screens in heterologous species systems must be evaluated with caution

The Reelin Receptors Apoer2 and Vldlr Coordinate the Patterning of Purkinje Cell Topography in the Developing Mouse Cerebellum

The adult cerebellar cortex is comprised of reproducible arrays of transverse zones and parasagittal stripes of Purkinje cells. Adult stripes are created through the perinatal rostrocaudal dispersion of embryonic Purkinje cell clusters, triggered by signaling through the Reelin pathway. Reelin is secreted by neurons in the external granular layer and deep cerebellar nuclei and binds to two high affinity extracellular receptors on Purkinje cells-the Very low density lipoprotein receptor (Vldlr) and apolipoprotein E receptor 2 (Apoer2). In mice null for either Reelin or double null for Vldlr and Apoer2, Purkinje cell clusters fail to disperse. Here we report that animals null for either Vldlr or Apoer2 individually, exhibit specific and parasagittally-restricted Purkinje cell ectopias. For example, in mice lacking Apoer2 function immunostaining reveals ectopic Purkinje cells that are largely restricted to the zebrin II-immunonegative population of the anterior vermis. In contrast, mice null for Vldlr have a much larger population of ectopic Purkinje cells that includes members from both the zebrin II-immunonegative and -immunopositive phenotypes. HSP25 immunoreactivity reveals that in Vldlr null animals a large portion of zebrin II-immunopositive ectopic cells are probably destined to become stripes in the central zone (lobules VI–VII). A small population of ectopic zebrin II-immunonegative Purkinje cells is also observed in animals heterozygous for both receptors (Apoer2+/−: Vldlr+/−), but no ectopia is present in mice heterozygous for either receptor alone. These results indicate that Apoer2 and Vldlr coordinate the dispersal of distinct, but overlapping subsets of Purkinje cells in the developing cerebellum