Mini-review: The neurobiology of treating substance use disorders with classical psychedelics

The potential of psychedelics to persistently treat substance use disorders is known since the 1960s. However, the biological mechanisms responsible for their therapeutic effects have not yet been fully elucidated. While it is known that serotonergic hallucinogens induce changes in gene expression and neuroplasticity, particularly in prefrontal regions, theories on how specifically this counteracts the alterations that occur in neuronal circuitry throughout the course of addiction are largely unknown. This narrative mini-review endeavors to synthesize well-established knowledge from addiction research with findings and theories regarding the neurobiological effects of psychedelics to give an overview of the potential mechanisms that underlie the treatment of substance use disorders with classical hallucinogenic compounds and point out gaps in the current understanding

Sierpinski signal generates 1/fα1/f^\alpha spectra

We investigate the row sum of the binary pattern generated by the Sierpinski automaton: Interpreted as a time series we calculate the power spectrum of this Sierpinski signal analytically and obtain a unique rugged fine structure with underlying power law decay with an exponent of approximately 1.15. Despite the simplicity of the model, it can serve as a model for $1/f^\alpha$ spectra in a certain class of experimental and natural systems like catalytic reactions and mollusc patterns.Comment: 4 pages (4 figs included). Accepted for publication in Physical Review

Existence and Stability of a Spike in the Central Component for a Consumer Chain Model

We study a three-component consumer chain model which is based on Schnakenberg type kinetics. In this model there is one consumer feeding on the producer and a second consumer feeding on the first consumer. This means that the first consumer (central component) plays a hybrid role: it acts both as consumer and producer. The model is an extension of the Schnakenberg model suggested in \cite{gm,schn1} for which there is only one producer and one consumer. It is assumed that both the producer and second consumer diffuse much faster than the central component. We construct single spike solutions on an interval for which the profile of the first consumer is that of a spike. The profiles of the producer and the second consumer only vary on a much larger spatial scale due to faster diffusion of these components. It is shown that there exist two different single spike solutions if the feed rates are small enough: a large-amplitude and a small-amplitude spike. We study the stability properties of these solutions in terms of the system parameters. We use a rigorous analysis for the linearized operator around single spike solutions based on nonlocal eigenvalue problems. The following result is established: If the time-relaxation constants for both producer and second consumer vanish, the large-amplitude spike solution is stable and the small-amplitude spike solution is unstable. We also derive results on the stability of solutions when these two time-relaxation constants are small. We show a new effect: if the time-relaxation constant of the second consumer is very small, the large-amplitude spike solution becomes unstable. To the best of our knowledge this phenomenon has not been observed before for the stability of spike patterns. It seems that this behavior is not possible for two-component reaction-diffusion systems but that at least three components are required. Our main motivation to study this system is mathematical since the novel interaction of a spike in the central component with two other components results in new types of conditions for the existence and stability of a spike. This model is realistic if several assumptions are made: (i) cooperation of consumers is prevalent in the system, (ii) the producer and the second consumer diffuse much faster than the first consumer, and (iii) there is practically an unlimited pool of producer. The first assumption has been proven to be correct in many types of consumer groups or populations, the second assumption occurs if the central component has a much smaller mobility than the other two, the third assumption is realistic if the consumers do not feel the impact of the limited amount of producer due to its large quantity. This chain model plays a role in population biology, where consumer and producer are often called predator and prey. This system can also be used as a model for a sequence of irreversible autocatalytic reactions in a container which is in contact with a well-stirred reservoir

Subbarrel patterns in somatosensory cortical barrels can emerge from local dynamic instabilities

Complex spatial patterning, common in the brain as well as in other biological systems, can emerge as a result of dynamic interactions that occur locally within developing structures. In the rodent somatosensory cortex, groups of neurons called "barrels" correspond to individual whiskers on the contralateral face. Barrels themselves often contain subbarrels organized into one of a few characteristic patterns. Here we demonstrate that similar patterns can be simulated by means of local growth-promoting and growth-retarding interactions within the circular domains of single barrels. The model correctly predicts that larger barrels contain more spatially complex subbarrel patterns, suggesting that the development of barrels and of the patterns within them may be understood in terms of some relatively simple dynamic processes. We also simulate the full nonlinear equations to demonstrate the predictive value of our linear analysis. Finally, we show that the pattern formation is robust with respect to the geometry of the barrel by simulating patterns on a realistically shaped barrel domain. This work shows how simple pattern forming mechanisms can explain neural wiring both qualitatively and quantitatively even in complex and irregular domains. © 2009 Ermentrout et al

Perturbative Linearization of Reaction-Diffusion Equations

We develop perturbative expansions to obtain solutions for the initial-value problems of two important reaction-diffusion systems, viz., the Fisher equation and the time-dependent Ginzburg-Landau (TDGL) equation. The starting point of our expansion is the corresponding singular-perturbation solution. This approach transforms the solution of nonlinear reaction-diffusion equations into the solution of a hierarchy of linear equations. Our numerical results demonstrate that this hierarchy rapidly converges to the exact solution.Comment: 13 pages, 4 figures, latex2

Stability of cluster solutions in a cooperative consumer chain model

This is the author's accepted manuscript. The final published article is available from the link below. Copyright @ Springer-Verlag Berlin Heidelberg 2012.We study a cooperative consumer chain model which consists of one producer and two consumers. It is an extension of the Schnakenberg model suggested in Gierer and Meinhardt [Kybernetik (Berlin), 12:30-39, 1972] and Schnakenberg (J Theor Biol, 81:389-400, 1979) for which there is only one producer and one consumer. In this consumer chain model there is a middle component which plays a hybrid role: it acts both as consumer and as producer. It is assumed that the producer diffuses much faster than the first consumer and the first consumer much faster than the second consumer. The system also serves as a model for a sequence of irreversible autocatalytic reactions in a container which is in contact with a well-stirred reservoir. In the small diffusion limit we construct cluster solutions in an interval which have the following properties: The spatial profile of the third component is a spike. The profile for the middle component is that of two partial spikes connected by a thin transition layer. The first component in leading order is given by a Green's function. In this profile multiple scales are involved: The spikes for the middle component are on the small scale, the spike for the third on the very small scale, the width of the transition layer for the middle component is between the small and the very small scale. The first component acts on the large scale. To the best of our knowledge, this type of spiky pattern has never before been studied rigorously. It is shown that, if the feedrates are small enough, there exist two such patterns which differ by their amplitudes.We also study the stability properties of these cluster solutions. We use a rigorous analysis to investigate the linearized operator around cluster solutions which is based on nonlocal eigenvalue problems and rigorous asymptotic analysis. The following result is established: If the time-relaxation constants are small enough, one cluster solution is stable and the other one is unstable. The instability arises through large eigenvalues of order O(1). Further, there are small eigenvalues of order o(1) which do not cause any instabilities. Our approach requires some new ideas: (i) The analysis of the large eigenvalues of order O(1) leads to a novel system of nonlocal eigenvalue problems with inhomogeneous Robin boundary conditions whose stability properties have been investigated rigorously. (ii) The analysis of the small eigenvalues of order o(1) needs a careful study of the interaction of two small length scales and is based on a suitable inner/outer expansion with rigorous error analysis. It is found that the order of these small eigenvalues is given by the smallest diffusion constant ε22.RGC of Hong Kon

Molecular Characterization of the Gastrula in the Turtle Emys orbicularis: An Evolutionary Perspective on Gastrulation

Due to the presence of a blastopore as in amphibians, the turtle has been suggested to exemplify a transition form from an amphibian- to an avian-type gastrulation pattern. In order to test this hypothesis and gain insight into the emergence of the unique characteristics of amniotes during gastrulation, we have performed the first molecular characterization of the gastrula in a reptile, the turtle Emys orbicularis. The study of Brachyury, Lim1, Otx2 and Otx5 expression patterns points to a highly conserved dynamic of expression with amniote model organisms and makes it possible to identify the site of mesoderm internalization, which is a long-standing issue in reptiles. Analysis of Brachyury expression also highlights the presence of two distinct phases, less easily recognizable in model organisms and respectively characterized by an early ring-shaped and a later bilateral symmetrical territory. Systematic comparisons with tetrapod model organisms lead to new insights into the relationships of the blastopore/blastoporal plate system shared by all reptiles, with the blastopore of amphibians and the primitive streak of birds and mammals. The biphasic Brachyury expression pattern is also consistent with recent models of emergence of bilateral symmetry, which raises the question of its evolutionary significance

Using Hyperfine Electron Paramagnetic Resonance Spectroscopy to Define the Proton-Coupled Electron Transfer Reaction at Fe-S Cluster N2 in Respiratory Complex I.

Energy-transducing respiratory complex I (NADH:ubiquinone oxidoreductase) is one of the largest and most complicated enzymes in mammalian cells. Here, we used hyperfine electron paramagnetic resonance (EPR) spectroscopic methods, combined with site-directed mutagenesis, to determine the mechanism of a single proton-coupled electron transfer reaction at one of eight iron-sulfur clusters in complex I, [4Fe-4S] cluster N2. N2 is the terminal cluster of the enzyme's intramolecular electron-transfer chain and the electron donor to ubiquinone. Because of its position and pH-dependent reduction potential, N2 has long been considered a candidate for the elusive "energy-coupling" site in complex I at which energy generated by the redox reaction is used to initiate proton translocation. Here, we used hyperfine sublevel correlation (HYSCORE) spectroscopy, including relaxation-filtered hyperfine and single-matched resonance transfer (SMART) HYSCORE, to detect two weakly coupled exchangeable protons near N2. We assign the larger coupling with A(1H) = [-3.0, -3.0, 8.7] MHz to the exchangeable proton of a conserved histidine and conclude that the histidine is hydrogen-bonded to N2, tuning its reduction potential. The histidine protonation state responds to the cluster oxidation state, but the two are not coupled sufficiently strongly to catalyze a stoichiometric and efficient energy transduction reaction. We thus exclude cluster N2, despite its proton-coupled electron transfer chemistry, as the energy-coupling site in complex I. Our work demonstrates the capability of pulse EPR methods for providing detailed information on the properties of individual protons in even the most challenging of energy-converting enzymes

Measurement and simulation of the 16/17 April 2010 Eyjafjallajökull volcanic ash layer dispersion in the northern Alpine region

The spatial structure and the progression speed of the first ash layer from the Icelandic Eyjafjallajökull volcano which reached Germany on 16/17 April is investigated from remote sensing data and numerical simulations. The ceilometer network of the German Meteorological Service was able to follow the progression of the ash layer over the whole of Germany. This first ash layer turned out to be a rather shallow layer of only several hundreds of metres thickness which was oriented slantwise in the middle troposphere and which was brought downward by large-scale sinking motion over Southern Germany and the Alps. Special Raman lidar measurements, trajectory analyses and in-situ observations from mountain observatories helped to confirm the volcanic origin of the detected aerosol layer. Ultralight aircraft measurements permitted the detection of the arrival of a second major flush of volcanic material in Southern Germany. Numerical simulations with the Eulerian meso-scale model MCCM were able to reproduce the temporal and spatial structure of the ash layer. Comparisons of the model results with the ceilometer network data on 17 April and with the ultralight aircraft data on 19 April were satisfying. This is the first example of a model validation study from this ceilometer network data