Cryo-EM structure and B-factor refinement with ensemble representation.

Cryo-EM experiments produce images of macromolecular assemblies that are combined to produce three-dimensional density maps. Typically, atomic models of the constituent molecules are fitted into these maps, followed by a density-guided refinement. We introduce TEMPy-ReFF, a method for atomic structure refinement in cryo-EM density maps. Our method represents atomic positions as components of a Gaussian mixture model, utilising their variances as B-factors, which are used to derive an ensemble description. Extensively tested on a substantial dataset of 229 cryo-EM maps from EMDB ranging in resolution from 2.1-4.9 Å with corresponding PDB and CERES atomic models, our results demonstrate that TEMPy-ReFF ensembles provide a superior representation of cryo-EM maps. On a single-model basis, it performs similarly to the CERES re-refinement protocol, although there are cases where it provides a better fit to the map. Furthermore, our method enables the creation of composite maps free of boundary artefacts. TEMPy-ReFF is useful for better interpretation of flexible structures, such as those involving RNA, DNA or ligands.

RIBFIND2: Identifying rigid bodies in protein and nucleic acid structures.

Molecular structures are often fitted into cryo-EM maps by flexible fitting. When this requires large conformational changes, identifying rigid bodies can help optimize the model-map fit. Tools for identifying rigid bodies in protein structures exist, however an equivalent for nucleic acid structures is lacking. With the increase in cryo-EM maps containing RNA and progress in RNA structure prediction, there is a need for such tools. We previously developed RIBFIND, a program for clustering protein secondary structures into rigid bodies. In RIBFIND2, this approach is extended to nucleic acid structures. RIBFIND2 can identify biologically relevant rigid bodies in important groups of complex RNA structures, capturing a wide range of dynamics, including large rigid-body movements. The usefulness of RIBFIND2-assigned rigid bodies in cryo-EM model refinement was demonstrated on three examples, with two conformations each: Group II Intron complexed IEP, Internal Ribosome Entry Site and the Processome, using cryo-EM maps at 2.7-5 Å resolution. A hierarchical refinement approach, performed on progressively smaller sets of RIBFIND2 rigid bodies, was clearly shown to have an advantage over classical all-atom refinement. RIBFIND2 is available via a web server with structure visualization and as a standalone tool.

Linking Gas-Phase and Solution-Phase Protein Unfolding via Mobile Proton Simulations.

Native mass spectrometry coupled to ion mobility (IM-MS) combined with collisional activation (CA) of ions in the gas phase (in vacuo) is an important method for the study of protein unfolding. It has advantages over classical biophysical and structural techniques as it can be used to analyze small volumes of low-concentration heterogeneous mixtures while maintaining solution-like behavior and does not require labeling with fluorescent or other probes. It is unclear, however, whether the unfolding observed during collision activation experiments mirrors solution-phase unfolding. To bridge the gap between in vacuo and in-solution behavior, we use unbiased molecular dynamics (MD) to create in silico models of in vacuo unfolding of a well-studied protein, the N-terminal domain of ribosomal L9 (NTL9) protein. We utilize a mobile proton algorithm (MPA) to create 100 thermally unfolded and coulombically unfolded in silico models for observed charge states of NTL9. The unfolding behavior in silico replicates the behavior in-solution and is in line with the in vacuo observations; however, the theoretical collision cross section (CCS) of the in silico models was lower compared to that of the in vacuo data, which may reflect reduced sampling.

Evaluation of acquisition modes for semi-quantitative analysis by targeted and untargeted mass spectrometry.

RATIONALE: Analyte quantitation by mass spectrometry underpins a diverse range of scientific endeavors. The fast-growing field of mass spectrometer development has resulted in several targeted and untargeted acquisition modes suitable for these applications. By characterizing the acquisition methods available on an ion mobility (IM)-enabled orthogonal acceleration time-of-flight (oa-ToF) instrument, the optimum modes for analyte semi-quantitation can be deduced. METHODS: Serial dilutions of commercial metabolite, peptide, or cross-linked peptide analytes were prepared in matrices of human urine or Escherichia coli digest. Each analyte dilution was introduced into an IM separation-enabled oa-ToF mass spectrometer by reversed-phase liquid chromatography and electrospray ionization. Data were acquired for each sample in duplicate using nine different acquisition modes, including four IM-enabled acquisitions modes, available on the mass spectrometer. RESULTS: Five (metabolite) or seven (peptide/cross-linked peptide) point calibration curves were prepared for analytes across each of the acquisition modes. A nonlinear response was observed at high concentrations for some modes, attributed to saturation effects. Two correction methods, one MS1 isotope-correction and one MS2 ion intensity-correction, were applied to address this observation, resulting in an up to twofold increase in dynamic range. By averaging the semi-quantitative results across analyte classes, two parameters, linear dynamic range (LDR) and lower limit of quantification (LLOQ), were determined to evaluate each mode. CONCLUSION: A comparison of the acquisition modes revealed that data-independent acquisition and parallel reaction monitoring methods are most robust for semi-quantitation when considering achievable LDR and LLOQ. IM-enabled modes exhibited sensitivity increases, but a simultaneous reduction in dynamic range required correction methods to recover. These findings will assist users in identifying the optimum acquisition mode for their analyte quantitation needs, supporting a diverse range of applications and providing guidance for future acquisition mode developments.

Integration of Mass Spectrometry Data for Structural Biology.

Mass spectrometry (MS) is increasingly being used to probe the structure and dynamics of proteins and the complexes they form with other macromolecules. There are now several specialized MS methods, each with unique sample preparation, data acquisition, and data processing protocols. Collectively, these methods are referred to as structural MS and include cross-linking, hydrogen-deuterium exchange, hydroxyl radical footprinting, native, ion mobility, and top-down MS. Each of these provides a unique type of structural information, ranging from composition and stoichiometry through to residue level proximity and solvent accessibility. Structural MS has proved particularly beneficial in studying protein classes for which analysis by classic structural biology techniques proves challenging such as glycosylated or intrinsically disordered proteins. To capture the structural details for a particular system, especially larger multiprotein complexes, more than one structural MS method with other structural and biophysical techniques is often required. Key to integrating these diverse data are computational strategies and software solutions to facilitate this process. We provide a background to the structural MS methods and briefly summarize other structural methods and how these are combined with MS. We then describe current state of the art approaches for the integration of structural MS data for structural biology. We quantify how often these methods are used together and provide examples where such combinations have been fruitful. To illustrate the power of integrative approaches, we discuss progress in solving the structures of the proteasome and the nuclear pore complex. We also discuss how information from structural MS, particularly pertaining to protein dynamics, is not currently utilized in integrative workflows and how such information can provide a more accurate picture of the systems studied. We conclude by discussing new developments in the MS and computational fields that will further enable in-cell structural studies.

Cryo-EM targets in CASP14.

Structures of seven CASP14 targets were determined using cryo-electron microscopy (cryo-EM) technique with resolution between 2.1 and 3.8 Å. We provide an evaluation of the submitted models versus the experimental data (cryo-EM density maps) and experimental reference structures built into the maps. The accuracy of models is measured in terms of coordinate-to-density and coordinate-to-coordinate fit. A-posteriori refinement of the most accurate models in their corresponding cryo-EM density resulted in structures that are close to the reference structure, including some regions with better fit to the density. Regions that were found to be less "refineable" correlate well with regions of high diversity between the CASP models and low goodness-of-fit to density in the reference structure.

Cyclic Ion Mobility-Collision Activation Experiments Elucidate Protein Behavior in the Gas Phase.

Ion mobility coupled to mass spectrometry (IM-MS) is widely used to study protein dynamics and structure in the gas phase. Increasing the energy with which the protein ions are introduced to the IM cell can induce them to unfold, providing information on the comparative energetics of unfolding between different proteoforms. Recently, a high-resolution cyclic IM-mass spectrometer (cIM-MS) was introduced, allowing multiple, consecutive tandem IM experiments (IMn) to be carried out. We describe a tandem IM technique for defining detailed protein unfolding pathways and the dynamics of disordered proteins. The method involves multiple rounds of IM separation and collision activation (CA): IM-CA-IM and CA-IM-CA-IM. Here, we explore its application to studies of a model protein, cytochrome C, and dimeric human islet amyloid polypeptide (hIAPP), a cytotoxic and amyloidogenic peptide involved in type II diabetes. In agreement with prior work using single stage IM-MS, several unfolding events are observed for cytochrome C. IMn-MS experiments also show evidence of interconversion between compact and extended structures. IMn-MS data for hIAPP shows interconversion prior to dissociation, suggesting that the certain conformations have low energy barriers between them and transition between compact and extended forms.

TEMPy2: a Python library with improved 3D electron microscopy density-fitting and validation workflows.

Structural determination of molecular complexes by cryo-EM requires large, often complex processing of the image data that are initially obtained. Here, TEMPy2, an update of the TEMPy package to process, optimize and assess cryo-EM maps and the structures fitted to them, is described. New optimization routines, comprehensive automated checks and workflows to perform these tasks are described.

Automated Modeling and Validation of Protein Complexes in Cryo-EM Maps.

The resolving power of cryo-EM experiments has dramatically improved in recent years. However, many cryo-EM maps may still not achieve a resolution that is sufficiently high to allow model building directly from the map. Instead, it is common practice to fit an initial atomic model to the map and refine this model. Depending on the resolution and whether the structure suffers from inherent flexibility or experimental limitations, different methods can be applied, to obtain high-quality, well-fitted atomic model of the macromolecular assembly represented by the map, and to assess its properties. In this review, we describe some of these methods, with the main focus on those that have been developed in our group over the last decade.

Cryo-electron microscopy targets in CASP13: Overview and evaluation of results.

Structures of seven CASP13 targets were determined using cryo-electron microscopy (cryo-EM) technique with resolution between 3.0 and 4.0 Å. We provide an overview of the experimentally derived structures and describe results of the numerical evaluation of the submitted models. The evaluation is carried out by comparing coordinates of models to those of reference structures (CASP-style evaluation), as well as checking goodness-of-fit of modeled structures to the cryo-EM density maps. The performance of contributing research groups in the CASP-style evaluation is measured in terms of backbone accuracy, all-atom local geometry and similarity of inter-subunit interfaces. The results on the cryo-EM targets are compared with those on the whole set of eighty CASP13 targets. A posteriori refinement of the best models in their corresponding cryo-EM density maps resulted in structures that are very close to the reference structure, including some regions with better fit to the density.