cryo-EM Density Research Articles

DAQ-refine: Protein structure model evaluation and refinement for cryo-EM maps.

As cryo-electron microscopy (cryo-EM) has become a more commonly used method in structural biology, many structure models of biomolecules have been determined from cryo-EM density maps. However, our recent research revealed that many protein structure models determined with cryo-EM and deposited in the Protein Data Bank (PDB) have potential errors. Such errors include cases where modeled local conformation is not accurate and situations where misassignment of amino acid types. Thus, establishing a computational protocol that evaluates the model accuracy and corrects low-quality regions in the model is crucial to ensuring the quality of structure models deposited to the public database, PDB. Here, we present a new protocol (DAQ-refine) for evaluating a protein model built from a cryo-EM map and for applying local structure refinement in case the model has potential errors. In the DAQ-refine protocol, model evaluation is performed with a deep learning-based model-local map assessment score, DAQ score, which we developed recently. Then, the subsequent local refinement is performed by a modified procedure of AlphaFold2. In the local refinement step, we used a trimmed template and trimmed multiple sequence alignment as input of AlphaFold2 to control which structure regions to refine while leaving other more confident regions in the model intact. DAQ-refine showed overall the highest accuracy in structure refinement when compared with four other methods on protein structures in PDB that were built from cryo-EM maps determined at a 3-4.5 Å resolution. DAQ score was also able to identify the most accurate model from about 20 alternative models generated for each target protein. The results demonstrated that our protocol, DAQ-refine, consistently improves low-quality regions of initial models.

Improving Protein-Ligand Interaction Modeling with cryo-EM Data, Templates, and Deep Learning in 2021 Ligand Model Challenge.

Elucidating protein-ligand interaction is crucial for studying the function of proteins and compounds in an organism and critical for drug discovery and design. The problem of protein-ligand interaction is traditionally tackled by molecular docking and simulation, which is based on physical forces and statistical potentials and cannot effectively leverage cryo-EM data and existing protein structural information in the protein-ligand modeling process. In this work, we developed a deep learning bioinformatics pipeline (DeepProLigand) to predict protein-ligand interactions from cryo-EM density maps of proteins and ligands. DeepProLigand first uses a deep learning method to predict the structure of proteins from cryo-EM maps, which is averaged with a reference (template) structure of the proteins to produce a combined structure to add ligands. The ligands are then identified and added into the structure to generate a protein-ligand complex structure, which is further refined. The method based on the deep learning prediction and template-based modeling was blindly tested in the 2021 EMDataResource Ligand Challenge and was ranked first in fitting ligands to cryo-EM density maps. These results demonstrate that the deep learning bioinformatics approach is a promising direction for modeling protein-ligand interactions on cryo-EM data using prior structural information.

Protein model refinement for cryo-EM maps using AlphaFold2 and the DAQ score.

As more protein structure models have been determined from cryogenic electron microscopy (cryo-EM) density maps, establishing how to evaluate the model accuracy and how to correct models in cases where they contain errors is becoming crucial to ensure the quality of the structural models deposited in the public database, the PDB. Here, a new protocol is presented for evaluating a protein model built from a cryo-EM map and applying local structure refinement in the case where the model has potential errors. Firstly, model evaluation is performed using a deep-learning-based model-local map assessment score, DAQ, that has recently been developed. The subsequent local refinement is performed by a modified AlphaFold2 procedure, in which a trimmed template model and a trimmed multiple sequence alignment are provided as input to control which structure regions to refine while leaving other more confident regions of the model intact. A benchmark study showed that this protocol, DAQ-refine, consistently improves low-quality regions of the initial models. Among 18 refined models generated for an initial structure, DAQ shows a high correlation with model quality and can identify the best accurate model for most of the tested cases. The improvements obtained by DAQ-refine were on average larger than other existing methods.

Protein Structure Predictions, Atomic Model Building, and Validation Using a Cryo-EM Density Map from Hepatitis B Virus Spherical Subviral Particle.

Hepatitis B virus (HBV) infection is a global public health concern. During chronic infection, the HBV small-surface antigen is expressed in large excess as non-infectious spherical subviral particles (SVPs), which possess strong immunogenicity. To date, attempts at understanding the structure of HBV spherical SVP have been restricted to 12-30 Å with contradictory conclusions regarding its architecture. We have used cryo-electron microscopy (cryo-EM) and 3D image reconstruction to solve the HBV spherical SVP to 6.3 Å. Here, we present an extended protocol on combining AlphaFold2 prediction with a moderate-resolution cryo-EM density map to build a reliable 3D model. This protocol utilizes multiple software packages that are routinely used in the cryo-EM community. The workflow includes 3D model prediction, model evaluation, rigid-body fitting, flexible fitting, real-space refinement, model validation, and model adjustment. Finally, the described protocol can also be applied to high-resolution cryo-EM datasets (2-4 Å).

Quantitative analysis of protein dynamics using a deep learning technique combined with experimental cryo-EM density data and MD simulations.

Protein functions associated with biological activity are precisely regulated by both tertiary structure and dynamic behavior. Thus, elucidating the high-resolution structures and quantitative information on in-solution dynamics is essential for understanding the molecular mechanisms. The main experimental approaches for determining tertiary structures include nuclear magnetic resonance (NMR), X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Among these procedures, recent remarkable advances in the hardware and analytical techniques of cryo-EM have increasingly determined novel atomic structures of macromolecules, especially those with large molecular weights and complex assemblies. In addition to these experimental approaches, deep learning techniques, such as AlphaFold 2, accurately predict structures from amino acid sequences, accelerating structural biology research. Meanwhile, the quantitative analyses of the protein dynamics are conducted using experimental approaches, such as NMR and hydrogen-deuterium mass spectrometry, and computational approaches, such as molecular dynamics (MD) simulations. Although these procedures can quantitatively explore dynamic behavior at high resolution, the fundamental difficulties, such as signal crowding and high computational cost, greatly hinder their application to large and complex biological macromolecules. In recent years, machine learning techniques, especially deep learning techniques, have been actively applied to structural data to identify features that are difficult for humans to recognize from big data. Here, we review our approach to accurately estimate dynamic properties associated with local fluctuations from three-dimensional cryo-EM density data using a deep learning technique combined with MD simulations.

Uncovering structural ensembles from single-particle cryo-EM data using cryoDRGN.

Single-particle cryogenic electron microscopy (cryo-EM) has emerged as a powerful technique to visualize the structural landscape sampled by a protein complex. However, algorithmic and computational bottlenecks in analyzing heterogeneous cryo-EM datasets have prevented the full realization of this potential. CryoDRGN is a machine learning system for heterogeneous cryo-EM reconstruction of proteins and protein complexes from single-particle cryo-EM data. Central to this approach is a deep generative model for heterogeneous cryo-EM density maps, which we empirically find is effective in modeling both discrete and continuous forms of structural variability. Once trained, cryoDRGN is capable of generating an arbitrary number of 3D density maps, and thus interpreting the resulting ensemble is a challenge. Here, we showcase interactive and automated processing approaches for analyzing cryoDRGN results. Specifically, we detail a step-by-step protocol for the analysis of an existing assembling 50S ribosome dataset, including preparation of inputs, network training and visualization of the resulting ensemble of density maps. Additionally, we describe and implement methods to comprehensively analyze and interpret the distribution of volumes with the assistance of an associated atomic model. This protocol is appropriate for structural biologists familiar with processing single-particle cryo-EM datasets and with moderate experience navigating Python and Jupyter notebooks. It requires 3-4 days to complete. CryoDRGN is open source software that is freely available.

Conformational motions and ligand-binding underlying gating and regulation in IP3R channel

Inositol-1,4,5-trisphosphate receptors (IP3Rs) are activated by IP3 and Ca2+ and their gating is regulated by various intracellular messengers that finely tune the channel activity. Here, using single particle cryo-EM analysis we determined 3D structures of the nanodisc-reconstituted IP3R1 channel in two ligand-bound states. These structures provide unprecedented details governing binding of IP3, Ca2+ and ATP, revealing conformational changes that couple ligand-binding to channel opening. Using a deep-learning approach and 3D variability analysis we extracted molecular motions of the key protein domains from cryo-EM density data. We find that IP3 binding relies upon intrinsic flexibility of the ARM2 domain in the tetrameric channel. Our results highlight a key role of dynamic side chains in regulating gating behavior of IP3R channels. This work represents a stepping-stone to developing mechanistic understanding of conformational pathways underlying ligand-binding, activation and regulation of the channel.

StarMap: a user-friendly workflow for Rosetta-driven molecular structure refinement.

Cryogenic electron microscopy (cryo-EM) data represent density maps of macromolecular systems at atomic or near-atomic resolution. However, building and refining 3D atomic models by using data from cryo-EM maps is not straightforward and requires significant hands-on experience and manual intervention. We recently developed StarMap, an easy-to-use interface between the popular structural display program ChimeraX and Rosetta, a powerful molecular modeling engine. StarMap offers a general approach for refining structural models of biological macromolecules into cryo-EM density maps by combining Monte Carlo sampling with local density-guided optimization, Rosetta-based all-atom refinement and real-space B-factor calculations in a straightforward workflow. StarMap includes options for structural symmetry, local refinements and independent model validation. The overall quality of the refinement and the structure resolution is then assessed via analytical outputs, such as magnification calibration (pixel size calibration) and Fourier shell correlations. Z-scores reported by StarMap provide an easily interpretable indicator of the goodness of fit for each residue and can be plotted to evaluate structural models and improve local residue refinements, as well as to identify flexible regions and potentially functional sites in large macromolecular complexes. The protocol requires general computer skills, without the need for coding expertise, because most parts of the workflow can be operated by clicking tabs within the ChimeraX graphical user interface. Time requirements for the model refinement depend on the size and quality of the input data; however, this step can typically be completed within 1 d. The analytical parts of the workflow are completed within minutes.

Auto-DRRAFTER: Automated RNA Modeling Based on Cryo-EM Density.

RNA three-dimensional structures provide rich and vital information for understanding their functions. Recent advances in cryogenic electron microscopy (cryo-EM) allow structure determination of RNAs and ribonucleoprotein (RNP) complexes. However, limited global and local resolutions of RNA cryo-EM mapspose great challenges in tracing RNA coordinates. The Rosetta-based "auto-DRRAFTER" method builds RNA models into moderate-resolution RNA cryo-EM density as part of the Ribosolve pipeline. Here, we describe a step-by-step protocol for auto-DRRAFTER using a glycine riboswitch from Fusobacterium nucleatum as an example. Successful implementation of this protocol allows automated RNA modeling into RNA cryo-EM density, accelerating our understanding of RNA structure-function relationships. Input and output files are being made available at https://github.com/auto-DRRAFTER/springer-chapter .

Building Protein Atomic Models from Cryo-EM Density Maps and Residue Co-Evolution.

Electron cryo-microscopy (cryo-EM) has emerged as a powerful method by which to obtain three-dimensional (3D) structures of macromolecular complexes at atomic or near-atomic resolution. However, de novo building of atomic models from near-atomic resolution (3–5 Å) cryo-EM density maps is a challenging task, in particular because poorly resolved side-chain densities hamper sequence assignment by automatic procedures at a lower resolution. Furthermore, segmentation of EM density maps into individual subunits remains a difficult problem when the structure of the subunits is not known, or when significant conformational rearrangement occurs between the isolated and associated form of the subunits. To tackle these issues, we have developed a graph-based method to thread most of the C- trace of the protein backbone into the EM density map. The EM density is described as a weighted graph such that the resulting minimum spanning tree encompasses the high-density regions of the map. A pruning algorithm cleans the tree and finds the most probable positions of the C- atoms, by using side-chain density when available, as a collection of C- trace fragments. By complementing experimental EM maps with contact predictions from sequence co-evolutionary information, we demonstrate that this approach can correctly segment EM maps into individual subunits and assign amino acid sequences to backbone traces to generate atomic models.

Residue-wise local quality estimation for protein models from cryo-EM maps.

An increasing number of protein structures are being determined by cryogenic electron microscopy (cryo-EM). Although the resolution of determined cryo-EM density maps is improving in general, there are still many cases where amino acids of a protein are assigned with different levels of confidence. Here we developed a method that identifies potential misassignment of residues in the map, including residue shifts along an otherwise correct main-chain trace. The score, named DAQ, computes the likelihood that the local density corresponds to different amino acids, atoms, and secondary structures, estimated via deep learning, and assesses the consistency of the amino acid assignment in the protein structure model with that likelihood. When DAQ was applied to different versions of model structures in the Protein Data Bank that were derived from the same density maps, a clear improvement in the DAQ score was observed in the newer versions of the models. DAQ also found potential misassignment errors in a substantial number of deposited protein structure models built into cryo-EM maps.

Beyond the Backbone: The Next Generation of Pathwalking Utilities for Model Building in CryoEM Density Maps.

Single-particle electron cryomicroscopy (cryoEM) has become an indispensable tool for studying structure and function in macromolecular assemblies. As an integral part of the cryoEM structure determination process, computational tools have been developed to build atomic models directly from a density map without structural templates. Nearly a decade ago, we created Pathwalking, a tool for de novo modeling of protein structure in near-atomic resolution cryoEM density maps. Here, we present the latest developments in Pathwalking, including the addition of probabilistic models, as well as a companion tool for modeling waters and ligands. This software was evaluated on the 2021 CryoEM Ligand Challenge density maps, in addition to identifying ligands in three IP3R1 density maps at ~3 Å to 4.1 Å resolution. The results clearly demonstrate that the Pathwalking de novo modeling pipeline can construct accurate protein structures and reliably localize and identify ligand density directly from a near-atomic resolution map.

AlphaFold2 and CryoEM: Revisiting CryoEM modeling in near-atomic resolution density maps

SummaryWith the advent of new artificial intelligence and machine learning algorithms, predictive modeling can, in some cases, produce structures on par with experimental methods. The combination of predictive modeling and experimental structure determination by electron cryomicroscopy (cryoEM) offers a tantalizing approach for producing robust atomic models of macromolecular assemblies. Here, we apply AlphaFold2 to a set of community standard data sets and compare the results with the corresponding reference maps and models. Moreover, we present three unique case studies from previously determined cryoEM density maps of viruses. Our results show that AlphaFold2 can not only produce reasonably accurate models for analysis and additional hypotheses testing, but can also potentially yield incorrect structures if not properly validated with experimental data. Whereas we outline numerous shortcomings and potential pitfalls of predictive modeling, the obvious synergy between predictive modeling and cryoEM will undoubtedly result in new computational modeling tools.

Assessing the Mobility of Severe Acute Respiratory Syndrome Coronavirus-2 Spike Protein Glycans by Structural and Computational Methods.

Two years after its emergence, the coronavirus disease-2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) remains difficult to control despite the availability of several vaccines. The extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates host cell entry by binding to the angiotensin converting enzyme 2 (ACE2) through its receptor binding domain (RBD), is the major target of neutralizing antibodies. Like to many other viral fusion proteins, the SARS-CoV-2 spike protein utilizes a glycan shield to thwart the host immune response. To grasp the influence of chemical signatures on carbohydrate mobility and reconcile the cryo-EM density of specific glycans we combined our cryo-EM map of the S ectodomain to 4.1 Å resolution, reconstructed from a limited number of particles, and all-atom molecular dynamics simulations. Chemical modifications modeled on representative glycans (defucosylation, sialylation and addition of terminal LacNAc units) show no significant influence on either protein shielding or glycan flexibility. By estimating at selected sites the local correlation between the full density map and atomic model-based maps derived from molecular dynamics simulations, we provide insight into the geometries of the α-Man-(1→3)-[α-Man-(1→6)-]-β-Man-(1→4)-β-GlcNAc(1→4)-β-GlcNAc core common to all N-glycosylation sites.

Progressive assembly of multi-domain protein structures from cryo-EM density maps.

Progress in cryo-electron microscopy has provided the potential for large-size protein structure determination. However, the success rate for solving multi-domain proteins remains low because of the difficulty in modelling inter-domain orientations. Here we developed domain enhanced modeling using cryo-electron microscopy (DEMO-EM), an automatic method to assemble multi-domain structures from cryo-electron microscopy maps through a progressive structural refinement procedure combining rigid-body domain fitting and flexible assembly simulations with deep-neural-network inter-domain distance profiles. The method was tested on a large-scale benchmark set of proteins containing up to 12 continuous and discontinuous domains with medium- to low-resolution density maps, where DEMO-EM produced models with correct inter-domain orientations (template modeling score (TM-score) >0.5) for 97% of cases and outperformed state-of-the-art methods. DEMO-EM was applied to the severe acute respiratory syndrome coronavirus 2 genome and generated models with average TM-score and root-mean-square deviation of 0.97 and 1.3 Å, respectively, with respect to the deposited structures. These results demonstrate an efficient pipeline that enables automated and reliable large-scale multi-domain protein structure modelling from cryo-electron microscopy maps.

LPTD: a novel linear programming-based topology determination method for cryo-EM maps.

Topology determination is one of the most important intermediate steps toward building the atomic structure of proteins from their medium-resolution cryo-electron microscopy (cryo-EM) map. The main goal in the topology determination is to identify correct matches (i.e. assignment and direction) between secondary structure elements (SSEs) (α-helices and β-sheets) detected in a protein sequence and cryo-EM density map. Despite many recent advances in molecular biology technologies, the problem remains a challenging issue. To overcome the problem, this article proposes a linear programming-based topology determination (LPTD) method to solve the secondary structure topology problem in three-dimensional geometrical space. Through modeling of the protein's sequence with the aid of extracting highly reliable features and a distance-based scoring function, the secondary structure matching problem is transformed into a complete weighted bipartite graph matching problem. Subsequently, an algorithm based on linear programming is developed as a decision-making strategy to extract the true topology (native topology) between all possible topologies. The proposed automatic framework is verified using 12 experimental and 15 simulated α-β proteins. Results demonstrate that LPTD is highly efficient and extremely fast in such a way that for 77% of cases in the dataset, the native topology has been detected in the first rank topology in <2 s. Besides, this method is able to successfully handle large complex proteins with as many as 65 SSEs. Such a large number of SSEs have never been solved with current tools/methods. The LPTD package (source code and data) is publicly available at https://github.com/B-Behkamal/LPTD. Moreover, two test samples as well as the instruction of utilizing the graphical user interface have been provided in the shared readme file. Supplementary data are available at Bioinformatics online.

Structural Comparative Modeling of Multi-Domain F508del CFTR.

Cystic fibrosis (CF) is a rare genetic disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR), an epithelial anion channel expressed in several vital organs. Absence of functional CFTR results in imbalanced osmotic equilibrium and subsequent mucus build up in the lungs-which increases the risk of infection and eventually causes death. CFTR is an ATP-binding cassette (ABC) transporter family protein composed of two transmembrane domains (TMDs), two nucleotide binding domains (NBDs), and an unstructured regulatory domain. The most prevalent patient mutation is the deletion of F508 (F508del), making F508del CFTR the primary target for current FDA approved CF therapies. However, no experimental multi-domain F508del CFTR structure has been determined and few studies have modeled F508del using multi-domain WT CFTR structures. Here, we used cryo-EM density data and Rosetta comparative modeling (RosettaCM) to compare a F508del model with published experimental data on CFTR NBD1 thermodynamics. We then apply this modeling method to generate multi-domain WT and F508del CFTR structural models. These models demonstrate the destabilizing effects of F508del on NBD1 and the NBD1/TMD interface in both the inactive and active conformation of CFTR. Furthermore, we modeled F508del/R1070W and F508del bound to the CFTR corrector VX-809. Our models reveal the stabilizing effects of VX-809 on multi-domain models of F508del CFTR and pave the way for rational design of additional drugs that target F508del CFTR for treatment of CF.

Source code and data for the paper 'Progressive assembly of multi-domain protein structures from cryo-EM density maps'

Source code and data for the paper 'Progressive assembly of multi-domain protein structures from cryo-EM density maps'

Progressive assembly of multi-domain protein structures from cryo-EM density maps

Progressive assembly of multi-domain protein structures from cryo-EM density maps

Progressive assembly of multi-domain protein structures from cryo-EM density maps

Progressive assembly of multi-domain protein structures from cryo-EM density maps