Nature Methods 12, 1058 (2015). doi:10.1038/nmeth.3578

Authors: Tyler Garvin, Robert Aboukhalil, Jude Kendall, Timour Baslan, Gurinder S Atwal, James Hicks, Michael Wigler & Michael C Schatz

We present Ginkgo (http://qb.cshl.edu/ginkgo), a user-friendly, open-source web platform for the analysis of single-cell copy-number variations (CNVs). Ginkgo automatically constructs copy-number profiles of cells from mapped reads and constructs phylogenetic trees of related cells. We validated Ginkgo by reproducing the results of five major studies. After comparing three commonly used single-cell amplification techniques, we concluded that degenerate oligonucleotide-primed PCR is the most consistent for CNV analysis.

Nature Methods 12, 1065 (2015). doi:10.1038/nmeth.3579

Authors: Florian Levet, Eric Hosy, Adel Kechkar, Corey Butler, Anne Beghin, Daniel Choquet & Jean-Baptiste Sibarita

Structure of mammalian eIF3 in the context of the 43S preinitiation complex

Nature 525, 7570 (2015). doi:10.1038/nature14891

Authors: Amedee des Georges, Vidya Dhote, Lauriane Kuhn, Christopher U. T. Hellen, Tatyana V. Pestova, Joachim Frank & Yaser Hashem

During eukaryotic translation initiation, 43S complexes, comprising a 40S ribosomal subunit, initiator transfer RNA and initiation factors (eIF) 2, 3, 1 and 1A, attach to the 5′-terminal region of messenger RNA and scan along it to the initiation codon. Scanning on structured mRNAs also requires

Crystal structures of a double-barrelled fluoride ion channel

Nature 525, 7570 (2015). doi:10.1038/nature14981

Authors: Randy B. Stockbridge, Ludmila Kolmakova-Partensky, Tania Shane, Akiko Koide, Shohei Koide, Christopher Miller & Simon Newstead

To contend with hazards posed by environmental fluoride, microorganisms export this anion through F−-specific ion channels of the Fluc family. Since the recent discovery of Fluc channels, numerous idiosyncratic features of these proteins have been unearthed, including strong selectivity for F− over Cl− and dual-topology dimeric assembly. To understand the chemical basis for F− permeation and how the antiparallel subunits convene to form a F−-selective pore, here we solve the crystal structures of two bacterial Fluc homologues in complex with three different monobody inhibitors, with and without F− present, to a maximum resolution of 2.1 Å. The structures reveal a surprising ‘double-barrelled’ channel architecture in which two F− ion pathways span the membrane, and the dual-topology arrangement includes a centrally coordinated cation, most likely Na+. F− selectivity is proposed to arise from the very narrow pores and an unusual anion coordination that exploits the quadrupolar edges of conserved phenylalanine rings.

Programs that control immune cell homeostasis are orchestrated through the coordinated action of a number of regulatory cell populations, including regulatory T cells, regulatory B cells, myeloid-derived suppressor cells, alternatively-activated macrophages and tolerogenic dendritic cells. These regulatory cell populations can prevent harmful inflammation following completion of protective responses and thwart the development of autoimmune pathology. However, they also have a detrimental role in cancer by favoring escape from immune surveillance.

by Elizabeth Ortiz-Gutiérrez, Karla García-Cruz, Eugenio Azpeitia, Aaron Castillo, María de la Paz Sánchez, Elena R. Álvarez-Buylla

Cell cycle control is fundamental in eukaryotic development. Several modeling efforts have been used to integrate the complex network of interacting molecular components involved in cell cycle dynamics. In this paper, we aimed at recovering the regulatory logic upstream of previously known components of cell cycle control, with the aim of understanding the mechanisms underlying the emergence of the cyclic behavior of such components. We focus on Arabidopsis thaliana, but given that many components of cell cycle regulation are conserved among eukaryotes, when experimental data for this system was not available, we considered experimental results from yeast and animal systems. We are proposing a Boolean gene regulatory network (GRN) that converges into only one robust limit cycle attractor that closely resembles the cyclic behavior of the key cell-cycle molecular components and other regulators considered here. We validate the model by comparing our in silico configurations with data from loss- and gain-of-function mutants, where the endocyclic behavior also was recovered. Additionally, we approximate a continuous model and recovered the temporal periodic expression profiles of the cell-cycle molecular components involved, thus suggesting that the single limit cycle attractor recovered with the Boolean model is not an artifact of its discrete and synchronous nature, but rather an emergent consequence of the inherent characteristics of the regulatory logic proposed here. This dynamical model, hence provides a novel theoretical framework to address cell cycle regulation in plants, and it can also be used to propose novel predictions regarding cell cycle regulation in other eukaryotes.

Abstract

Three implicit solvent models, namely GBMVII, FACTS, and SCPISM, were evaluated for their abilities to emulate an explicit solvent environment by comparing the simulated conformational ensembles, dynamics, and electrostatic interactions of the Src SH2 domain and the Lyn kinase domain. This assessment in terms of structural features in folded proteins expands upon the use of hydration energy as a metric for comparison. All-against-all rms coordinate deviation, average positional fluctuations, and ion-pair distance distribution were used to compare the implicit solvent models with the TIP3P explicit solvent model. Our study shows that the Src SH2 domains solvated with TIP3P, GBMVII, and FACTS sample similar global conformations. Additionally, the Src SH2 ion-pair distance distributions of solvent-exposed side chains corresponding to TIP3P, GBMVII, and FACTS do not differ substantially, indicating that GBMVII and FACTS are capable of modeling these electrostatic interactions. The ion-pair distance distributions of SCPISM are distinct from others, demonstrating that these electrostatic interactions are not adequately reproduced with the SCPISM model. On the other hand, for the Lyn kinase domain, a non-globular protein with bilobal structure and a large concavity on the surface, implicit solvent does not accurately model solvation to faithfully reproduce partially buried electrostatic interactions and lobe-lobe conformations. Our work reveals that local structure and dynamics of small, globular proteins are modeled well using FACTS and GBMVII. Nonetheless, global conformations and electrostatic interactions in concavities of multi-lobal proteins resulting from simulations with implicit solvent models do not match those obtained from explicit water simulations.

Abstract

Salt bridges are essential to protein stability and dynamics. Despite the importance, there has been scarce of detailed discussion on how salt bridge partners interact with each other in distinct solvent exposed environments. In this study, employing a recent generalized orthogonal space tempering (gOST) method, we enabled efficient molecular dynamics simulation of repetitive breaking and reforming of salt bridge structures within a minimalist salt-bridge model, the Asp-Arg dipeptide and thereby were able to map its detailed free energy landscape in aqueous solution. Free energy surface analysis shows that although individually-solvated states are more favorable, salt-bridge states still occupy a noticeable portion of the overall population. Notably, the competing forces, e.g. intercharge attractions that drive the formation of salt bridges and solvation forces that pull the charged groups away from each other, are energetically comparable. As the result, the salt bridge stability is highly tunable by local environments; for instance when local water molecules are perturbed to interact more strongly with each other, the population of the salt-bridge states is likely to increase. Our results reveal the critical role of local solvent structures in modulating salt-bridge partner interactions and imply the importance of water fluctuations on conformational dynamics that involves solvent accessible salt bridge formations.

ABSTRACT

Chemokines form a family of signaling proteins mainly responsible for directing the traffic of leukocytes, where their biological activity can be modulated by their oligomerization state. We characterize the dynamics and thermodynamic stability of monomer and homodimer structures of CXCL7, one of the most abundant platelet chemokines, using experimental methods that include circular dichroism (CD) and nuclear magnetic resonance (NMR) spectroscopy, and computational methods that include the anisotropic network model (ANM), molecular dynamics (MD) simulations and the distance constraint model (DCM). A consistent picture emerges for the effects of dimerization and Cys5-Cys31 and Cys7-Cys47 disulfide bonds formation. The presence of disulfide bonds is not critical for maintaining structural stability in the monomer or dimer, but the monomer is destabilized more than the dimer upon removal of disulfide bonds. Disulfide bonds play a key role in shaping the characteristics of native state dynamics. The combined analysis shows that upon dimerization flexibly correlated motions are induced between the 30s and 50s loop within each monomer and across the dimer interface. Interestingly, the greatest gain in flexibility upon dimerization occurs when both disulfide bonds are present, and the homodimer is least stable relative to its two monomers. These results suggest that the highly conserved disulfide bonds in chemokines facilitate a structural mechanism that is tuned to optimally distinguish functional characteristics between monomer and dimer. Proteins 2015; 83:1987–2007. © 2015 Wiley Periodicals, Inc.

Our PLOS Comp Biol paper is finally out, here .

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Author(s): M.S. Fernández, J.I. Arias, A. Neira-Carrillo, J.L. Arias

Comparative analyzes of biomineralization models have being crucial for the understanding of the functional properties of biominerals and the elucidation of the processes through which biomacromolecules control the synthesis and structural organization of inorganic mineral-based biomaterials. Among calcium carbonate-containing bioceramics, egg, mollusk and echinoderm shells, and crustacean carapaces, have being fairly well characterized. However, Thoraceca barnacles, although being crustacea, showing molting cycle, build a quite stable and heavily mineralized shell that completely surround the animal, which is for life firmly cemented to the substratum. This makes barnacles an interesting model for studying processes of biomineralization. Here we studied the main microstructural and ultrastructural features of Austromegabalanus psittacus barnacle shell, characterize the occurrence of specific proteoglycans (keratan-, dermatan- and chondroitin-6-sulfate proteoglycans) in different soluble and insoluble organic fractions extracted from the shell, and tested them for their ability to crystallize calcium carbonate in vitro. Our results indicate that, in the barnacle model, proteoglycans are good candidates for the modification of the calcite crystal morphology, although the cooperative effect of some additional proteins in the shell could not be excluded.

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Author(s): Paraskevi L. Tsiolaki, Nikolaos N. Louros, Stavros J. Hamodrakas, Vassiliki A. Iconomidou

Amyloidogenic proteins like human Cystatin C (hCC) have been shown to form dimers and oligomers by exchange of subdomains of the monomeric proteins. Normally, the hCC monomer, a low molecular type 2 Cystatin, consists of 120 amino acid residues and functions as an inhibitor of cysteine proteases. The oligomerization of hCC is involved in the pathophysiology of a rare form of amyloidosis namely Icelandic hereditary cerebral amyloid angiopathy, in which an L68Q mutant is deposited as amyloid in brain arteries of young adults. In order to find the shortest stretch responsible to drive the fibril formation of hCC, we have previously demonstrated that the LQVVR peptide forms amyloid fibrils, in vitro (Tsiolaki et al., 2015). Predictions by AMYLPRED, an amyloidogenic determinant prediction algorithm developed in our lab, led us to synthesize and experimentally study two additional predicted peptides derived from hCC. Along with our previous findings, in this work, we reveal that these peptides self-assemble, in a similar way, into amyloid-like fibrils in vitro, as electron microscopy, X-ray fiber diffraction, ATR FT-IR spectroscopy and Congo red staining studies have shown. Further to our experimental results, all three peptides seem to have a fundamental contribution in forming the “aggregation-prone” core of human Cystatin C.

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Author(s): Gurmeet Kaur, Srikrishna Subramanian

Two putative oxidative-stress sensor proteins from Pseudomonas aeruginosa, PA1607 and PA1374, belong to the MarR family of transcription regulators and possess a unique mode of dimerization. In these proteins, in addition to the α-helices involved in dimerization, inter-subunit contacts are strengthened by additional C-terminal β-strands. Using sequence and structure analysis we show that these β-strands constitute a novel segment-swapped zinc ribbon domain. We detect the presence of the zinc ribbon domain in MarR proteins from many bacterial homologs. While the metal-chelating residues of the zinc ribbons are absent in most members of this family, we could however identify several species of Proteobacteria, Actinobacteria and Firmicutes that possess intact zinc-chelating sites. Conservation pattern of metal-chelating residues together with the extensive structural resemblance to zinc ribbons, in particular to the bridge-region of the dsDNA break repair protein Ku, suggests that the C-terminal β-rich region of these proteins is a zinc ribbon. Sequence analysis also supports a distant evolutionary connection between the zinc ribbons of the MarR and Ku families. However, unlike Ku where the segment-swapped zinc ribbons play a role in DNA-binding and obligate dimerization, their primary role in MarR appears to be in dimerization and strengthening of inter-subunit contacts.

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Author(s): Jochen Wuerges, Lorenzo Caputi, Michele Cianci, Stephane Boivin, Rob Meijers, Stefano Benini

Levansucrases are members of the glycoside hydrolase family and catalyse both the hydrolysis of the substrate sucrose and the transfer of fructosyl units to acceptor molecules. In the presence of sufficient sucrose, this may either lead to the production of fructooligosaccharides or fructose polymers. Aim of this study is to rationalise the differences in the polymerisation properties of bacterial levansucrases and in particular to identify structural features that determine different product spectrum in the levansucrase of the Gram-negative bacterium Erwinia amylovora (Ea Lsc, EC 2.4.1.10) as compared to Gram-positive bacteria such as Bacillus subtilis levansucrase. Ea is an enterobacterial pathogen responsible for the Fire Blight disease in rosaceous plants (e.g., apple and pear) with considerable interest for the agricultural industry. The crystal structure of Ea Lsc was solved at 2.77Å resolution and compared to those of other fructosyltransferases from Gram-positive and Gram-negative bacteria. We propose the structural features, determining the different reaction products, to reside in just a few loops at the rim of the active site funnel. Moreover we propose that loop 8 may have a role in product length determination in Gluconacetobacter diazotrophicus LsdA and Microbacterium saccharophilum FFase. The Ea Lsc structure shows for the first time the products of sucrose hydrolysis still bound in the active site.

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Author(s): Gregory P. Kishchenko, Radostin Danev, Rebecca Fisher, Jie He, Chyongere Hsieh, Michael Marko, Haixin Sui

Zernike phase-plate (ZPP) imaging greatly increases contrast in cryo-electron microscopy, however fringe artifacts appear in the images. A computational de-fringing method has been proposed, but it has not been widely employed, perhaps because the importance of de-fringing has not been clearly demonstrated. For testing purposes, we employed Zernike phase-plate imaging in a cryo-electron tomographic study of radial-spoke complexes attached to microtubule doublets. We found that the contrast enhancement by ZPP imaging made nonlinear denoising insensitive to the filtering parameters, such that simple low-frequency band-pass filtering made the same improvement in map quality. We employed sub-tomogram averaging, which compensates for the effect of the “missing wedge” and considerably improves map quality. We found that fringes (caused by the abrupt cut-on of the central hole in the phase plate) can lead to incorrect representation of a structure that is well-known from the literature. The expected structure was restored by amplitude scaling, as proposed in the literature. Our results show that de-fringing is an important part of image-processing for cryo-electron tomography of macromolecular complexes with ZPP imaging.

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Author(s): Beat Vögeli, Simon Olsson, Roland Riek, Peter Güntert

The study of the spatial sampling of biomolecules is essential to understanding the structure–dynamics–function relationship. We have established a protocol for the determination of multiple-state ensembles based on exact measurements of the nuclear Overhauser effect (eNOE). The protocol is practical since it does not require any additional data, while all other NMR data sets must be supplemented by NOE restraints. The question arises as to how much structural and dynamics information is shared between the eNOEs and other NMR probes. We compile one of the largest and most diverse NMR data sets of a protein to date consisting of eNOEs, RDCs and J couplings for GB3. We show that the eNOEs improve the back-prediction of RDCs and J couplings, either upon use of more than one state, or in comparison to conventional NOEs. Our findings indicate that the eNOE data is self-consistent, consistent with other data, and that the structural representation with multiple states is warranted.

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Author(s): Toshio Moriya, Erman Acar, R. Holland Cheng, Ulla Ruotsalainen

In the single particle reconstruction, the initial 3D structure often suffers from the limited angular sampling artifact. Selecting 2D class averages of particle images generally improves the accuracy and efficiency of the reference-free 3D angle estimation, but causes an insufficient angular sampling to fill the information of the target object in the 3D frequency space. Similarly, the initial 3D structure by the random-conical tilt reconstruction has the well-known “missing cone” artifact. Here, we attempted to solve the limited angular sampling problem by sequentially applying maximum a posteriori estimate with expectation maximization algorithm (sMAP-EM). Using both simulated and experimental cryo-electron microscope images, the sMAP-EM was compared to the direct Fourier method on the basis of reconstruction error and resolution. To establish selection criteria of the final regularization weight for the sMAP-EM, the effects of noise level and sampling sparseness on the reconstructions were examined with evenly distributed sampling simulations. The frequency information filled in the missing cone of the conical tilt sampling simulations was assessed by developing new quantitative measurements. All the results of visual and numerical evaluations showed the sMAP-EM performed better than the direct Fourier method, regardless of the sampling method, noise level, and sampling sparseness. Furthermore, the frequency domain analysis demonstrated that the sMAP-EM can fill the meaningful information in the unmeasured angular space without detailed a priori knowledge of the objects. The current research demonstrated that the sMAP-EM has a high potential to facilitate the determination of 3D protein structures at near atomic-resolution.

Publication date: September 2015
Source:Journal of Structural Biology, Volume 191, Issue 3

Author(s): Jademilson Celestino dos Santos, Amanda Bernardes, Letizia Giampietro, Alessandra Ammazzalorso, Barbara De Filippis, Rosa Amoroso, Igor Polikarpov

Peroxisome Proliferator-Activated Receptors (PPARs) are ligand-dependent transcription factors that control various functions in human organism, including the control of glucose and lipid metabolism. PPARγ is a target of TZD agonists, clinically used to improve insulin sensitivity whereas fibrates, PPARα ligands, lower serum triglyceride levels. We report here the structural studies of GL479, a synthetic dual PPARα/γ agonist, designed by a combination of clofibric acid skeleton and a phenyldiazenyl moiety, as bioisosteric replacement of stilbene group, in complex with both PPARα and PPARγ receptors. GL479 was previously reported as a partial agonist of PPARγ and a full agonist of PPARα with high affinity for both PPARs. Our structural studies reveal different binding modes of GL479 to PPARα and PPARγ, which may explain the distinct activation behaviors observed for each receptor. In both cases the ligand interacts with a Tyr located at helix 12 (H12), resulting in the receptor active conformation. In the complex with PPARα, GL479 occupies the same region of the ligand-binding pocket (LBP) observed for other full agonists, whereas GL479 bound to PPARγ displays a new binding mode. Our results indicate a novel region of PPARs LBP that may be explored for the design of partial agonists as well dual PPARα/γ agonists that combine, simultaneously, the therapeutic effects of the treatment of insulin resistance and dyslipidemia.