Higher resolution in cryo-EM by the combination of macromolecular prior knowledge and image-processing tools.

Single-particle cryo-electron microscopy has become a powerful technique for the 3D structure determination of biological molecules. The last decade has seen an astonishing development of both hardware and software, and an exponential growth of new structures obtained at medium-high resolution. However, the knowledge accumulated in this field over the years has hardly been utilized as feedback in the reconstruction of new structures. In this context, this article explores the use of the deep-learning approach deepEMhancer as a regularizer in the RELION refinement process. deepEMhancer introduces prior information derived from macromolecular structures, and contributes to noise reduction and signal enhancement, as well as a higher degree of isotropy. These features have a direct effect on image alignment and reduction of overfitting during iterative refinement. The advantages of this combination are demonstrated for several membrane proteins, for which it is especially useful because of their high disorder and flexibility.

Advances in Xmipp for Cryo-Electron Microscopy: From Xmipp to Scipion.

Xmipp is an open-source software package consisting of multiple programs for processing data originating from electron microscopy and electron tomography, designed and managed by the Biocomputing Unit of the Spanish National Center for Biotechnology, although with contributions from many other developers over the world. During its 25 years of existence, Xmipp underwent multiple changes and updates. While there were many publications related to new programs and functionality added to Xmipp, there is no single publication on the Xmipp as a package since 2013. In this article, we give an overview of the changes and new work since 2013, describe technologies and techniques used during the development, and take a peek at the future of the package.

Structural basis for cooperativity of human monoclonal antibodies to meningococcal factor H-binding protein.

Monoclonal antibody (mAb) cooperativity is a phenomenon triggered when mAbs couples promote increased bactericidal killing compared to individual partners. Cooperativity has been deeply investigated among mAbs elicited by factor H-binding protein (fHbp), a Neisseria meningitidis surface-exposed lipoprotein and one of the key antigens included in both serogroup B meningococcus vaccine Bexsero and Trumenba. Here we report the structural and functional characterization of two cooperative mAbs pairs isolated from Bexsero vaccines. The 3D electron microscopy structures of the human mAb-fHbp-mAb cooperative complexes indicate that the angle formed between the antigen binding fragments (fAbs) assume regular angle and that fHbp is able to bind simultaneously and stably the cooperative mAbs pairs and human factor H (fH) in vitro. These findings shed light on molecular basis of the antibody-based mechanism of protection driven by simultaneous recognition of the different epitopes of the fHbp and underline that cooperativity is crucial in vaccine efficacy.

Validation of electron microscopy initial models via small angle X-ray scattering curves.

MOTIVATION: Cryo electron microscopy (EM) is currently one of the main tools to reveal the structural information of biological macromolecules. The re-construction of three-dimensional (3D) maps is typically carried out following an iterative process that requires an initial estimation of the 3D map to be refined in subsequent steps. Therefore, its determination is key in the quality of the final results, and there are cases in which it is still an open issue in single particle analysis (SPA). Small angle X-ray scattering (SAXS) is a well-known technique applied to structural biology. It is useful from small nanostructures up to macromolecular ensembles for its ability to obtain low resolution information of the biological sample measuring its X-ray scattering curve. These curves, together with further analysis, are able to yield information on the sizes, shapes and structures of the analyzed particles. RESULTS: In this paper, we show how the low resolution structural information revealed by SAXS is very useful for the validation of EM initial 3D models in SPA, helping the following refinement process to obtain more accurate 3D structures. For this purpose, we approximate the initial map by pseudo-atoms and predict the SAXS curve expected for this pseudo-atomic structure. The match between the predicted and experimental SAXS curves is considered as a good sign of the correctness of the EM initial map. AVAILABILITY AND IMPLEMENTATION: The algorithm is freely available as part of the Scipion 1.2 software at http://scipion.i2pc.es/.

Gene signature associated with benign neurofibroma transformation to malignant peripheral nerve sheath tumors.

Benign neurofibromas, the main phenotypic manifestations of the rare neurological disorder neurofibromatosis type 1, degenerate to malignant tumors associated to poor prognosis in about 10% of patients. Despite efforts in the field of (epi)genomics, the lack of prognostic biomarkers with which to predict disease evolution frustrates the adoption of appropriate early therapeutic measures. To identify potential biomarkers of malignant neurofibroma transformation, we integrated four human experimental studies and one for mouse, using a gene score-based meta-analysis method, from which we obtained a score-ranked signature of 579 genes. Genes with the highest absolute scores were classified as promising disease biomarkers. By grouping genes with similar neurofibromatosis-related profiles, we derived panels of potential biomarkers. The addition of promoter methylation data to gene profiles indicated a panel of genes probably silenced by hypermethylation. To identify possible therapeutic treatments, we used the gene signature to query drug expression databases. Trichostatin A and other histone deacetylase inhibitors, as well as cantharidin and tamoxifen, were retrieved as putative therapeutic means to reverse the aberrant regulation that drives to malignant cell proliferation and metastasis. This in silico prediction corroborated reported experimental results that suggested the inclusion of these compounds in clinical trials. This experimental validation supported the suitability of the meta-analysis method used to integrate several sources of public genomic information, and the reliability of the gene signature associated to the malignant evolution of neurofibromas to generate working hypotheses for prognostic and drug-responsive biomarkers or therapeutic measures, thus showing the potential of this in silico approach for biomarker discovery.

XTEND: Extending the depth of field in cryo soft X-ray tomography.

We have developed a new data collection method and processing framework in full field cryo soft X-ray tomography to computationally extend the depth of field (DOF) of a Fresnel zone plate lens. Structural features of 3D-reconstructed eukaryotic cells that are affected by DOF artifacts in standard reconstruction are now recovered. This approach, based on focal series projections, is easily applicable with closed expressions to select specific data acquisition parameters.

Electron Microscopy Structural Insights into CPAP Oligomeric Behavior: A Plausible Assembly Process of a Supramolecular Scaffold of the Centrosome.

Centrosomal P4.1-associated protein (CPAP) is a cell cycle regulated protein fundamental for centrosome assembly and centriole elongation. In humans, the region between residues 897-1338 of CPAP mediates interactions with other proteins and includes a homodimerization domain. CPAP mutations cause primary autosomal recessive microcephaly and Seckel syndrome. Despite of the biological/clinical relevance of CPAP, its mechanistic behavior remains unclear and its C-terminus (the G-box/TCP domain) is the only part whose structure has been solved. This situation is perhaps due in part to the challenges that represent obtaining the protein in a soluble, homogeneous state for structural studies. Our work constitutes a systematic structural analysis on multiple oligomers of HsCPAP897-1338, using single-particle electron microscopy (EM) of negatively stained (NS) samples. Based on image classification into clearly different regular 3D maps (putatively corresponding to dimers and tetramers) and direct observation of individual images representing other complexes of HsCPAP897-1338 (i.e., putative flexible monomers and higher-order multimers), we report a dynamic oligomeric behavior of this protein, where different homo-oligomers coexist in variable proportions. We propose that dimerization of the putative homodimer forms a putative tetramer which could be the structural unit for the scaffold that either tethers the pericentriolar material to centrioles or promotes procentriole elongation. A coarse fitting of atomic models into the NS 3D maps at resolutions around 20 Å is performed only to complement our experimental data, allowing us to hypothesize on the oligomeric composition of the different complexes. In this way, the current EM work represents an initial step toward the structural characterization of different oligomers of CPAP, suggesting further insights to understand how this protein works, contributing to the elucidation of control mechanisms for centriole biogenesis.

Characterization of transfer function, resolution and depth of field of a soft X-ray microscope applied to tomography enhancement by Wiener deconvolution.

Full field soft X-ray microscopy is becoming a powerful imaging technique to analyze whole cells preserved under cryo conditions. Images obtained in these X-ray microscopes can be combined by tomographic reconstruction to quantitatively estimate the three-dimensional (3D) distribution of absorption coefficients inside the cell. The impulse response of an imaging system is one of the factors that limits the quality of the X-ray microscope reconstructions. The main goal of this work is to experimentally measure the 3D impulse response and to assess the optical resolution and depth of field of the Mistral microscope at ALBA synchrotron (Barcelona, Spain). To this end we measure the microscope apparent transfer function (ATF) and we use it to design a deblurring Wiener filter, obtaining an increase in the image quality when applied to experimental datasets collected at ALBA.

The Electron Microscopy eXchange (EMX) initiative.

Three-dimensional electron microscopy (3DEM) of ice-embedded samples allows the structural analysis of large biological macromolecules close to their native state. Different techniques have been developed during the last forty years to process cryo-electron microscopy (cryo-EM) data. Not surprisingly, success in analysis and interpretation is highly correlated with the continuous development of image processing packages. The field has matured to the point where further progress in data and methods sharing depends on an agreement between the packages on how to describe common image processing tasks. Such standardization will facilitate the use of software as well as seamless collaboration, allowing the sharing of rich information between different platforms. Our aim here is to describe the Electron Microscopy eXchange (EMX) initiative, launched at the 2012 Instruct Image Processing Center Developer Workshop, with the intention of developing a first set of standard conventions for the interchange of information for single-particle analysis (EMX version 1.0). These conventions cover the specification of the metadata for micrograph and particle images, including contrast transfer function (CTF) parameters and particle orientations. EMX v1.0 has already been implemented in the Bsoft, EMAN, Xmipp and Scipion image processing packages. It has been and will be used in the CTF and EMDataBank Validation Challenges respectively. It is also being used in EMPIAR, the Electron Microscopy Pilot Image Archive, which stores raw image data related to the 3DEM reconstructions in EMDB.

Structures of Adenovirus Incomplete Particles Clarify Capsid Architecture and Show Maturation Changes of Packaging Protein L1 52/55k.

UNLABELLED: Adenovirus is one of the most complex icosahedral, nonenveloped viruses. Even after its structure was solved at near-atomic resolution by both cryo-electron microscopy and X-ray crystallography, the location of minor coat proteins is still a subject of debate. The elaborated capsid architecture is the product of a correspondingly complex assembly process, about which many aspects remain unknown. Genome encapsidation involves the concerted action of five virus proteins, and proteolytic processing by the virus protease is needed to prime the virion for sequential uncoating. Protein L1 52/55k is required for packaging, and multiple cleavages by the maturation protease facilitate its release from the nascent virion. Light-density particles are routinely produced in adenovirus infections and are thought to represent assembly intermediates. Here, we present the molecular and structural characterization of two different types of human adenovirus light particles produced by a mutant with delayed packaging. We show that these particles lack core polypeptide V but do not lack the density corresponding to this protein in the X-ray structure, thereby adding support to the adenovirus cryo-electron microscopy model. The two types of light particles present different degrees of proteolytic processing. Their structures provide the first glimpse of the organization of L1 52/55k protein inside the capsid shell and of how this organization changes upon partial maturation. Immature, full-length L1 52/55k is poised beneath the vertices to engage the virus genome. Upon proteolytic processing, L1 52/55k disengages from the capsid shell, facilitating genome release during uncoating. IMPORTANCE: Adenoviruses have been extensively characterized as experimental systems in molecular biology, as human pathogens, and as therapeutic vectors. However, a clear picture of many aspects of their basic biology is still lacking. Two of these aspects are the location of minor coat proteins in the capsid and the molecular details of capsid assembly. Here, we provide evidence supporting one of the two current models for capsid architecture. We also show for the first time the location of the packaging protein L1 52/55k in particles lacking the virus genome and how this location changes during maturation. Our results contribute to clarifying standing questions in adenovirus capsid architecture and provide new details on the role of L1 52/55k protein in assembly.

Measurement of the modulation transfer function of an X-ray microscope based on multiple Fourier orders analysis of a Siemens star.

Soft X-ray tomography (SXT) is becoming a powerful imaging technique to analyze eukaryotic whole cells close to their native state. Central to the analysis of the quality of SXT 3D reconstruction is the estimation of the spatial resolution and Depth of Field of the X-ray microscope. In turn, the characterization of the Modulation Transfer Function (MTF) of the optical system is key to calculate both parameters. Consequently, in this work we introduce a fully automated technique to accurately estimate the transfer function of such an optical system. Our proposal is based on the preprocessing of the experimental images to obtain an estimate of the input pattern, followed by the analysis in Fourier space of multiple orders of a Siemens Star test sample, extending in this way its measured frequency range.

CTF Challenge: Result summary.

Image formation in bright field electron microscopy can be described with the help of the contrast transfer function (CTF). In this work the authors describe the "CTF Estimation Challenge", called by the Madrid Instruct Image Processing Center (I2PC) in collaboration with the National Center for Macromolecular Imaging (NCMI) at Houston. Correcting for the effects of the CTF requires accurate knowledge of the CTF parameters, but these have often been difficult to determine. In this challenge, researchers have had the opportunity to test their ability in estimating some of the key parameters of the electron microscope CTF on a large micrograph data set produced by well-known laboratories on a wide set of experimental conditions. This work presents the first analysis of the results of the CTF Estimation Challenge, including an assessment of the performance of the different software packages under different conditions, so as to identify those areas of research where further developments would be desirable in order to achieve high-resolution structural information.

Efficient initial volume determination from electron microscopy images of single particles.

MOTIVATION: Structural information of macromolecular complexes provides key insights into the way they carry out their biological functions. The reconstruction process leading to the final 3D map requires an approximate initial model. Generation of an initial model is still an open and challenging problem in single-particle analysis. RESULTS: We present a fast and efficient approach to obtain a reliable, low-resolution estimation of the 3D structure of a macromolecule, without any a priori knowledge, addressing the well-known issue of initial volume estimation in the field of single-particle analysis. The input of the algorithm is a set of class average images obtained from individual projections of a biological object at random and unknown orientations by transmission electron microscopy micrographs. The proposed method is based on an initial non-lineal dimensionality reduction approach, which allows to automatically selecting representative small sets of class average images capturing the most of the structural information of the particle under study. These reduced sets are then used to generate volumes from random orientation assignments. The best volume is determined from these guesses using a random sample consensus (RANSAC) approach. We have tested our proposed algorithm, which we will term 3D-RANSAC, with simulated and experimental data, obtaining satisfactory results under the low signal-to-noise conditions typical of cryo-electron microscopy. AVAILABILITY: The algorithm is freely available as part of the Xmipp 3.1 package [http://xmipp.cnb.csic.es]. CONTACT: jvargas@cnb.csic.es SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structure and non-structure of centrosomal proteins.

Here we perform a large-scale study of the structural properties and the expression of proteins that constitute the human Centrosome. Centrosomal proteins tend to be larger than generic human proteins (control set), since their genes contain in average more exons (20.3 versus 14.6). They are rich in predicted disordered regions, which cover 57% of their length, compared to 39% in the general human proteome. They also contain several regions that are dually predicted to be disordered and coiled-coil at the same time: 55 proteins (15%) contain disordered and coiled-coil fragments that cover more than 20% of their length. Helices prevail over strands in regions homologous to known structures (47% predicted helical residues against 17% predicted as strands), and even more in the whole centrosomal proteome (52% against 7%), while for control human proteins 34.5% of the residues are predicted as helical and 12.8% are predicted as strands. This difference is mainly due to residues predicted as disordered and helical (30% in centrosomal and 9.4% in control proteins), which may correspond to alpha-helix forming molecular recognition features (α-MoRFs). We performed expression assays for 120 full-length centrosomal proteins and 72 domain constructs that we have predicted to be globular. These full-length proteins are often insoluble: Only 39 out of 120 expressed proteins (32%) and 19 out of 72 domains (26%) were soluble. We built or retrieved structural models for 277 out of 361 human proteins whose centrosomal localization has been experimentally verified. We could not find any suitable structural template with more than 20% sequence identity for 84 centrosomal proteins (23%), for which around 74% of the residues are predicted to be disordered or coiled-coils. The three-dimensional models that we built are available at http://ub.cbm.uam.es/centrosome/models/index.php.

Image formation in cellular X-ray microscopy.

Soft X-ray Tomographic (TomoX) microscopy has become a reality in the last years. The resolution range of this technique nicely fits between confocal and electron microscopies and will play a key role in the elucidation of the organization between the molecular and the organelle levels. In fact, it offers the possibility of imaging three-dimensional structures of hydrated biological specimens near their native state without chemical pre-treatment. Ideally, TomoX reconstructs the specimen absorption coefficients from projections of this specimen, but, unfortunately, X-ray micrographs are only an approximation to projections of the specimen, resulting in inaccuracies if a tomographic reconstruction is performed without explicitly incorporating these approximations. In an attempt to mitigate some of these inaccuracies, we develop in this work an image formation model within the approximation of assuming incoherent illumination.

PRIDEViewer: a novel user-friendly interface to visualize PRIDE XML files.

Current standardization initiatives have greatly contributed to share the information derived by proteomics experiments. One of these initiatives is the XML-based repository PRIDE (PRoteomics IDEntification database), although an XML-based document does not appear to present a user-friendly view at the first glance. PRIDEViewer is a novel Java-based application that presents the information available in a PRIDE XML file in a user-friendly manner, facilitating the interaction among end users as well as the understanding and evaluation of the compiled information. PRIDEViewer is freely available at: http://proteo.cnb.csic.es/prideviewer/.

MARQ: an online tool to mine GEO for experiments with similar or opposite gene expression signatures.

The enormous amount of data available in public gene expression repositories such as Gene Expression Omnibus (GEO) offers an inestimable resource to explore gene expression programs across several organisms and conditions. This information can be used to discover experiments that induce similar or opposite gene expression patterns to a given query, which in turn may lead to the discovery of new relationships among diseases, drugs or pathways, as well as the generation of new hypotheses. In this work, we present MARQ, a web-based application that allows researchers to compare a query set of genes, e.g. a set of over- and under-expressed genes, against a signature database built from GEO datasets for different organisms and platforms. MARQ offers an easy-to-use and integrated environment to mine GEO, in order to identify conditions that induce similar or opposite gene expression patterns to a given experimental condition. MARQ also includes additional functionalities for the exploration of the results, including a meta-analysis pipeline to find genes that are differentially expressed across different experiments. The application is freely available at http://marq.dacya.ucm.es.

Assessment of protein set coherence using functional annotations.

BACKGROUND: Analysis of large-scale experimental datasets frequently produces one or more sets of proteins that are subsequently mined for functional interpretation and validation. To this end, a number of computational methods have been devised that rely on the analysis of functional annotations. Although current methods provide valuable information (e.g. significantly enriched annotations, pairwise functional similarities), they do not specifically measure the degree of homogeneity of a protein set. RESULTS: In this work we present a method that scores the degree of functional homogeneity, or coherence, of a set of proteins on the basis of the global similarity of their functional annotations. The method uses statistical hypothesis testing to assess the significance of the set in the context of the functional space of a reference set. As such, it can be used as a first step in the validation of sets expected to be homogeneous prior to further functional interpretation. CONCLUSION: We evaluate our method by analysing known biologically relevant sets as well as random ones. The known relevant sets comprise macromolecular complexes, cellular components and pathways described for Saccharomyces cerevisiae, which are mostly significantly coherent. Finally, we illustrate the usefulness of our approach for validating 'functional modules' obtained from computational analysis of protein-protein interaction networks. Matlab code and supplementary data are available at http://www.cnb.csic.es/~monica/coherence/

Quasi-atomic model of bacteriophage t7 procapsid shell: insights into the structure and evolution of a basic fold.

The existence of similar folds among major structural subunits of viral capsids has shown unexpected evolutionary relationships suggesting common origins irrespective of the capsids' host life domain. Tailed bacteriophages are emerging as one such family, and we have studied the possible existence of the HK97-like fold in bacteriophage T7. The procapsid structure at approximately 10 A resolution was used to obtain a quasi-atomic model by fitting a homology model of the T7 capsid protein gp10 that was based on the atomic structure of the HK97 capsid protein. A number of fold similarities, such as the fitting of domains A and P into the L-shaped procapsid subunit, are evident between both viral systems. A different feature is related to the presence of the amino-terminal domain of gp10 found at the inner surface of the capsid that might play an important role in the interaction of capsid and scaffolding proteins.

GENECODIS: a web-based tool for finding significant concurrent annotations in gene lists.

We present GENECODIS, a web-based tool that integrates different sources of information to search for annotations that frequently co-occur in a set of genes and rank them by statistical significance. The analysis of concurrent annotations provides significant information for the biologic interpretation of high-throughput experiments and may outperform the results of standard methods for the functional analysis of gene lists. GENECODIS is publicly available at http://genecodis.dacya.ucm.es/.