Reconstruction of 1,000 Projection Neurons Reveals New Cell Types and Organization of Long-Range Connectivity in the Mouse Brain.

Neuronal cell types are the nodes of neural circuits that determine the flow of information within the brain. Neuronal morphology, especially the shape of the axonal arbor, provides an essential descriptor of cell type and reveals how individual neurons route their output across the brain. Despite the importance of morphology, few projection neurons in the mouse brain have been reconstructed in their entirety. Here we present a robust and efficient platform for imaging and reconstructing complete neuronal morphologies, including axonal arbors that span substantial portions of the brain. We used this platform to reconstruct more than 1,000 projection neurons in the motor cortex, thalamus, subiculum, and hypothalamus. Together, the reconstructed neurons constitute more than 85 meters of axonal length and are available in a searchable online database. Axonal shapes revealed previously unknown subtypes of projection neurons and suggest organizational principles of long-range connectivity.

Visually Guided Behavior and Optogenetically Induced Learning in Head-Fixed Flies Exploring a Virtual Landscape.

Studying the intertwined roles of sensation, experience, and directed action in navigation has been facilitated by the development of virtual reality (VR) environments for head-fixed animals, allowing for quantitative measurements of behavior in well-controlled conditions. VR has long featured in studies of Drosophila melanogaster, but these experiments have typically allowed the fly to change only its heading in a visual scene and not its position. Here we explore how flies move in two dimensions (2D) using a visual VR environment that more closely captures an animal's experience during free behavior. We show that flies' 2D interaction with landmarks cannot be automatically derived from their orienting behavior under simpler one-dimensional (1D) conditions. Using novel paradigms, we then demonstrate that flies in 2D VR adapt their behavior in response to optogenetically delivered appetitive and aversive stimuli. Much like free-walking flies after encounters with food, head-fixed flies exploring a 2D VR respond to optogenetic activation of sugar-sensing neurons by initiating a local search, which appears not to rely on visual landmarks. Visual landmarks can, however, help flies to avoid areas in VR where they experience an aversive, optogenetically generated heat stimulus. By coupling aversive virtual heat to the flies' presence near visual landmarks of specific shapes, we elicit selective learned avoidance of those landmarks. Thus, we demonstrate that head-fixed flies adaptively navigate in 2D virtual environments, but their reliance on visual landmarks is context dependent. These behavioral paradigms set the stage for interrogation of the fly brain circuitry underlying flexible navigation in complex multisensory environments.

OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular Simulation.

OpenMM is a software toolkit for performing molecular simulations on a range of high performance computing architectures. It is based on a layered architecture: the lower layers function as a reusable library that can be invoked by any application, while the upper layers form a complete environment for running molecular simulations. The library API hides all hardware-specific dependencies and optimizations from the users and developers of simulation programs: they can be run without modification on any hardware on which the API has been implemented. The current implementations of OpenMM include support for graphics processing units using the OpenCL and CUDA frameworks. In addition, OpenMM was designed to be extensible, so new hardware architectures can be accommodated and new functionality (e.g., energy terms and integrators) can be easily added.

Fast flexible modeling of RNA structure using internal coordinates.

Modeling the structure and dynamics of large macromolecules remains a critical challenge. Molecular dynamics (MD) simulations are expensive because they model every atom independently, and are difficult to combine with experimentally derived knowledge. Assembly of molecules using fragments from libraries relies on the database of known structures and thus may not work for novel motifs. Coarse-grained modeling methods have yielded good results on large molecules but can suffer from difficulties in creating more detailed full atomic realizations. There is therefore a need for molecular modeling algorithms that remain chemically accurate and economical for large molecules, do not rely on fragment libraries, and can incorporate experimental information. RNABuilder works in the internal coordinate space of dihedral angles and thus has time requirements proportional to the number of moving parts rather than the number of atoms. It provides accurate physics-based response to applied forces, but also allows user-specified forces for incorporating experimental information. A particular strength of RNABuilder is that all Leontis-Westhof basepairs can be specified as primitives by the user to be satisfied during model construction. We apply RNABuilder to predict the structure of an RNA molecule with 160 bases from its secondary structure, as well as experimental information. Our model matches the known structure to 10.2 Angstroms RMSD and has low computational expense.

Accelerating molecular dynamic simulation on graphics processing units.

We describe a complete implementation of all-atom protein molecular dynamics running entirely on a graphics processing unit (GPU), including all standard force field terms, integration, constraints, and implicit solvent. We discuss the design of our algorithms and important optimizations needed to fully take advantage of a GPU. We evaluate its performance, and show that it can be more than 700 times faster than a conventional implementation running on a single CPU core.

Structural basis for isozyme-specific regulation of electron transfer in nitric-oxide synthase.

Three nitric-oxide synthase (NOS) isozymes play crucial, but distinct, roles in neurotransmission, vascular homeostasis, and host defense, by catalyzing Ca(2+)/calmodulin-triggered NO synthesis. Here, we address current questions regarding NOS activity and regulation by combining mutagenesis and biochemistry with crystal structure determination of a fully assembled, electron-supplying, neuronal NOS reductase dimer. By integrating these results, we structurally elucidate the unique mechanisms for isozyme-specific regulation of electron transfer in NOS. Our discovery of the autoinhibitory helix, its placement between domains, and striking similarities with canonical calmodulin-binding motifs, support new mechanisms for NOS inhibition. NADPH, isozyme-specific residue Arg(1400), and the C-terminal tail synergistically repress NOS activity by locking the FMN binding domain in an electron-accepting position. Our analyses suggest that calmodulin binding or C-terminal tail phosphorylation frees a large scale swinging motion of the entire FMN domain to deliver electrons to the catalytic module in the holoenzyme.

Characterization of the electrophile binding site and substrate binding mode of the 26-kDa glutathione S-transferase from Schistosoma japonicum.

The 26-kDa glutathione S-transferase from Schistosoma japonicum (Sj26GST), a helminth worm that causes schistosomiasis, catalyzes the conjugation of glutathione with toxic secondary products of membrane lipid peroxidation. Crystal structures of Sj26GST in complex with glutathione sulfonate (Sj26GSTSLF), S-hexyl glutathione (Sj26GSTHEX), and S-2-iodobenzyl glutathione (Sj26GSTIBZ) allow characterization of the electrophile binding site (H site) of Sj26GST. The S-hexyl and S-2-iodobenzyl moieties of these product analogs bind in a pocket defined by side-chains from the beta1-alpha1 loop (Tyr7, Trp8, Ile10, Gly12, Leu13), helix alpha4 (Arg103, Tyr104, Ser107, Tyr111), and the C-terminal coil (Gln204, Gly205, Trp206, Gln207). Changes in the Ser107 and Gln204 dihedral angles make the H site more hydrophobic in the Sj26GSTHEX complex relative to the ligand-free structure. These structures, together with docking studies, indicate a possible binding mode of Sj26GST to its physiologic substrates 4-hydroxynon-2-enal (4HNE), trans-non-2-enal (NE), and ethacrynic acid (EA). In this binding mode, hydrogen bonds of Tyr111 and Gln207 to the carbonyl oxygen atoms of 4HNE, NE, and EA could orient the substrates and enhance their electrophilicity to promote conjugation with glutathione.

Design of a monomeric human glutathione transferase GSTP1, a structurally stable but catalytically inactive protein.

By the introduction of 10 site-specific mutations in the dimer interface of human glutathione transferase P1-1 (GSTP1-1), a stable monomeric protein variant, GSTP1, was obtained. The monomer had lost the catalytic activity but retained the affinity for a number of electrophilic compounds normally serving as substrates for GSTP1-1. Fluorescence and circular dichroism spectra of the monomer and wild-type proteins were similar, indicating that there are no large structural differences between the subunits of the respective proteins. The GSTs have potential as targets for in vitro evolution and redesign with the aim of developing proteins with novel properties. To this end, a monomeric GST variant may have distinct advantages.