Scipion3: A workflow engine for cryo-electron microscopy image processing and structural biology.

Image-processing pipelines require the design of complex workflows combining many different steps that bring the raw acquired data to a final result with biological meaning. In the image-processing domain of cryo-electron microscopy single-particle analysis (cryo-EM SPA), hundreds of steps must be performed to obtain the three-dimensional structure of a biological macromolecule by integrating data spread over thousands of micrographs containing millions of copies of allegedly the same macromolecule. The execution of such complicated workflows demands a specific tool to keep track of all these steps performed. Additionally, due to the extremely low signal-to-noise ratio (SNR), the estimation of any image parameter is heavily affected by noise resulting in a significant fraction of incorrect estimates. Although low SNR and processing millions of images by hundreds of sequential steps requiring substantial computational resources are specific to cryo-EM, these characteristics may be shared by other biological imaging domains. Here, we present Scipion, a Python generic open-source workflow engine specifically adapted for image processing. Its main characteristics are: (a) interoperability, (b) smart object model, (c) gluing operations, (d) comparison operations, (e) wide set of domain-specific operations, (f) execution in streaming, (g) smooth integration in high-performance computing environments, (h) execution with and without graphical capabilities, (i) flexible visualization, (j) user authentication and private access to private data, (k) scripting capabilities, (l) high performance, (m) traceability, (n) reproducibility, (o) self-reporting, (p) reusability, (q) extensibility, (r) software updates, and (s) non-restrictive software licensing.

3DCONS-DB: A Database of Position-Specific Scoring Matrices in Protein Structures.

Many studies have used position-specific scoring matrices (PSSM) profiles to characterize residues in protein structures and to predict a broad range of protein features. Moreover, PSSM profiles of Protein Data Bank (PDB) entries have been recalculated in many works for different purposes. Although the computational cost of calculating a single PSSM profile is affordable, many statistical studies or machine learning-based methods used thousands of profiles to achieve their goals, thereby leading to a substantial increase of the computational cost. In this work we present a new database compiling PSSM profiles for the proteins of the PDB. Currently, the database contains 333,532 protein chain profiles involving 123,135 different PDB entries.

Pea PSII-LHCII supercomplexes form pairs by making connections across the stromal gap.

In higher plant thylakoids, the heterogeneous distribution of photosynthetic protein complexes is a determinant for the formation of grana, stacks of membrane discs that are densely populated with Photosystem II (PSII) and its light harvesting complex (LHCII). PSII associates with LHCII to form the PSII-LHCII supercomplex, a crucial component for solar energy conversion. Here, we report a biochemical, structural and functional characterization of pairs of PSII-LHCII supercomplexes, which were isolated under physiologically-relevant cation concentrations. Using single-particle cryo-electron microscopy, we determined the three-dimensional structure of paired C2S2M PSII-LHCII supercomplexes at 14 Å resolution. The two supercomplexes interact on their stromal sides through a specific overlap between apposing LHCII trimers and via physical connections that span the stromal gap, one of which is likely formed by interactions between the N-terminal loops of two Lhcb4 monomeric LHCII subunits. Fast chlorophyll fluorescence induction analysis showed that paired PSII-LHCII supercomplexes are energetically coupled. Molecular dynamics simulations revealed that additional flexible physical connections may form between the apposing LHCII trimers of paired PSII-LHCII supercomplexes in appressed thylakoid membranes. Our findings provide new insights into how interactions between pairs of PSII-LHCII supercomplexes can link adjacent thylakoids to mediate the stacking of grana membranes.

Cryo-soft X-ray tomography as a quantitative three-dimensional tool to model nanoparticle:cell interaction.

BACKGROUND: Recent advances in nanoparticle design have generated new possibilities for nano-biotechnology and nano-medicine. Here we used cryo-soft X-ray tomography (cryo-SXT) to collect comprehensive three-dimensional (3D) data to characterise the interaction of superparamagnetic iron oxide nanoparticles (SPION) with a breast cancer cell line. RESULTS: We incubated MCF-7 (a human breast cancer cell line) from 0 to 24 h with SPION (15 nm average diameter, coated with dimercaptosuccinic acid), a system that has been studied previously using various microscopy and bulk techniques. This system facilitates the validation and contextualization of the new 3D data acquired using the cryo-SXT-based approach. After vitrification, samples tested by correlative cryo-epifluorescent microscopy showed SPION accumulation in acidic vesicles related to the endocytic pathway. Microscopy grids bearing MCF-7 cells were then analysed by cryo-SXT to generate whole cell volume 3D maps. Cryo-SXT is an emerging technique that benefits from high X-ray penetration into the biological material to image close-to-native vitrified cells at nanometric resolution with no chemical fixation or staining agents. This unique possibility of obtaining 3D information from whole cells allows quantitative statistical analysis of SPION-containing vesicle (SCV) accumulation inside cells, including vesicle number and size, distances between vesicles, and their distance from the nucleus. CONCLUSIONS: Correlation between fluorescent microscopy, cryo-SXT and transmission electron microscopy allowed us to identify SCV and to generate 3D data for statistical analysis of SPION:cell interaction. This study supports continuous transfer of the internalized SPION from the plasma membrane to an accumulation area near the cell nucleus. Statistical analysis showed SCV increase in number and size concomitant with longer incubation times, and therefore an increase in their accumulated volume within the cell. This cumulative effect expands the accumulation area and cell organelles such as mitochondria are consequently displaced to the periphery. Our 3D cryo-SXT approach demonstrates that a comprehensive quantitative description of SPION:cell interaction is possible, which will serve as a basis for metal-based nanoparticle design and for selection of those best suited for hyperthermia treatment, drug delivery and image diagnosis in nanobiomedicine.

Alignment of direct detection device micrographs using a robust Optical Flow approach.

The introduction of direct detection devices in cryo-EM has shown that specimens present beam-induced motion (BIM). Consequently, in this work, we develop a BIM correction method at the image level, resulting in an integrated image in which the in-plane BIM blurring is compensated prior to particle picking. The methodology is based on a robust Optical Flow (OF) approach that can efficiently correct for local movements in a rapid manner. The OF works particularly well if the BIM pattern presents a substantial degree of local movements, which occurs in our data sets for Falcon II data. However, for those cases in which the BIM pattern corresponds to global movements, we have found it advantageous to first run a global motion correction approach and to subsequently apply OF. Additionally, spatial analysis of the Optical Flow allows for quantitative analysis of the BIM pattern. The software that incorporates the new approach is available in XMIPP (http://xmipp.cnb.csic.es).

Chromothripsis: breakage-fusion-bridge over and over again.

The acquisition of massive but localized chromosome translocations, a phenomenon termed chromothripsis, has received widespread attention since its discovery over a year ago. Until recently, chromothripsis was believed to originate from a single catastrophic event, but the molecular mechanisms leading to this event are yet to be uncovered. Because a thorough interpretation of the data are missing, the phenomenon itself has wrongly acquired the status of a mechanism used to justify many kinds of complex rearrangements. Although the assumption that all translocations in chromothripsis originate from a single event has met with criticism, satisfactory explanations for the intense but localized nature of this phenomenon are still missing. Here, we show why the data used to describe massive catastrophic rearrangements are incompatible with a model comprising a single event only and propose a molecular mechanism in which a combination of known cellular pathways accounts for chromothripsis. Instead of a single traumatic event, the protection of undamaged chromosomes by telomeres can limit repetitive breakage-fusion-bridge events to a single chromosome arm. Ultimately, common properties of chromosomal instability, such as aneuploidy and centromere fission, might establish the complex genetic pattern observed in this genomic state.

Marker-free image registration of electron tomography tilt-series.

BACKGROUND: Tilt series are commonly used in electron tomography as a means of collecting three-dimensional information from two-dimensional projections. A common problem encountered is the projection alignment prior to 3D reconstruction. Current alignment techniques usually employ gold particles or image derived markers to correctly align the images. When these markers are not present, correlation between adjacent views is used to align them. However, sequential pairwise correlation is prone to bias and the resulting alignment is not always optimal. RESULTS: In this paper we introduce an algorithm to find regions of the tilt series which can be tracked within a subseries of the tilt series. These regions act as landmarks allowing the determination of the alignment parameters. We show our results with synthetic data as well as experimental cryo electron tomography. CONCLUSION: Our algorithm is able to correctly align a single-tilt tomographic series without the help of fiducial markers thanks to the detection of thousands of small image patches that can be tracked over a short number of images in the series.