Many-particle Brownian and Langevin Dynamics Simulations with the Brownmove package

<p>Abstract</p> <p>Background</p> <p>Brownian Dynamics (BD) is a coarse-grained implicit-solvent simulation method that is routinely used to investigate binary protein association dynamics, but due to its efficiency in handling large simulation volumes and particle numbers it is well suited to also describe many-protein scenarios as they often occur in biological cells.</p> <p>Results</p> <p>Here we introduce our "brownmove" simulation package which was designed to handle many-particle problems with varying particle numbers and allows for a very flexible definition of rigid and flexible protein and polymer models. Both a Brownian and a Langevin dynamics (LD) propagation scheme can be used and hydrodynamic interactions are treated efficiently with our recently introduced TEA-HI ansatz [Geyer, Winter, JCP 130 (2009) 114905]. With simulations of constrained polymers and flexible models of spherical proteins we demonstrate that it is crucial to include hydrodynamics when multi-bead models are used in BD or LD simulations. Only then both the translational and the rotational diffusion coefficients and the timescales of the internal dynamics can be reproduced correctly. In the third example project we show how constant density boundary conditions [Geyer et al, JCP 120 (2004) 4573] can be used to set up a non-equilibrium simulation of diffusional transport across an array of fixed obstacles. Finally, we demonstrate how the agglomeration dynamics of multiple particles with attractive patches can be analysed conveniently with the help of a dynamic interaction network.</p> <p>Conclusions</p> <p>Combining BD and LD propagation, fast hydrodynamics, a flexible protein model, and interfaces for "open" simulation settings, our freely available "brownmove" simulation package constitutes a new platform for coarse-grained many-particle simulations of biologically relevant diffusion and transport processes.</p

Intrinsically determined cell death of developing cortical interneurons

Cortical inhibitory circuits are formed by GABAergic interneurons, a cell population that originates far from the cerebral cortex in the embryonic ventral forebrain. Given their distant developmental origins, it is intriguing how the number of cortical interneurons is ultimately determined. One possibility, suggested by the neurotrophic hypothesis1-5, is that cortical interneurons are overproduced, and then following their migration into cortex, excess interneurons are eliminated through a competition for extrinsically derived trophic signals. Here we have characterized the developmental cell death of mouse cortical interneurons in vivo, in vitro, and following transplantation. We found that 40% of developing cortical interneurons were eliminated through Bax- (Bcl-2 associated X-) dependent apoptosis during postnatal life. When cultured in vitro or transplanted into the cortex, interneuron precursors died at a cellular age similar to that at which endogenous interneurons died during normal development. Remarkably, over transplant sizes that varied 200-fold, a constant fraction of the transplanted population underwent cell death. The death of transplanted neurons was not affected by the cell-autonomous disruption of TrkB (tropomyosin kinase receptor B), the main neurotrophin receptor expressed by central nervous system (CNS) neurons6-8. Transplantation expanded the cortical interneuron population by up to 35%, but the frequency of inhibitory synaptic events did not scale with the number of transplanted interneurons. Together, our findings indicate that interneuron cell death is intrinsically determined, either cell-autonomously, or through a population-autonomous competition for survival signals derived from other interneurons

Diffusion in crowded biological environments: applications of Brownian dynamics

Biochemical reactions in living systems occur in complex, heterogeneous media with total concentrations of macromolecules in the range of 50 - 400 mgml. Molecular species occupy a significant fraction of the immersing medium, up to 40% of volume. Such complex and volume-occupied environments are generally termed 'crowded' and/or 'confined'. In crowded conditions non-specific interactions between macromolecules may hinder diffusion - a major process determining metabolism, transport, and signaling. Also, the crowded media can alter, both qualitatively and quantitatively, the reactions in vivo in comparison with their in vitro counterparts. This review focuses on recent developments in particle-based Brownian dynamics algorithms, their applications to model diffusive transport in crowded systems, and their abilities to reproduce and predict the behavior of macromolecules under in vivo conditions

Neocortical Axon Arbors Trade-off Material and Conduction Delay Conservation

The brain contains a complex network of axons rapidly communicating information between billions of synaptically connected neurons. The morphology of individual axons, therefore, defines the course of information flow within the brain. More than a century ago, Ramón y Cajal proposed that conservation laws to save material (wire) length and limit conduction delay regulate the design of individual axon arbors in cerebral cortex. Yet the spatial and temporal communication costs of single neocortical axons remain undefined. Here, using reconstructions of in vivo labelled excitatory spiny cell and inhibitory basket cell intracortical axons combined with a variety of graph optimization algorithms, we empirically investigated Cajal's conservation laws in cerebral cortex for whole three-dimensional (3D) axon arbors, to our knowledge the first study of its kind. We found intracortical axons were significantly longer than optimal. The temporal cost of cortical axons was also suboptimal though far superior to wire-minimized arbors. We discovered that cortical axon branching appears to promote a low temporal dispersion of axonal latencies and a tight relationship between cortical distance and axonal latency. In addition, inhibitory basket cell axonal latencies may occur within a much narrower temporal window than excitatory spiny cell axons, which may help boost signal detection. Thus, to optimize neuronal network communication we find that a modest excess of axonal wire is traded-off to enhance arbor temporal economy and precision. Our results offer insight into the principles of brain organization and communication in and development of grey matter, where temporal precision is a crucial prerequisite for coincidence detection, synchronization and rapid network oscillations

Low-energy short-term cold atmospheric plasma ::controlling the inactivation efficacy of bacterial spores in powders

The present research work aims to elucidate kinetics and mechanisms of the inactivation of Bacillus subtilis spores by a surface micro-discharge (SMD) - cold atmospheric pressure plasma (CAPP). Regarding industrial applications, the inactivation of spores was also studied for a static layer of a biopolymer powder or film, with an air plasma and at ambient pressure. Close to 4 log10 cycles of inactivation of Bacillus subtilis spores were achieved when exposing spores on flat glass to the SMD-CAPP. This effect can be reached at a very low plasma power density of 5 mW/cm2 in 7 min exposure time. The maximum inactivation level of spores drops when treating corn-starch powder to 2.6 log10 cycles at 7 mW/cm2 plasma power density for 5 min and with a polymer load of 5 mg/cm2. Similar is true for films produced with hydroxymethyl cellulose (HMC). The inactivation efficacy can be tuned and is a function of applied surface energy (product of the plasma power density and the exposure time) and the polymer load. Plasma diagnostics reveal the fundamental importance of reactive nitrogen species (RNS) in the inactivation. Etching of spore hull is supposed to be triggered by the plasma density, while UV-C and UV-B radiation do not contribute directly and significantly to the inactivation effect at least in a biopolymer matrix. Fluidization of a fixed powder layer is supposed to overcome limitations of the inactivation efficacy by reducing the diffusion distance of active plasma species between the source and the sample. The combination of low plasma power density with short treatment time is supposed to reduce the risk of the formation of side-products from the matrix

Diffusional Encounter of Barnase and Barstar

We present an analysis of trajectories from Brownian dynamics simulations of diffusional protein-protein encounter for the well-studied system of barnase and barstar. This analysis reveals details about the optimal association pathways, the regions of the encounter complex, possible differences of the pathways for dissociation and association, the coupling of translational and rotation motion, and the effect of mutations on the trajectories. We found that a small free-energy barrier divides the energetically most favorable region into a region of the encounter complex above the barnase binding interface and a region around a second energy minimum near the RNA binding loop. When entering the region of the encounter complex from the region near the RNA binding loop, barstar has to change its orientation to increase the electrostatic attraction between the proteins. By concentrating the analysis on the successful binding trajectories, we found that the region of the second minimum is not essential for the binding of barstar to barnase. Nevertheless, this region may be helpful to steer barstar into the region of the encounter complex. When applying the same analysis to several barnase mutants, we found that single mutations may drastically change the free-energy landscape and may significantly alter the population of the two minima. Therefore, certain protein-protein pairs may require careful adaptation of the positions of encounter and transition states when interpreting mutation effects on kinetic rates of association and/or dissociation

A key mechanism underlying sensory experience-dependent maturation of neocortical GABAergic circuits in vivo

Mechanisms underlying experience-dependent refinement of cortical connections, especially GABAergic inhibitory circuits, are unknown. By using a line of mutant mice that lack activity-dependent BDNF expression (bdnf-KIV), we show that experience regulation of cortical GABAergic network is mediated by activity-driven BDNF expression. Levels of endogenous BDNF protein in the barrel cortex are strongly regulated by sensory inputs from whiskers. There is a severe alteration of excitation and inhibition balance in the barrel cortex of bdnf-KIV mice as a result of reduced inhibitory but not excitatory conductance. Within the inhibitory circuits, the mutant barrel cortex exhibits significantly reduced levels of GABA release only from the parvalbumin-expressing fast-spiking (FS) interneurons, but not other interneuron subtypes. Postnatal deprivation of sensory inputs markedly decreased perisomatic inhibition selectively from FS cells in wild-type but not bdnf-KIV mice. These results suggest that postnatal experience, through activity-driven BDNF expression, controls cortical development by regulating FS cell-mediated perisomatic inhibition in vivo