Journal Articles

Correction: Time-Course Gene Set Analysis for Longitudinal Gene Expression Data

PLoS Computational Biology - Fri, 08/21/2015 - 16:00

by The PLOS Computational Biology Staff

Categories: Journal Articles

Self-Organization of Microcircuits in Networks of Spiking Neurons with Plastic Synapses

PLoS Computational Biology - Thu, 08/20/2015 - 16:00

by Gabriel Koch Ocker, Ashok Litwin-Kumar, Brent Doiron

The synaptic connectivity of cortical networks features an overrepresentation of certain wiring motifs compared to simple random-network models. This structure is shaped, in part, by synaptic plasticity that promotes or suppresses connections between neurons depending on their joint spiking activity. Frequently, theoretical studies focus on how feedforward inputs drive plasticity to create this network structure. We study the complementary scenario of self-organized structure in a recurrent network, with spike timing-dependent plasticity driven by spontaneous dynamics. We develop a self-consistent theory for the evolution of network structure by combining fast spiking covariance with a slow evolution of synaptic weights. Through a finite-size expansion of network dynamics we obtain a low-dimensional set of nonlinear differential equations for the evolution of two-synapse connectivity motifs. With this theory in hand, we explore how the form of the plasticity rule drives the evolution of microcircuits in cortical networks. When potentiation and depression are in approximate balance, synaptic dynamics depend on weighted divergent, convergent, and chain motifs. For additive, Hebbian STDP these motif interactions create instabilities in synaptic dynamics that either promote or suppress the initial network structure. Our work provides a consistent theoretical framework for studying how spiking activity in recurrent networks interacts with synaptic plasticity to determine network structure.
Categories: Journal Articles

Experimentally Verified Parameter Sets for Modelling Heterogeneous Neocortical Pyramidal-Cell Populations

PLoS Computational Biology - Thu, 08/20/2015 - 16:00

by Paul M. Harrison, Laurent Badel, Mark J. Wall, Magnus J. E. Richardson

Models of neocortical networks are increasingly including the diversity of excitatory and inhibitory neuronal classes. Significant variability in cellular properties are also seen within a nominal neuronal class and this heterogeneity can be expected to influence the population response and information processing in networks. Recent studies have examined the population and network effects of variability in a particular neuronal parameter with some plausibly chosen distribution. However, the empirical variability and covariance seen across multiple parameters are rarely included, partly due to the lack of data on parameter correlations in forms convenient for model construction. To addess this we quantify the heterogeneity within and between the neocortical pyramidal-cell classes in layers 2/3, 4, and the slender-tufted and thick-tufted pyramidal cells of layer 5 using a combination of intracellular recordings, single-neuron modelling and statistical analyses. From the response to both square-pulse and naturalistic fluctuating stimuli, we examined the class-dependent variance and covariance of electrophysiological parameters and identify the role of the h current in generating parameter correlations. A byproduct of the dynamic I-V method we employed is the straightforward extraction of reduced neuron models from experiment. Empirically these models took the refractory exponential integrate-and-fire form and provide an accurate fit to the perisomatic voltage responses of the diverse pyramidal-cell populations when the class-dependent statistics of the model parameters were respected. By quantifying the parameter statistics we obtained an algorithm which generates populations of model neurons, for each of the four pyramidal-cell classes, that adhere to experimentally observed marginal distributions and parameter correlations. As well as providing this tool, which we hope will be of use for exploring the effects of heterogeneity in neocortical networks, we also provide the code for the dynamic I-V method and make the full electrophysiological data set available.
Categories: Journal Articles

Computational scheme for pH-dependent binding free energy calculation with explicit solvent

Protein Science - Thu, 08/20/2015 - 13:17
Abstract

We present a computational scheme to compute the pH-dependence of binding free energy with explicit solvent. Despite the importance of pH, the effect of pH has been generally neglected in binding free energy calculations because of a lack of accurate methods to model it. To address this limitation, we use a constant-pH methodology to obtain a true ensemble of multiple protonation states of a titratable system at a given pH and analyze the ensemble using the Bennett acceptance ratio (BAR) method. The constant pH method is based on the combination of enveloping distribution sampling (EDS) with the Hamiltonian replica exchange method (HREM), which yields an accurate semi-grand canonical ensemble of a titratable system. By considering the free energy change of constraining multiple protonation states to a single state or releasing a single protonation state to multiple states, the pH dependent binding free energy profile can be obtained. We perform benchmark simulations of a host-guest system: cucurbit[7]uril (CB[7]) and benzimidazole (BZ). BZ experiences a large pKa shift upon complex formation. The pH-dependent binding free energy profiles of the benchmark system are obtained with three different long-range interaction calculation schemes: a cutoff, the particle mesh Ewald (PME), and the isotropic periodic sum (IPS) method. Our scheme captures the pH-dependent behavior of binding free energy successfully. Absolute binding free energy values obtained with the PME and IPS methods are consistent, while cutoff method results are off by 2 kcal mol−1. We also discuss the characteristics of three long-range interaction calculation methods for constant-pH simulations.

Categories: Journal Articles

North Pacific ‘blob’ stirs up fisheries management

Nature - Wed, 08/19/2015 - 23:00

North Pacific ‘blob’ stirs up fisheries management

Nature 524, 7566 (2015). http://www.nature.com/doifinder/10.1038/nature.2015.18218

Author: Virginia Gewin

Unusually warm ocean strengthens calls to consider ecosystem variables in setting catch limits.

Categories: Journal Articles

Structural characterization of the carbohydrate-binding module of NanA sialidase, a pneumococcal virulence factor

BMC Structural Biology - Wed, 08/19/2015 - 19:00
Background: Streptococcus pneumoniae Neuraminidase A (NanA) is a multi-domain protein anchored to the bacterial surface. Upstream of the catalytic domain of NanA is a domain that conforms to the sialic acid-recognising CBM40 family of the CAZY (carbohydrate-active enzymes) database. This domain has been identified to play a critical role in allowing the bacterium to promote adhesion and invasion of human brain microvascular endothelial cells, and hence may play a key role in promoting bacterial meningitis. In addition, the CBM40 domain has also been reported to activate host chemokines and neutrophil recruitment during infection. Results: Crystal structures of both apo- and holo- forms of the NanA CBM40 domain (residues 121 to 305), have been determined to 1.8 Å resolution. The domain shares the fold of other CBM40 domains that are associated with sialidases. When in complex with α2,3- or α2,6-sialyllactose, the domain is shown to interact only with the terminal sialic acid. Significantly, a deep acidic pocket adjacent to the sialic acid-binding site is identified, which is occupied by a lysine from a symmetry-related molecule in the crystal. This pocket is adjacent to a region that is predicted to be involved in protein-protein interactions. Conclusions: The structural data provide the details of linkage-independent sialyllactose binding by NanA CBM40 and reveal striking surface features that may hold the key to recognition of binding partners on the host cell surface. The structure also suggests that small molecules or sialic acid analogues could be developed to fill the acidic pocket and hence provide a new therapeutic avenue against meningitis caused by S. pneumoniae.
Categories: Journal Articles

Age-related cleavages of crystallins in human lens cortical fiber cells generate a plethora of endogenous peptides and high molecular weight complexes

ABSTRACT

Low molecular weight peptides derived from the breakdown of crystallins have been reported in adult human lenses. The proliferation of these LMW peptides coincides with the earliest stages of cataract formation, suggesting that the protein cleavages involved may contribute to the aggregation and insolubilization of crystallins. This study reports the identification of 238 endogenous LMW crystallin peptides from the cortical extracts of four human lenses representing young, middle and old-age human lenses. Analysis of the peptide terminal amino acids showed that Lys and Arg were situated at the C-terminus with significantly higher frequency compared to other residues, suggesting that trypsin-like proteolysis may be active in the lens cortical fiber cells. Selected reaction monitoring analysis of an endogenous αA-crystallin peptide (αA57-65) showed that the concentration of this peptide in the human lens increased gradually to middle age, after which the rate of αA57-65 formation escalated significantly. Using 2D gel electrophoresis/nanoLC-ESI-MS/MS, 12 protein complexes of 40–150 kDa consisting of multiple crystallin components were characterized from the water soluble cortical extracts of an adult human lens. The detection of these protein complexes suggested the possibility of crystallin cross-linking, with these complexes potentially acting to stabilize degraded crystallins by sequestration into water soluble complexes. Proteins 2015; 83:1878–1886. © 2015 Wiley Periodicals, Inc.

Categories: Journal Articles

Conference tweeting rule frustrates ecologists

Nature - Tue, 08/18/2015 - 23:00

Conference tweeting rule frustrates ecologists

Nature 524, 7566 (2015). doi:10.1038/524391f

Author: Chris Woolston

Complaints ensued when attendees at an ecology meeting were asked to get permission before live-tweeting.

Categories: Journal Articles

Evolution: Gene transfer in complex cells

Nature - Tue, 08/18/2015 - 23:00

Evolution: Gene transfer in complex cells

Nature 524, 7566 (2015). doi:10.1038/nature15205

Authors: John M. Archibald

A comparative genomic study shows that, during evolution, nucleus-containing cells acquired DNA from bacteria primarily by endosymbiosis — the uptake and integration of one cell by another. See Article p.427

Categories: Journal Articles

Endosymbiotic origin and differential loss of eukaryotic genes

Nature - Tue, 08/18/2015 - 23:00

Endosymbiotic origin and differential loss of eukaryotic genes

Nature 524, 7566 (2015). doi:10.1038/nature14963

Authors: Chuan Ku, Shijulal Nelson-Sathi, Mayo Roettger, Filipa L. Sousa, Peter J. Lockhart, David Bryant, Einat Hazkani-Covo, James O. McInerney, Giddy Landan & William F. Martin

Chloroplasts arose from cyanobacteria, mitochondria arose from proteobacteria. Both organelles have conserved their prokaryotic biochemistry, but their genomes are reduced, and most organelle proteins are encoded in the nucleus. Endosymbiotic theory posits that bacterial genes in eukaryotic genomes entered the eukaryotic lineage via organelle ancestors.

Categories: Journal Articles

Sidekick 2 directs formation of a retinal circuit that detects differential motion

Nature - Tue, 08/18/2015 - 23:00

Sidekick 2 directs formation of a retinal circuit that detects differential motion

Nature 524, 7566 (2015). doi:10.1038/nature14682

Authors: Arjun Krishnaswamy, Masahito Yamagata, Xin Duan, Y. Kate Hong & Joshua R. Sanes

In the mammalian retina, processes of approximately 70 types of interneurons form specific synapses on roughly 30 types of retinal ganglion cells (RGCs) in a neuropil called the inner plexiform layer. Each RGC type extracts salient features from visual input, which are sent deeper into the brain for further processing. The specificity and stereotypy of synapses formed in the inner plexiform layer account for the feature-detecting ability of RGCs. Here we analyse the development and function of synapses on one mouse RGC type, called the W3B-RGC. These cells have the remarkable property of responding when the timing of the movement of a small object differs from that of the background, but not when they coincide. Such cells, known as local edge detectors or object motion sensors, can distinguish moving objects from a visual scene that is also moving. We show that W3B-RGCs receive strong and selective input from an unusual excitatory amacrine cell type known as VG3-AC (vesicular glutamate transporter 3). Both W3B-RGCs and VG3-ACs express the immunoglobulin superfamily recognition molecule sidekick 2 (Sdk2), and both loss- and gain-of-function studies indicate that Sdk2-dependent homophilic interactions are necessary for the selectivity of the connection. The Sdk2-specified synapse is essential for visual responses of W3B-RGCs: whereas bipolar cells relay visual input directly to most RGCs, the W3B-RGCs receive much of their input indirectly, via the VG3-ACs. This non-canonical circuit introduces a delay into the pathway from photoreceptors in the centre of the receptive field to W3B-RGCs, which could improve their ability to judge the synchrony of local and global motion.

Categories: Journal Articles

Cell mixing induced by myc is required for competitive tissue invasion and destruction

Nature - Tue, 08/18/2015 - 23:00

Cell mixing induced by myc is required for competitive tissue invasion and destruction

Nature 524, 7566 (2015). doi:10.1038/nature14684

Authors: Romain Levayer, Barbara Hauert & Eduardo Moreno

Cell–cell intercalation is used in several developmental processes to shape the normal body plan. There is no clear evidence that intercalation is involved in pathologies. Here we use the proto-oncogene myc to study a process analogous to early phase of tumour expansion: myc-induced cell competition. Cell competition is a conserved mechanism driving the elimination of slow-proliferating cells (so-called ‘losers’) by faster-proliferating neighbours (so-called ‘winners’) through apoptosis and is important in preventing developmental malformations and maintain tissue fitness. Here we show, using long-term live imaging of myc-driven competition in the Drosophila pupal notum and in the wing imaginal disc, that the probability of elimination of loser cells correlates with the surface of contact shared with winners. As such, modifying loser–winner interface morphology can modulate the strength of competition. We further show that elimination of loser clones requires winner–loser cell mixing through cell–cell intercalation. Cell mixing is driven by differential growth and the high tension at winner–winner interfaces relative to winner–loser and loser–loser interfaces, which leads to a preferential stabilization of winner–loser contacts and reduction of clone compactness over time. Differences in tension are generated by a relative difference in F-actin levels between loser and winner junctions, induced by differential levels of the membrane lipid phosphatidylinositol (3,4,5)-trisphosphate. Our results establish the first link between cell–cell intercalation induced by a proto-oncogene and how it promotes invasiveness and destruction of healthy tissues.

Categories: Journal Articles

Topological and sequence information predict that foldons organize a partially overlapped and hierarchical structure

ABSTRACT

It has been suggested that proteins have substructures, called foldons, which can cooperatively fold into the native structure. However, several prior investigations define foldons in various ways, citing different foldon characteristics, thereby making the concept of a foldon ambiguous. In this study, we perform a Gō model simulation and analyze the characteristics of substructures that cooperatively fold into the native-like structure. Although some results do not agree well with the experimental evidence due to the simplicity of our coarse-grained model, our results strongly suggest that cooperatively folding units sometimes organize a partially overlapped and hierarchical structure. This view makes us easy to interpret some different proposal about the foldon as a difference of the hierarchical structure. On the basis of this finding, we present a new method to assign foldons and their hierarchy, using structural and sequence information. The results show that the foldons assigned by our method correspond to the intermediate structures identified by some experimental techniques. The new method makes it easy to predict whether a protein folds sequentially into the native structure or whether some foldons fold into the native structure in parallel. Proteins 2015; 83:1900–1913. © 2015 Wiley Periodicals, Inc.

Categories: Journal Articles

Protein structure refinement via molecular-dynamics simulations: What works and what does not?

ABSTRACT

Protein structure refinement during CASP11 by the Feig group was described. Molecular dynamics simulations were used in combination with an improved selection and averaging protocol. On average, modest refinement was achieved with some targets improved significantly. Analysis of the CASP submission from our group focused on refinement success versus amount of sampling, refinement of different secondary structure elements and whether refinement varied as a function of which group provided initial models. The refinement of local stereochemical features was examined via the MolProbity score and an updated protocol was developed that can generate high-quality structures with very low MolProbity scores for most starting structures with modest computational effort. Proteins 2015. © 2015 Wiley Periodicals, Inc.

Categories: Journal Articles

Multimetallic catalysed cross-coupling of aryl bromides with aryl triflates

Nature - Sun, 08/16/2015 - 23:00

Multimetallic catalysed cross-coupling of aryl bromides with aryl triflates

Nature 524, 7566 (2015). doi:10.1038/nature14676

Authors: Laura K. G. Ackerman, Matthew M. Lovell & Daniel J. Weix

The advent of transition-metal catalysed strategies for forming new carbon-carbon bonds has revolutionized the field of organic chemistry, enabling the efficient synthesis of ligands, materials, and biologically active molecules. In cases where a single metal fails to promote a selective or efficient transformation, the synergistic cooperation of two distinct catalysts—multimetallic catalysis—can be used instead. Many important reactions rely on multimetallic catalysis, such as the Wacker oxidation of olefins and the Sonogashira coupling of alkynes with aryl halides, but this approach has largely been limited to the use of metals with distinct reactivities, with only one metal catalyst undergoing oxidative addition. Here, we demonstrate that cooperativity between two group 10 metal catalysts—(bipyridine)nickel and (1,3-bis(diphenylphosphino)propane)palladium—enables a general cross-Ullmann reaction (the cross-coupling of two different aryl electrophiles). Our method couples aryl bromides with aryl triflates directly, eliminating the use of arylmetal reagents and avoiding the challenge of differentiating between multiple carbon–hydrogen bonds that is required for direct arylation methods. Selectivity can be achieved without an excess of either substrate and originates from the orthogonal reactivity of the two catalysts and the relative stability of the two arylmetal intermediates. While (1,3-bis(diphenylphosphino)propane)palladium reacts preferentially with aryl triflates to afford a persistent intermediate, (bipyridine)nickel reacts preferentially with aryl bromides to form a transient, reactive intermediate. Although each catalyst forms less than 5 per cent cross-coupled product in isolation, together they are able to achieve a yield of up to 94 per cent. Our results reveal a new method for the synthesis of biaryls, heteroaryls, and dienes, as well as a general mechanism for the selective transfer of ligands between two metal catalysts. We anticipate that this reaction will simplify the synthesis of pharmaceuticals, many of which are currently made with pre-formed organometallic reagents, and lead to the discovery of new multimetallic reactions.

Categories: Journal Articles

CASP11 refinement experiments with ROSETTA

ABSTRACT

We report new Rosetta-based approaches to tackling the major issues that confound protein structure refinement, and the testing of these approaches in the CASP11 experiment. Automated refinement protocols were developed that integrate a range of sampling methods using parallel computation and multiobjective optimization. In CASP11, we used a more aggressive large-scale structure rebuilding approach for poor starting models, and a less aggressive local rebuilding plus core refinement approach for starting models likely to be closer to the native structure. The more incorrectly modeled a structure was predicted to be, the more it was allowed to vary during refinement. The CASP11 experiment revealed strengths and weaknesses of the approaches: the high-resolution strategy incorporating local rebuilding with core refinement consistently improved starting structures, while the low-resolution strategy incorporating the reconstruction of large parts of the structures improved starting models in some cases but often considerably worsened them, largely because of model selection issues. Overall, the results suggest the high-resolution refinement protocol is a promising method orthogonal to other approaches, while the low-resolution refinement method clearly requires further development. Proteins 2015. © 2015 Wiley Periodicals, Inc.

Categories: Journal Articles

Accurate contact predictions using covariation techniques and machine learning

ABSTRACT

Here we present the results of residue–residue contact predictions achieved in CASP11 by the CONSIP2 server, which is based around our MetaPSICOV contact prediction method. On a set of 40 target domains with a median family size of around 40 effective sequences, our server achieved an average top-L/5 long-range contact precision of 27%. MetaPSICOV method bases on a combination of classical contact prediction features, enhanced with three distinct covariation methods embedded in a two-stage neural network predictor. Some unique features of our approach are (1) the tuning between the classical and covariation features depending on the depth of the input alignment and (2) a hybrid approach to generate deepest possible multiple-sequence alignments by combining jackHMMer and HHblits. We discuss the CONSIP2 pipeline, our results and show that where the method underperformed, the major factor was relying on a fixed set of parameters for the initial sequence alignments and not attempting to perform domain splitting as a preprocessing step. Proteins 2015. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Categories: Journal Articles

Molecular dynamics simulations indicate that tyrosineB10 limits motions of distal histidine to regulate CO binding in soybean leghemoglobin

ABSTRACT

Myoglobin (Mb) uses strong electrostatic interaction in its distal heme pocket to regulate ligand binding. The mechanism of regulation of ligand binding in soybean leghemoglobin a (Lba) has been enigmatic and more so due to the absence of gaseous ligand bound atomic resolution three-dimensional structure of the plant globin. While the 20-fold higher oxygen affinity of Lba compared with Mb is required for its dual physiological function, the mechanism by which this high affinity is achieved is only emerging. Extensive mutational analysis combined with kinetic and CO-FT-IR spectroscopic investigation led to the hypothesis that Lba depended on weakened electrostatic interaction between distal HisE7 and bound ligand achieved by invoking B10Tyr, which itself hydrogen bonds with HisE7 thus restricting it in a single conformation detrimental to Mb-like strong electrostatic interaction. Such theory has been re-assessed here using CO-Lba in silico model and molecular dynamics simulation. The investigation supports the presence of at least two major conformations of HisE7 in Lba brought about by imidazole ring flip, one of which makes hydrogen bonds effectively with B10Tyr affecting the former's ability to stabilize bound ligand, while the other does not. However, HisE7 in Lba has limited conformational freedom unlike high frequency of imidazole ring flips observed in Mb and in TyrB10Leu mutant of Lba. Thus, it appears that TyrB10 limits the conformational freedom of distal His in Lba, tuning down ligand dissociation rate constant by reducing the strength of hydrogen bonding to bound ligand, which the freedom of distal His of Mb allows. Proteins 2015; 83:1836–1848. © 2015 Wiley Periodicals, Inc.

Categories: Journal Articles

Crystal structure of the PepSY-containing domain of the YpeB protein involved in germination of bacillus spores

ABSTRACT

The crystal structure of the C-terminal domain of the Bacillus megaterium YpeB protein has been solved by X-ray crystallography to 1.80-Å resolution. The full-length protein is essential in stabilising the SleB cortex lytic enzyme in Bacillus spores, and may have a role in regulating SleB activity during spore germination. The YpeB-C crystal structure comprises three tandemly repeated PepSY domains, which are aligned to form an extended laterally compressed molecule. A predominantly positively charged region located in the second PepSY domain may provide a site for protein interactions that are important in stabilising SleB and YpeB within the spore. Proteins 2015; 83:1914–1921. © 2015 Wiley Periodicals, Inc.

Categories: Journal Articles

New insights into the molecular mechanism of the Rab GTPase Sec4p activation

BMC Structural Biology - Wed, 08/12/2015 - 07:00
Background: Sec4p is a small monomeric Ras-related GTP-binding protein (23 kDa) that regulates polarized exocytosis in S. cerevisiae. In this study we examine the structural effects of a conserved serine residue in the P-loop corresponding to G12 in Ras. Results: We show that the Sec4p residue serine 29 forms a hydrogen bond with the nucleotide. Mutations of this residue have a different impact than equivalent mutations in Ras and can form stable associations with the exchange factor allowing us to elucidate the structure of a complex of Sec4p bound to the exchange factor Sec2p representing an early stage of the exchange reaction. Conclusions: Our structural investigation of the Sec4p-Sec2p complex reveals the role of the Sec2p coiled-coil domain in facilitating the fast kinetics of the exchange reaction. For Ras-family GTPases, single point mutations that impact the signaling state of the molecule have been well described however less structural information is available for equivalent mutations in the case of Rab proteins. Understanding the structural properties of mutants such as the one described here, provides useful insights into unique aspects of Rab GTPase function.
Categories: Journal Articles
Syndicate content