Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Master corepressor inactivation through multivalent SLiM-induced polymerization mediated by the oncogene suppressor RAI2.

While the elucidation of regulatory mechanisms of folded proteins is facilitated due to their amenability to high-resolution structural characterization, investigation of these mechanisms in disordered proteins is more challenging due to their structural heterogeneity, which can be captured by a variety of biophysical approaches. Here, we used the transcriptional master corepressor CtBP, which binds the putative metastasis suppressor RAI2 through repetitive SLiMs, as a model system. Using cryo-electron microscopy embedded in an integrative structural biology approach, we show that RAI2 unexpectedly induces CtBP polymerization through filaments of stacked tetrameric CtBP layers. These filaments lead to RAI2-mediated CtBP nuclear foci and relieve its corepressor function in RAI2-expressing cancer cells. The impact of RAI2-mediated CtBP loss-of-function is illustrated by the analysis of a diverse cohort of prostate cancer patients, which reveals a substantial decrease in RAI2 in advanced treatment-resistant cancer subtypes. As RAI2-like SLiM motifs are found in a wide range of organisms, including pathogenic viruses, our findings serve as a paradigm for diverse functional effects through multivalent interaction-mediated polymerization by disordered proteins in healthy and diseased conditions.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Predicted models and CCP4.

In late 2020, the results of CASP14, the 14th event in a series of competitions to assess the latest developments in computational protein structure-prediction methodology, revealed the giant leap forward that had been made by Google's Deepmind in tackling the prediction problem. The level of accuracy in their predictions was the first instance of a competitor achieving a global distance test score of better than 90 across all categories of difficulty. This achievement represents both a challenge and an opportunity for the field of experimental structural biology. For structure determination by macromolecular X-ray crystallography, access to highly accurate structure predictions is of great benefit, particularly when it comes to solving the phase problem. Here, details of new utilities and enhanced applications in the CCP4 suite, designed to allow users to exploit predicted models in determining macromolecular structures from X-ray diffraction data, are presented. The focus is mainly on applications that can be used to solve the phase problem through molecular replacement.

RNA 3D structure modeling by fragment assembly with small-angle X-ray scattering restraints.

SUMMARY: Structure determination is a key step in the functional characterization of many non-coding RNA molecules. High-resolution RNA 3D structure determination efforts, however, are not keeping up with the pace of discovery of new non-coding RNA sequences. This increases the importance of computational approaches and low-resolution experimental data, such as from the small-angle X-ray scattering experiments. We present RNA Masonry, a computer program and a web service for a fully automated modeling of RNA 3D structures. It assemblies RNA fragments into geometrically plausible models that meet user-provided secondary structure constraints, restraints on tertiary contacts, and small-angle X-ray scattering data. We illustrate the method description with detailed benchmarks and its application to structural studies of viral RNAs with SAXS restraints. AVAILABILITY AND IMPLEMENTATION: The program web server is available at http://iimcb.genesilico.pl/rnamasonry. The source code is available at https://gitlab.com/gchojnowski/rnamasonry.

DoubleHelix: nucleic acid sequence identification, assignment and validation tool for cryo-EM and crystal structure models.

Sequence assignment is a key step of the model building process in both cryogenic electron microscopy (cryo-EM) and macromolecular crystallography (MX). If the assignment fails, it can result in difficult to identify errors affecting the interpretation of a model. There are many model validation strategies that help experimentalists in this step of protein model building, but they are virtually non-existent for nucleic acids. Here, I present doubleHelix-a comprehensive method for assignment, identification, and validation of nucleic acid sequences in structures determined using cryo-EM and MX. The method combines a neural network classifier of nucleobase identities and a sequence-independent secondary structure assignment approach. I show that the presented method can successfully assist sequence-assignment step in nucleic-acid model building at lower resolutions, where visual map interpretation is very difficult. Moreover, I present examples of sequence assignment errors detected using doubleHelix in cryo-EM and MX structures of ribosomes deposited in the Protein Data Bank, which escaped the scrutiny of available model-validation approaches. The doubleHelix program source code is available under BSD-3 license at https://gitlab.com/gchojnowski/doublehelix.

Sequence-assignment validation in protein crystal structure models with checkMySequence.

Sequence-register shifts remain one of the most elusive errors in experimental macromolecular models. They may affect model interpretation and propagate to newly built models from older structures. In a recent publication, it was shown that register shifts in cryo-EM models of proteins can be detected using a systematic reassignment of short model fragments to the target sequence. Here, it is shown that the same approach can be used to detect register shifts in crystal structure models using standard, model-bias-corrected electron-density maps (2mFo - DFc). Five register-shift errors in models deposited in the PDB detected using this method are described in detail.

The CCP4 suite: integrative software for macromolecular crystallography.

The Collaborative Computational Project No. 4 (CCP4) is a UK-led international collective with a mission to develop, test, distribute and promote software for macromolecular crystallography. The CCP4 suite is a multiplatform collection of programs brought together by familiar execution routines, a set of common libraries and graphical interfaces. The CCP4 suite has experienced several considerable changes since its last reference article, involving new infrastructure, original programs and graphical interfaces. This article, which is intended as a general literature citation for the use of the CCP4 software suite in structure determination, will guide the reader through such transformations, offering a general overview of the new features and outlining future developments. As such, it aims to highlight the individual programs that comprise the suite and to provide the latest references to them for perusal by crystallographers around the world.

Asymmetric horseshoe-like assembly of peroxisomal yeast oxalyl-CoA synthetase.

Oxalyl-CoA synthetase from Saccharomyces cerevisiae is one of the most abundant peroxisomal proteins in yeast and hence has become a model to study peroxisomal translocation. It contains a C-terminal Peroxisome Targeting Signal 1, which however is partly dispensable, suggesting additional receptor bindings sites. To unravel any additional features that may contribute to its capacity to be recognized as peroxisomal target, we determined its assembly and overall architecture by an integrated structural biology approach, including X-ray crystallography, single particle cryo-electron microscopy and small angle X-ray scattering. Surprisingly, it assembles into mixture of concentration-dependent dimers, tetramers and hexamers by dimer self-association. Hexameric particles form an unprecedented asymmetric horseshoe-like arrangement, which considerably differs from symmetric hexameric assembly found in many other protein structures. A single mutation within the self-association interface is sufficient to abolish any higher-level oligomerization, resulting in a homogenous dimeric assembly. The small C-terminal domain of yeast Oxalyl-CoA synthetase is connected by a partly flexible hinge with the large N-terminal domain, which provides the sole basis for oligomeric assembly. Our data provide a basis to mechanistically study peroxisomal translocation of this target.

AlphaPulldown-a python package for protein-protein interaction screens using AlphaFold-Multimer.

SUMMARY: The artificial intelligence-based structure prediction program AlphaFold-Multimer enabled structural modelling of protein complexes with unprecedented accuracy. Increasingly, AlphaFold-Multimer is also used to discover new protein-protein interactions (PPIs). Here, we present AlphaPulldown, a Python package that streamlines PPI screens and high-throughput modelling of higher-order oligomers using AlphaFold-Multimer. It provides a convenient command-line interface, a variety of confidence scores and a graphical analysis tool. AVAILABILITY AND IMPLEMENTATION: AlphaPulldown is freely available at https://www.embl-hamburg.de/AlphaPulldown. SUPPLEMENTARY INFORMATION: Supplementary note is available at Bioinformatics online.

Using deep-learning predictions of inter-residue distances for model validation.

Determination of protein structures typically entails building a model that satisfies the collected experimental observations and its deposition in the Protein Data Bank. Experimental limitations can lead to unavoidable uncertainties during the process of model building, which result in the introduction of errors into the deposited model. Many metrics are available for model validation, but most are limited to consideration of the physico-chemical aspects of the model or its match to the experimental data. The latest advances in the field of deep learning have enabled the increasingly accurate prediction of inter-residue distances, an advance which has played a pivotal role in the recent improvements observed in the field of protein ab initio modelling. Here, new validation methods are presented based on the use of these precise inter-residue distance predictions, which are compared with the distances observed in the protein model. Sequence-register errors are particularly clearly detected and the register shifts required for their correction can be reliably determined. The method is available in the ConKit package (https://www.conkit.org).

CCP4 Cloud for structure determination and project management in macromolecular crystallography.

Nowadays, progress in the determination of three-dimensional macromolecular structures from diffraction images is achieved partly at the cost of increasing data volumes. This is due to the deployment of modern high-speed, high-resolution detectors, the increased complexity and variety of crystallographic software, the use of extensive databases and high-performance computing. This limits what can be accomplished with personal, offline, computing equipment in terms of both productivity and maintainability. There is also an issue of long-term data maintenance and availability of structure-solution projects as the links between experimental observations and the final results deposited in the PDB. In this article, CCP4 Cloud, a new front-end of the CCP4 software suite, is presented which mitigates these effects by providing an online, cloud-based environment for crystallographic computation. CCP4 Cloud was developed for the efficient delivery of computing power, database services and seamless integration with web resources. It provides a rich graphical user interface that allows project sharing and long-term storage for structure-solution projects, and can be linked to data-producing facilities. The system is distributed with the CCP4 software suite version 7.1 and higher, and an online publicly available instance of CCP4 Cloud is provided by CCP4.

Sequence-assignment validation in cryo-EM models with checkMySequence.

The availability of new artificial intelligence-based protein-structure-prediction tools has radically changed the way that cryo-EM maps are interpreted, but it has not eliminated the challenges of map interpretation faced by a microscopist. Models will continue to be locally rebuilt and refined using interactive tools. This inevitably results in occasional errors, among which register shifts remain one of the most difficult to identify and correct. Here, checkMySequence, a fast, fully automated and parameter-free method for detecting register shifts in protein models built into cryo-EM maps, is introduced. It is shown that the method can assist model building in cases where poorer map resolution hinders visual interpretation. It is also shown that checkMySequence could have helped to avoid a widely discussed sequence-register error in a model of SARS-CoV-2 RNA-dependent RNA polymerase that was originally detected thanks to a visual residue-by-residue inspection by members of the structural biology community. The software is freely available at https://gitlab.com/gchojnowski/checkmysequence.

Crystal structures of the molecular class A ß-lactamase TEM-171 and its complexes with tazobactam.

The resistance of bacteria to ß-lactam antibiotics is primarily caused by the production of ß-lactamases. Here, novel crystal structures of the native ß-lactamase TEM-171 and two complexes with the widely used inhibitor tazobactam are presented, alongside complementary data from UV spectroscopy and fluorescence quenching. The six chemically identical ß-lactamase molecules in the crystallographic asymmetric unit displayed different degrees of disorder. The tazobactam intermediate was covalently bound to the catalytic Ser70 in the trans-enamine configuration. While the conformation of tazobactam in the first complex resembled that in published ß-lactamase-tazobactam structures, in the second complex, which was obtained after longer soaking of the native crystals in the inhibitor solution, a new and previously unreported tazobactam conformation was observed. It is proposed that the two complexes correspond to different stages along the deacylation path of the acyl-enzyme intermediate. The results provide a novel structural basis for the rational design of new ß-lactamase inhibitors.

findMySequence: a neural-network-based approach for identification of unknown proteins in X-ray crystallography and cryo-EM.

Although experimental protein-structure determination usually targets known proteins, chains of unknown sequence are often encountered. They can be purified from natural sources, appear as an unexpected fragment of a well characterized protein or appear as a contaminant. Regardless of the source of the problem, the unknown protein always requires characterization. Here, an automated pipeline is presented for the identification of protein sequences from cryo-EM reconstructions and crystallographic data. The method's application to characterize the crystal structure of an unknown protein purified from a snake venom is presented. It is also shown that the approach can be successfully applied to the identification of protein sequences and validation of sequence assignments in cryo-EM protein structures.

Cryo-EM and artificial intelligence visualize endogenous protein community members.

Cellular function is underlined by megadalton assemblies organizing in proximity, forming communities. Metabolons are protein communities involving metabolic pathways such as protein, fatty acid, and thioesters of coenzyme-A synthesis. Metabolons are highly heterogeneous due to their function, making their analysis particularly challenging. Here, we simultaneously characterize metabolon-embedded architectures of a 60S pre-ribosome, fatty acid synthase, and pyruvate/oxoglutarate dehydrogenase complex E2 cores de novo. Cryo-electron microscopy (cryo-EM) 3D reconstructions are resolved at 3.84-4.52 Å resolution by collecting <3,000 micrographs of a single cellular fraction. After combining cryo-EM with artificial intelligence-based atomic modeling and de novo sequence identification methods, at this resolution range, polypeptide hydrogen bonding patterns are discernible. Residing molecular components resemble their purified counterparts from other eukaryotes but also exhibit substantial conformational variation with potential functional implications. Our results propose an integrated tool, boosted by machine learning, that opens doors for structural systems biology spearheaded by cryo-EM characterization of native cell extracts.

Structure of the mycobacterial ESX-5 type VII secretion system pore complex.

The ESX-5 type VII secretion system is a membrane-spanning protein complex key to the virulence of mycobacterial pathogens. However, the overall architecture of the fully assembled translocation machinery and the composition of the central secretion pore have remained unknown. Here, we present the high-resolution structure of the 2.1-megadalton ESX-5 core complex. Our structure captured a dynamic, secretion-competent conformation of the pore within a well-defined transmembrane section, sandwiched between two flexible protein layers at the cytosolic entrance and the periplasmic exit. We propose that this flexibility endows the ESX-5 machinery with large conformational plasticity required to accommodate targeted protein secretion. Compared to known secretion systems, a highly dynamic state of the pore may represent a fundamental principle of bacterial secretion machineries.

Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge.

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1 Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

The accuracy of protein models automatically built into cryo-EM maps with ARP/wARP.

Recent developments in cryogenic electron microscopy (cryo-EM) have enabled structural studies of large macromolecular complexes at resolutions previously only attainable using macromolecular crystallography. Although a number of methods can already assist in de novo building of models into high-resolution cryo-EM maps, automated and reliable map interpretation remains a challenge. Presented here is a systematic study of the accuracy of models built into cryo-EM maps using ARP/wARP. It is demonstrated that the local resolution is a good indicator of map interpretability, and for the majority of the test cases ARP/wARP correctly builds 90% of main-chain fragments in regions where the local resolution is 4.0 Šor better. It is also demonstrated that the coordinate accuracy for models built into cryo-EM maps is comparable to that of X-ray crystallographic models at similar local cryo-EM and crystallographic resolutions. The model accuracy also correlates with the refined atomic displacement parameters.

The use of local structural similarity of distant homologues for crystallographic model building from a molecular-replacement solution.

The performance of automated protein model building usually decreases with resolution, mainly owing to the lower information content of the experimental data. This calls for a more elaborate use of the available structural information about macromolecules. Here, a new method is presented that uses structural homologues to improve the quality of protein models automatically constructed using ARP/wARP. The method uses local structural similarity between deposited models and the model being built, and results in longer main-chain fragments that in turn can be more reliably docked to the protein sequence. The application of the homology-based model extension method to the example of a CFA synthase at 2.7 Šresolution resulted in a more complete model with almost all of the residues correctly built and docked to the sequence. The method was also evaluated on 1493 molecular-replacement solutions at a resolution of 4.0 Šand better that were submitted to the ARP/wARP web service for model building. A significant improvement in the completeness and sequence coverage of the built models has been observed.