Quantitative assessment of computational models for retinotopic map formation.

Molecular and activity-based cues acting together are thought to guide retinal axons to their terminal sites in vertebrate optic tectum or superior colliculus (SC) to form an ordered map of connections. The details of mechanisms involved, and the degree to which they might interact, are still not well understood. We have developed a framework within which existing computational models can be assessed in an unbiased and quantitative manner against a set of experimental data curated from the mouse retinocollicular system. Our framework facilitates comparison between models, testing new models against known phenotypes and simulating new phenotypes in existing models. We have used this framework to assess four representative models that combine Eph/ephrin gradients and/or activity-based mechanisms and competition. Two of the models were updated from their original form to fit into our framework. The models were tested against five different phenotypes: wild type, Isl2-EphA3(ki/ki), Isl2-EphA3(ki/+), ephrin-A2,A3,A5 triple knock-out (TKO), and Math5(-/-) (Atoh7). Two models successfully reproduced the extent of the Math5(-/-) anteromedial projection, but only one of those could account for the collapse point in Isl2-EphA3(ki/+). The models needed a weak anteroposterior gradient in the SC to reproduce the residual order in the ephrin-A2,A3,A5 TKO phenotype, suggesting either an incomplete knock-out or the presence of another guidance molecule. Our article demonstrates the importance of testing retinotopic models against as full a range of phenotypes as possible, and we have made available MATLAB software, we wrote to facilitate this process.Contract grant sponsors: Wellcome Trust; contract grant numbers: 083205, Engineering and Physical Sciences Research Council (EPSRC) (CSC).This is the published version. It's also available from Wiley at http://onlinelibrary.wiley.com/doi/10.1002/dneu.22241/abstract

Estimating the location and size of retinal injections from orthogonal images of an intact retina.

BACKGROUND: To study the mapping from the retina to the brain, typically a small region of the retina is injected with a dye, which then propagates to the retina's target structures. To determine the location of the injection, usually the retina is dissected out of the eye, flattened and then imaged, causing tears and stretching of the retina. The location of the injection is then estimated from the image of the flattened retina. Here we propose a new method that avoids dissection of the retina. RESULTS: We have developed IntactEye, a software package that uses two orthogonal images of the intact retina to locate focal injections of a dye. The two images are taken while the retina is still inside the eye. This bypasses the dissection step, avoiding unnecessary damage to the retina, and speeds up data acquisition. By using the native spherical coordinates of the eye, we avoid distortions caused by interpreting a curved structure in a flat coordinate system. Our method compares well to the projection method and to the Retistruct package, which both use the flattened retina as a starting point. We have tested the method also on synthetic data, where the injection location is known. Our method has been designed for analysing mouse retinas, where there are no visible landmarks for discerning retinal orientation, but can also be applied to retinas from other species. CONCLUSIONS: IntactEye allows the user to precisely specify the location and size of a retinal injection from two orthogonal images taken of the eye. We are solving the abstract problem of locating a point on a spherical object from two orthogonal images, which might have applications outside the field of neuroscience.SJE and MR gratefully acknowledge the support of the University of Strasbourg Institute for Advanced Study (USIAS). SJE and JJJH were supported by the Wellcome Trust (grant number 083205). The authors wish to thank Ellese Cotterill for analysing synthetic data for verification of accuracy

Human Synapses Show a Wide Temporal Window for Spike-Timing-Dependent Plasticity

Throughout our lifetime, activity-dependent changes in neuronal connection strength enable the brain to refine neural circuits and learn based on experience. Synapses can bi-directionally alter strength and the magnitude and sign depend on the millisecond timing of presynaptic and postsynaptic action potential firing. Recent findings on laboratory animals have shown that neurons can show a variety of temporal windows for spike-timing-dependent plasticity (STDP). It is unknown what synaptic learning rules exist in human synapses and whether similar temporal windows for STDP at synapses hold true for the human brain. Here, we directly tested in human slices cut from hippocampal tissue removed for surgical treatment of deeper brain structures in drug-resistant epilepsy patients, whether adult human synapses can change strength in response to millisecond timing of pre- and postsynaptic firing. We find that adult human hippocampal synapses can alter synapse strength in response to timed pre- and postsynaptic activity. In contrast to rodent hippocampal synapses, the sign of plasticity does not sharply switch around 0-ms timing. Instead, both positive timing intervals, in which presynaptic firing preceded the postsynaptic action potential, and negative timing intervals, in which postsynaptic firing preceded presynaptic activity down to −80 ms, increase synapse strength (tLTP). Negative timing intervals between −80 to −130 ms induce a lasting reduction of synapse strength (tLTD). Thus, similar to rodent synapses, adult human synapses can show spike-timing-dependent changes in strength. The timing rules of STDP in human hippocampus, however, seem to differ from rodent hippocampus, and suggest a less strict interpretation of Hebb's predictions

Synapto-protective drugs evaluation in reconstructed neuronal network

Chronic neurodegenerative syndromes such as Alzheimer’s and Parkinson’s diseases, or acute syndromes such as ischemic stroke or traumatic brain injuries are characterized by early synaptic collapse which precedes axonal and neuronal cell body degeneration and promotes early cognitive impairment in patients. Until now, neuroprotective strategies have failed to impede the progression of neurodegenerative syndromes. Drugs preventing the loss of cell body do not prevent the cognitive decline, probably because they lack synapto-protective effects. The absence of physiologically realistic neuronal network models which can be easily handled has hindered the development of synapto-protective drugs suitable for therapies. Here we describe a new microfluidic platform which makes it possible to study the consequences of axonal trauma of reconstructed oriented mouse neuronal networks. Each neuronal population and sub-compartment can be chemically addressed individually. The somatic, mid axon, presynaptic and postsynaptic effects of local pathological stresses or putative protective molecules can thus be evaluated with the help of this versatile “brain on chip” platform. We show that presynaptic loss is the earliest event observed following axotomy of cortical fibers, before any sign of axonal fragmentation or post-synaptic spine alteration. This platform can be used to screen and evaluate the synapto-protective potential of several drugs. For instance, NAD+ and the Rho-kinase inhibitor Y27632 can efficiently prevent synaptic disconnection, whereas the broad-spectrum caspase inhibitor zVAD-fmk and the stilbenoid resveratrol do not prevent presynaptic degeneration. Hence, this platform is a promising tool for fundamental research in the field of developmental and neurodegenerative neurosciences, and also offers the opportunity to set up pharmacological screening of axon-protective and synapto-protective drugs

Density functional theory based screening of ternary alkali-transition metal borohydrides: A computational material design project

The dissociation of molecules, even the most simple hydrogen molecule, cannot be described accurately within density functional theory because none of the currently available functionals accounts for strong on-site correlation. This problem led to a discussion of properties that the local Kohn-Sham potential has to satisfy in order to correctly describe strongly correlated systems. We derive an analytic expression for the nontrivial form of the Kohn-Sham potential in between the two fragments for the dissociation of a single bond. We show that the numerical calculations for a one-dimensional two-electron model system indeed approach and reach this limit. It is shown that the functional form of the potential is universal, i.e., independent of the details of the two fragments.We acknowledge funding by the Spanish MEC (Grant No. FIS2007-65702-C02-01), “Grupos Consolidados UPV/EHU del Gobierno Vasco” (Grant No. IT-319-07), and the European Community through e-I3 ETSF project (Grant Agreement No. 211956).Peer reviewe

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light–matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).This work was supported by the European Research Council (Grant No. ERC-2015-AdG694097), the Cluster of Excellence “Advanced Imaging of Matter” (AIM), Grupos Consolidados (IT1249-19), and SFB925. The Flatiron Institute is a division of the Simons Foundation. X.A., A.W., and A.C. acknowledge that part of this work was performed under the auspices of the U.S. Department of Energy at Lawrence Livermore National Laboratory under Contract No. DE-AC52-07A27344. J.J.-S. gratefully acknowledges the funding from the European Union Horizon 2020 Research and Innovation Program under the Marie Sklodowska-Curie Grant Agreement No. 795246-StrongLights. J.F. acknowledges financial support from the Deutsche Forschungsgemeinschaft (DFG Forschungsstipendium FL 997/1-1). D.A.S. acknowledges University of California, Merced start-up funding.Peer reviewe

Competition between neurite branches in a complex morphology.

<p>(A) Example morphology of a reconstructed pyramidal neuron with apical and basal dendrites. (B) In the control case, starting from the reconstructed morphology, the neuron was allowed to grow out for 10 hours in the model. The simulation was then repeated with the same initial conditions, but with increased polymerization rate for one of the growth cones. The dendritic morphology obtained in this last simulation is represented by a dendrogram, colored according to the tubulin concentration in the branches. The gray vertical lines at the terminal segments indicate the starting morphology, and the black vertical lines show the neurite length after 10 hours in the control case. The black dot marks the growth cone with increased polymerization rate. (C) The competition between branches increases with increasing path distance to the soma. The graph shows the total retraction of all neurites, divided by the growth of the modified growth cone, as a function of path length between the modified growth cone and the soma.</p

Default and optimized values of the model parameters.

<p>Default and optimized values of the model parameters.</p