Modeling tricks and fitting techniques for multiresolution structures

Observation of Bacteriophage Ultrastructure by Cryo-Electron Microscopy

Transmission electron microscopy (TEM) is an ideal method to observe and determine the structure of bacteriophages. From early studies by negative staining to the present atomic structure models derived from cryo-TEM, bacteriophage detection, classification, and structure determination have been mostly done by electron microscopy. Although embedding in metal salts has been a routine method for virus observation for many years, the preservation of bacteriophages in a thin layer of fast frozen buffer has proven to be the most convenient preparation method for obtaining images using cryo-electron microscopy (cryo-EM). In this technique, frozen samples are observed at liquid nitrogen temperature, and the images are acquired using different recording media. The incorporation of direct electron detectors has been a fundamental step in achieving atomic resolution images of a number of viruses. These projection images can be numerically combined using different approaches to render a three-dimensional model of the virus. For those viral components exhibiting any symmetry, averaging can nowadays achieve atomic structures in most cases. Image processing methods have also evolved to improve the resolution in asymmetric viral components or regions showing different types of symmetries (symmetry mismatch).

Compounds derived from Humulus lupulus inhibit SARS-CoV-2 papain-like protease and virus replication

BACKGROUND

Selected natural compounds exhibit very good antiviral properties. Especially, the medicinal plant Humulus lupulus (hop) contains several secondary plant metabolites some of which have previously shown antiviral activities. Among them, the prenylated chalcone xanthohumol (XN) demonstrated to be a potent inhibitor of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro).

HYPOTHESIS/PURPOSE

Following the finding that xanthohumol (XN) is a potent inhibitor of SARS-CoV-2 Mpro, the effect of XN and its major derivatives isoxanthohumol (IXN), 6-prenylnaringenin (6-PN), and 8-prenylnaringenin (8-PN) from hops on SARS-CoV-2 papain-like protease (PLpro) were investigated.

STUDY DESIGN

The modulatory effect of the hop compounds on PLpro were studied first in silico and then in vitro. In addition, the actual effect of hop compounds on the replication of SARS-CoV-2 in host cells was investigated.

METHODS

In silico docking analysis was used to predict the binding affinity of hop compounds to the active site of PLpro. A recombinant PLpro was cloned, purified, characterized, and analyzed by small-angle X-ray scattering (SAXS), deISGylation assays, and kinetic analyses. Antiviral activity of hop compounds was assessed using the fluorescently labeled wildtype SARS-CoV-2 (icSARS-CoV-2-mNG) in Caco-2 host cells.

RESULTS

Our in silico docking suggests that the purified hop compounds bind to the active site of SARS-CoV-2 PLpro blocking the access of its natural substrates. The hop-derived compounds inhibit SARS-CoV-2 PLpro with half maximal inhibitory concentration (IC50) values in the range of 59-162 µM. Furthermore, we demonstrate that XN and 6-PN, in particular, impede viral replication with IC50 values of 3.3 µM and 7.3 µM, respectively.

CONCLUSION

In addition to the already known inhibition of Mpro by XN, our results show, for the first time, that hop-derived compounds target also SARS-CoV-2 PLpro which is a promising therapeutic target as it contributes to both viral replication and modulation of the immune system. These findings support the possibility to develop new hop-derived antiviral drugs targeting human coronaviruses.

Structural dynamics of SARS-CoV-2 nucleocapsid protein induced by RNA binding

The nucleocapsid (N) protein of the SARS-CoV-2 virus, the causal agent of COVID-19, is a multifunction phosphoprotein that plays critical roles in the virus life cycle, including transcription and packaging of the viral RNA. To play such diverse roles, the N protein has two globular RNA-binding modules, the N- (NTD) and C-terminal (CTD) domains, which are connected by an intrinsically disordered region. Despite the wealth of structural data available for the isolated NTD and CTD, how these domains are arranged in the full-length protein and how the oligomerization of N influences its RNA-binding activity remains largely unclear. Herein, using experimental data from electron microscopy and biochemical/biophysical techniques combined with molecular modeling and molecular dynamics simulations, we show that, in the absence of RNA, the N protein formed structurally dynamic dimers, with the NTD and CTD arranged in extended conformations. However, in the presence of RNA, the N protein assumed a more compact conformation where the NTD and CTD are packed together. We also provided an octameric model for the full-length N bound to RNA that is consistent with electron microscopy images of the N protein in the presence of RNA. Together, our results shed new light on the dynamics and higher-order oligomeric structure of this versatile protein.

The nucleocapsid (N) protein of the SARS-CoV-2 virus plays an essential role in virus particle assembly as it specifically binds to and wraps the virus genomic RNA into a well-organized structure known as the ribonucleoprotein. Understanding how the N protein wraps around the virus RNA is critical for the development of strategies to inhibit virus assembly within host cells. One of the limitations regarding the molecular structure of the ribonucleoprotein, however, is that the N protein has several unstructured and mobile regions that preclude the resolution of its full atomic structure. Moreover, the N protein can form higher-order oligomers, both in the presence and absence of RNA. Here we employed computational methods, supported by experimental data, to simulate the N protein structural dynamics in the absence and presence of RNA. Our data suggest that the N protein forms structurally dynamic dimers in the absence of RNA, with its structured N- and C-terminal domains oriented in extended conformations. In the presence of RNA, however, the N protein assumes a more compact conformation. Our model for the oligomeric structure of the N protein bound to RNA helps to understand how N protein dimers interact to each other to form the ribonucleoprotein.

Three-Dimensional Graph Matching to Identify Secondary Structure Correspondence of Medium-Resolution Cryo-EM Density Maps

Cryo-electron microscopy (cryo-EM) is a structural technique that has played a significant role in protein structure determination in recent years. Compared to the traditional methods of X-ray crystallography and NMR spectroscopy, cryo-EM is capable of producing images of much larger protein complexes. However, cryo-EM reconstructions are limited to medium-resolution (~4-10 Å) for some cases. At this resolution range, a cryo-EM density map can hardly be used to directly determine the structure of proteins at atomic level resolutions, or even at their amino acid residue backbones. At such a resolution, only the position and orientation of secondary structure elements (SSEs) such as α-helices and β-sheets are observable. Consequently, finding the mapping of the secondary structures of the modeled structure (SSEs-A) to the cryo-EM map (SSEs-C) is one of the primary concerns in cryo-EM modeling. To address this issue, this study proposes a novel automatic computational method to identify SSEs correspondence in three-dimensional (3D) space. Initially, through a modeling of the target sequence with the aid of extracting highly reliable features from a generated 3D model and map, the SSEs matching problem is formulated as a 3D vector matching problem. Afterward, the 3D vector matching problem is transformed into a 3D graph matching problem. Finally, a similarity-based voting algorithm combined with the principle of least conflict (PLC) concept is developed to obtain the SSEs correspondence. To evaluate the accuracy of the method, a testing set of 25 experimental and simulated maps with a maximum of 65 SSEs is selected. Comparative studies are also conducted to demonstrate the superiority of the proposed method over some state-of-the-art techniques. The results demonstrate that the method is efficient, robust, and works well in the presence of errors in the predicted secondary structures of the cryo-EM images.

Simulation-Based Methods for Model Building and Refinement in Cryoelectron Microscopy

Advances in cryoelectron microscopy (cryo-EM) have revolutionized the structural investigation of large macromolecular assemblies. In this review, we first provide a broad overview of modeling methods used for flexible fitting of molecular models into cryo-EM density maps. We give special attention to approaches rooted in molecular simulations-atomistic molecular dynamics and Monte Carlo. Concise descriptions of the methods are given along with discussion of their advantages, limitations, and most popular alternatives. We also describe recent extensions of the widely used molecular dynamics flexible fitting (MDFF) method and discuss how different model-building techniques could be incorporated into new hybrid modeling schemes and simulation workflows. Finally, we provide two illustrative examples of model-building and refinement strategies employing MDFF, cascade MDFF, and RosettaCM. These examples come from recent cryo-EM studies that elucidated transcription preinitiation complexes and shed light on the functional roles of these assemblies in gene expression and gene regulation.

Structure and transformation of bacteriophage A511 baseplate and tail upon infection of Listeria cells

Contractile injection systems (bacteriophage tails, type VI secretions system, R-type pyocins, etc.) utilize a rigid tube/contractile sheath assembly for breaching the envelope of bacterial and eukaryotic cells. Among contractile injection systems, bacteriophages that infect Gram-positive bacteria represent the least understood members. Here, we describe the structure of Listeria bacteriophage A511 tail in its pre- and post-host attachment states (extended and contracted, respectively) using cryo-electron microscopy, cryo-electron tomography, and X-ray crystallography. We show that the structure of the tube-baseplate complex of A511 is similar to that of phage T4, but the A511 baseplate is decorated with different receptor-binding proteins, which undergo a large structural transformation upon host attachment and switch the symmetry of the baseplate-tail fiber assembly from threefold to sixfold. For the first time under native conditions, we show that contraction of the phage tail sheath assembly starts at the baseplate and propagates through the sheath in a domino-like motion.

Acceleration of cryo-EM Flexible Fitting for Large Biomolecular Systems by Efficient Space Partitioning

Flexible fitting is a powerful technique to build the 3D structures of biomolecules from cryoelectron microscopy (cryo-EM) density maps. One popular method is a cross-correlation coefficient-based approach, where the molecular dynamics (MD) simulation is carried out with the biasing potential that includes the cross-correlation coefficient between the experimental and simulated density maps. Here, we propose efficient parallelization schemes for the calculation of the cross-correlation coefficient to accelerate flexible fitting. Our schemes are tested for small, medium, and large biomolecules using CPU and hybrid CPU + GPU architectures. The scheme for the atomic decomposition MD is suitable for small proteins such as Ca²⁺-ATPase with the all-atom Go model, while that for the domain decomposition MD is better for larger systems such as ribosome with the all-atom Go or the all-atom explicit solvent models. Our methods allow flexible fitting for various biomolecules with reasonable computational cost. This approach also connects high-resolution structure refinements with investigation of protein structure-function relationship.

Structures of ubiquitin-like (Ubl) and Hsp90-like domains of sacsin provide insight into pathological mutations

Autosomal recessive spastic ataxia of Charlevoix-Saguenay (ARSACS) is a neurodegenerative disease that is caused by mutations in the SACS gene. The product of this gene is a very large 520-kDa cytoplasmic protein, sacsin, with a ubiquitin-like (Ubl) domain at the N terminus followed by three large sacsin internal repeat (SIRPT) supradomains and C-terminal J and HEPN domains. The SIRPTs are predicted to contain Hsp90-like domains, suggesting a potential chaperone activity. In this work, we report the structures of the Hsp90-like Sr1 domain of SIRPT1 and the N-terminal Ubl domain determined at 1.55- and 2.1-Å resolutions, respectively. The Ubl domain crystallized as a swapped dimer that could be relevant in the context of full-length protein. The Sr1 domain displays the Bergerat protein fold with a characteristic nucleotide-binding pocket, although it binds nucleotides with very low affinity. The Sr1 structure reveals that ARSACS-causing missense mutations (R272H, R272C, and T201K) disrupt protein folding, most likely leading to sacsin degradation. This work lends structural support to the view of sacsin as a molecular chaperone and provides a framework for future studies of this protein.

CryoEM-based hybrid modeling approaches for structure determination

Recent advances in cryo-electron microscopy (cryoEM) have dramatically improved the resolutions at which vitrified biological specimens can be studied, revealing new structural and mechanistic insights over a broad range of spatial scales. Bolstered by these advances, much effort has been directed toward the development of hybrid modeling methodologies for the construction and refinement of high-fidelity atomistic models from cryoEM data. In this brief review, we will survey the key elements of cryoEM-based hybrid modeling, providing an overview of available computational tools and strategies as well as several recent applications.

Analytical Approaches to Improve Accuracy in Solving the Protein Topology Problem

To take advantage of recent advances in genomics and proteomics it is critical that the three-dimensional physical structure of biological macromolecules be determined. Cryo-Electron Microscopy (cryo-EM) is a promising and improving method for obtaining this data, however resolution is often not sufficient to directly determine the atomic scale structure. Despite this, information for secondary structure locations is detectable. De novo modeling is a computational approach to modeling these macromolecular structures based on cryo-EM derived data. During de novo modeling a mapping between detected secondary structures and the underlying amino acid sequence must be identified. DP-TOSS (Dynamic Programming for determining the Topology Of Secondary Structures) is one tool that attempts to automate the creation of this mapping. By treating the correspondence between the detected structures and the structures predicted from sequence data as a constraint graph problem DP-TOSS achieved good accuracy in its original iteration. In this paper, we propose modifications to the scoring methodology of DP-TOSS to improve its accuracy. Three scoring schemes were applied to DP-TOSS and tested: (i) a skeleton-based scoring function; (ii) a geometry-based analytical function; and (iii) a multi-well potential energy-based function. A test of 25 proteins shows that a combination of these schemes can improve the performance of DP-TOSS to solve the topology determination problem for macromolecule proteins.

Observation of Bacteriophage Ultrastructure by Cryo-electron Microscopy

Transmission Electron Microscopy (TEM) is an ideal method to observe and determine the structure of bacteriophages. From early studies by negative staining to the present atomic structure models derived from cryo-TEM, bacteriophage detection, classification, and structure determination has been mostly done by electron microscopy. Although embedding in metal salts has been a routine method for virus observation for many years, preservation of bacteriophages in a thin layer of fast frozen buffer has proven to be a most convenient preparation method for obtaining images using cryo-electron microscopy (cryo-EM). In this technique, frozen samples are observed at liquid nitrogen temperature and the images are acquired using different recording media. The incorporation of direct electron detectors has been a fundamental step to achieve atomic resolution images of a number of viruses. These projection images can be numerically combined using different approaches to render a three-dimensional model of the virus. For those viral components exhibiting any symmetry, averaging procedures help to render near-atomic resolution structures.

PEM-fitter: A Coarse-Grained Method to Validate Protein Candidate Models

The volumetric images produced by Cryo-Electron Microscopy (cryo-EM) technique are used to model macromolecular assemblies and machines. De novo protein modeling uses these images to computationally model the structure of the molecules. Many candidate conformations are usually generated during the intermediate step. Conventionally, each of these candidates is evaluated by time-consuming approaches such as potential energy. We introduce an initial version of a geometrical screening method that uses the skeleton of the cryo-EM images to evaluate candidate structures. The aim of this method is to reduce the number of native-like candidate conformations and, therefore, reduce the time required for structural evaluation by energy calculations. A test of two datasets was performed. The first dataset contains 10 proteins and shows that our method can successfully detect the correct native structure for the given skeleton among a set of different protein structures. The second dataset contains 12 proteins and shows that our method can filter slightly modified decoy conformations of the same protein. The efficiency of the method is also reported.

Flexible Fitting of Atomic Models into Cryo-EM Density Maps Guided by Helix Correspondences

Although electron cryo-microscopy (cryo-EM) has recently achieved resolutions of better than 3 Å, at which point molecular modeling can be done directly from the density map, analysis and annotation of a cryo-EM density map still primarily rely on fitting atomic or homology models to the density map. In this article, we present, to our knowledge, a new method for flexible fitting of known or modeled protein structures into cryo-EM density maps. Unlike existing methods that are guided by local density gradients, our method is guided by correspondences between the α-helices in the density map and model, and does not require an initial rigid-body fitting step. Compared with current methods on both simulated and experimental density maps, our method not only achieves greater accuracy for proteins with large deformations but also runs as fast or faster than many of the other flexible fitting routines.

A novel fast Fourier transform accelerated off-grid exhaustive search method for cryo-electron microscopy fitting

This paper presents a novel fast Fourier transform (FFT)-based exhaustive search method extended to off-grid translational and rotational degrees of freedom. The method combines the advantages of the FFT-based exhaustive search, which samples all the conformations of a system under study on a grid, with a local optimization technique that guarantees to find the nearest optimal off-grid conformation. The method is demonstrated on a fitting problem and can be readily applied to a docking problem. The algorithm first samples a scoring function on a six-dimensional grid of sizeN6using the FFT. This operation has an asymptotic complexity ofO(N6logN). Then, the method performs the off-grid search using a local quadratic approximation of the cost function and the trust-region optimization algorithm. The computation of the quadratic approximation is also accelerated by FFT at the same additional asymptotic cost ofO(N6logN). The method is demonstrated by fitting atomic protein models into several simulated and experimental maps from cryo-electron microscopy. The method is available at https://team.inria.fr/nano-d/software/offgridfit.

Variability of Protein Structure Models from Electron Microscopy

An increasing number of biomolecular structures are solved by electron microscopy (EM). However, the quality of structure models determined from EM maps vary substantially. To understand to what extent structure models are supported by information embedded in EM maps, we used two computational structure refinement methods to examine how much structures can be refined using a dataset of 49 maps with accompanying structure models. The extent of structure modification as well as the disagreement between refinement models produced by the two computational methods scaled inversely with the global and the local map resolutions. A general quantitative estimation of deviations of structures for particular map resolutions are provided. Our results indicate that the observed discrepancy between the deposited map and the refined models is due to the lack of structural information present in EM maps and thus these annotations must be used with caution for further applications.

SAXS Structural Studies of Dps from Deinococcus radiodurans Highlights the Conformation of the Mobile N-Terminal Extensions

The radiation-resistant bacterium Deinococcus radiodurans contains two DNA-binding proteins from starved cells (Dps): Dps1 (DR2263) and Dps2 (DRB0092). These are suggested to play a role in DNA interaction and manganese and iron storage. The proteins assemble as a conserved dodecameric structure with structurally uncharacterised N-terminal extensions. In the case of DrDps1, these extensions have been proposed to be involved in DNA interactions, while in DrDps2, their function has yet to be established. The reported data reveal the relative position of the N-terminal extensions to the dodecameric sphere in solution for both Dps. The low-resolution small angle X-ray scattering (SAXS) results show that the N-terminal extensions protrude from the spherical shell of both proteins. The SAXS envelope of a truncated form of DrDps1 without the N-terminal extensions appears as a dodecameric sphere, contrasting strongly with the protrusions observed in the full-length models. The effect of iron incorporation into DrDps2 was investigated by static and stopped-flow SAXS measurements, revealing dynamic structural changes upon iron binding and core formation, as reflected by a quick alteration of its radius of gyration. The truncated and full-length versions of DrDps were also compared on the basis of their interaction with DNA to analyse functional roles of the N-terminal extensions. DrDps1 N-terminal protrusions appear to be directly involved with DNA, whilst those from DrDps2 are indirectly associated with DNA binding. Furthermore, detection of DrDps2 in the D. radiodurans membrane fraction suggests that the N-terminus of the protein interacts with the membrane.

Honey Bee Deformed Wing Virus Structures Reveal that Conformational Changes Accompany Genome Release

The picornavirus-like deformed wing virus (DWV) has been directly linked to colony collapse; however, little is known about the mechanisms of host attachment or entry for DWV or its molecular and structural details. Here we report the three-dimensional (3-D) structures of DWV capsids isolated from infected honey bees, including the immature procapsid, the genome-filled virion, the putative entry intermediate (A-particle), and the empty capsid that remains after genome release. The capsids are decorated by large spikes around the 5-fold vertices. The 5-fold spikes had an open flower-like conformation for the procapsid and genome-filled capsids, whereas the putative A-particle and empty capsids that had released the genome had a closed tube-like spike conformation. Between the two conformations, the spikes undergo a significant hinge-like movement that we predicted using a Robetta model of the structure comprising the spike. We conclude that the spike structures likely serve a function during host entry, changing conformation to release the genome, and that the genome may escape from a 5-fold vertex to initiate infection. Finally, the structures illustrate that, similarly to picornaviruses, DWV forms alternate particle conformations implicated in assembly, host attachment, and RNA release.

IMPORTANCE

Honey bees are critical for global agriculture, but dramatic losses of entire hives have been reported in numerous countries since 2006. Deformed wing virus (DWV) and infestation with the ectoparasitic mite Varroa destructor have been linked to colony collapse disorder. DWV was purified from infected adult worker bees to pursue biochemical and structural studies that allowed the first glimpse into the conformational changes that may be required during transmission and genome release for DWV.

Near-Atomic Resolution Structure of a Highly Neutralizing Fab Bound to Canine Parvovirus

Canine parvovirus (CPV) is a highly contagious pathogen that causes severe disease in dogs and wildlife. Previously, a panel of neutralizing monoclonal antibodies (MAb) raised against CPV was characterized. An antibody fragment (Fab) of MAb E was found to neutralize the virus at low molar ratios. Using recent advances in cryo-electron microscopy (cryo-EM), we determined the structure of CPV in complex with Fab E to 4.1 Å resolution, which allowed de novo building of the Fab structure. The footprint identified was significantly different from the footprint obtained previously from models fitted into lower-resolution maps. Using single-chain variable fragments, we tested antibody residues that control capsid binding. The near-atomic structure also revealed that Fab binding had caused capsid destabilization in regions containing key residues conferring receptor binding and tropism, which suggests a mechanism for efficient virus neutralization by antibody. Furthermore, a general technical approach to solving the structures of small molecules is demonstrated, as binding the Fab to the capsid allowed us to determine the 50-kDa Fab structure by cryo-EM.

IMPORTANCE

Using cryo-electron microscopy and new direct electron detector technology, we have solved the 4 Å resolution structure of a Fab molecule bound to a picornavirus capsid. The Fab induced conformational changes in regions of the virus capsid that control receptor binding. The antibody footprint is markedly different from the previous one identified by using a 12 Å structure. This work emphasizes the need for a high-resolution structure to guide mutational analysis and cautions against relying on older low-resolution structures even though they were interpreted with the best methodology available at the time.

Tools for Model Building and Optimization into Near-Atomic Resolution Electron Cryo-Microscopy Density Maps

Electron cryo-microscopy (cryoEM) has advanced dramatically to become a viable tool for high-resolution structural biology research. The ultimate outcome of a cryoEM study is an atomic model of a macromolecule or its complex with interacting partners. This chapter describes a variety of algorithms and software to build a de novo model based on the cryoEM 3D density map, to optimize the model with the best stereochemistry restraints and finally to validate the model with proper protocols. The full process of atomic structure determination from a cryoEM map is described. The tools outlined in this chapter should prove extremely valuable in revealing atomic interactions guided by cryoEM data.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling

Molecular Dynamics Flexible Fitting (MDFF) is an established technique for fitting all-atom structures of molecules into corresponding cryo-electron microscopy (cryo-EM) densities. The practical application of MDFF is simple but requires a user to be aware of and take measures against a variety of possible challenges presented by each individual case. Some of these challenges arise from the complexity of a molecular structure or the limited quality of available structural models and densities to be interpreted, while others stem from the intricacies of MDFF itself. The current article serves as an overview of the strategies that have been developed since MDFF's inception to overcome common challenges and successfully perform MDFF simulations.

Conformational Disorder of the Most Immature Cu, Zn-Superoxide Dismutase Leading to Amyotrophic Lateral Sclerosis

Misfolding of Cu,Zn-superoxide dismutase (SOD1) is a pathological change in the familial form of amyotrophic lateral sclerosis caused by mutations in the SOD1 gene. SOD1 is an enzyme that matures through the binding of copper and zinc ions and the formation of an intramolecular disulfide bond. Pathogenic mutations are proposed to retard the post-translational maturation, decrease the structural stability, and hence trigger the misfolding of SOD1 proteins. Despite this, a misfolded and potentially pathogenic conformation of immature SOD1 remains obscure. Here, we show significant and distinct conformational changes of apoSOD1 that occur only upon reduction of the intramolecular disulfide bond in solution. In particular, loop regions in SOD1 lose their restraint and become significantly disordered upon dissociation of metal ions and reduction of the disulfide bond. Such drastic changes in the solution structure of SOD1 may trigger misfolding and fibrillar aggregation observed as pathological changes in the familial form of amyotrophic lateral sclerosis.

Conserved Intramolecular Interactions Maintain Myosin Interacting-Heads Motifs Explaining Tarantula Muscle Super-Relaxed State Structural Basis

Tarantula striated muscle is an outstanding system for understanding the molecular organization of myosin filaments. Three-dimensional reconstruction based on cryo-electron microscopy images and single-particle image processing revealed that, in a relaxed state, myosin molecules undergo intramolecular head-head interactions, explaining why head activity switches off. The filament model obtained by rigidly docking a chicken smooth muscle myosin structure to the reconstruction was improved by flexibly fitting an atomic model built by mixing structures from different species to a tilt-corrected 2-nm three-dimensional map of frozen-hydrated tarantula thick filament. We used heavy and light chain sequences from tarantula myosin to build a single-species homology model of two heavy meromyosin interacting-heads motifs (IHMs). The flexibly fitted model includes previously missing loops and shows five intramolecular and five intermolecular interactions that keep the IHM in a compact off structure, forming four helical tracks of IHMs around the backbone. The residues involved in these interactions are oppositely charged, and their sequence conservation suggests that IHM is present across animal species. The new model, PDB 3JBH, explains the structural origin of the ATP turnover rates detected in relaxed tarantula muscle by ascribing the very slow rate to docked unphosphorylated heads, the slow rate to phosphorylated docked heads, and the fast rate to phosphorylated undocked heads. The conservation of intramolecular interactions across animal species and the presence of IHM in bilaterians suggest that a super-relaxed state should be maintained, as it plays a role in saving ATP in skeletal, cardiac, and smooth muscles.

Improving the Accuracy of Fitted Atomic Models in Cryo-EM Density Maps of Protein Assemblies Using Evolutionary Information from Aligned Homologous Proteins

Cryo-Electron Microscopy (cryo-EM) has become an important technique to obtain structural insights into large macromolecular assemblies. However the resolution of the density maps do not allow for its interpretation at atomic level. Hence they are combined with high resolution structures along with information from other experimental or bioinformatics techniques to obtain pseudo-atomic models. Here, we describe the use of evolutionary conservation of residues as obtained from protein structures and alignments of homologous proteins to detect errors in the fitting of atomic structures as well as improve accuracy of the protein-protein interfacial regions in the cryo-EM density maps.

Molecular modeling of Plasmodium falciparum peptide deformylase and structure-based pharmacophore screening for inhibitors

The role of metal coordination geometry and actinonin (inhibitor) binding was examined to develop pharmacophore-based inhibitor design strategy forPlasmodium falciparumpeptide deformylase.

PF2fit: Polar Fast Fourier Matched Alignment of Atomistic Structures with 3D Electron Microscopy Maps

There continue to be increasing occurrences of both atomistic structure models in the PDB (possibly reconstructed from X-ray diffraction or NMR data), and 3D reconstructed cryo-electron microscopy (3D EM) maps (albeit at coarser resolution) of the same or homologous molecule or molecular assembly, deposited in the EMDB. To obtain the best possible structural model of the molecule at the best achievable resolution, and without any missing gaps, one typically aligns (match and fits) the atomistic structure model with the 3D EM map. We discuss a new algorithm and generalized framework, named PF(2) fit (Polar Fast Fourier Fitting) for the best possible structural alignment of atomistic structures with 3D EM. While PF(2) fit enables only a rigid, six dimensional (6D) alignment method, it augments prior work on 6D X-ray structure and 3D EM alignment in multiple ways: Scoring. PF(2) fit includes a new scoring scheme that, in addition to rewarding overlaps between the volumes occupied by the atomistic structure and 3D EM map, rewards overlaps between the volumes complementary to them. We quantitatively demonstrate how this new complementary scoring scheme improves upon existing approaches. PF(2) fit also includes two scoring functions, the non-uniform exterior penalty and the skeleton-secondary structure score, and implements the scattering potential score as an alternative to traditional Gaussian blurring. Search. PF(2) fit utilizes a fast polar Fourier search scheme, whose main advantage is the ability to search over uniformly and adaptively sampled subsets of the space of rigid-body motions. PF(2) fit also implements a new reranking search and scoring methodology that considerably improves alignment metrics in results obtained from the initial search.

TEMPy: a Python library for assessment of three-dimensional electron microscopy density fits

TEMPy is an object-oriented Python library that provides the means to validate density fits in electron microscopy reconstructions. This article highlights several features of particular interest for this purpose and includes some customized examples.

Three-dimensional electron microscopy is currently one of the most promising techniques used to study macromolecular assemblies. Rigid and flexible fitting of atomic models into density maps is often essential to gain further insights into the assemblies they represent. Currently, tools that facilitate the assessment of fitted atomic models and maps are needed. TEMPy (template and electron microscopy comparison using Python) is a toolkit designed for this purpose. The library includes a set of methods to assess density fits in intermediate-to-low resolution maps, both globally and locally. It also provides procedures for single-fit assessment, ensemble generation of fits, clustering, and multiple and consensus scoring, as well as plots and output files for visualization purposes to help the user in analysing rigid and flexible fits. The modular nature of TEMPy helps the integration of scoring and assessment of fits into large pipelines, making it a tool suitable for both novice and expert structural biologists.

Building a pseudo-atomic model of the anaphase-promoting complex

The anaphase-promoting complex (APC/C) is a large E3 ubiquitin ligase that regulates progression through specific stages of the cell cycle by coordinating the ubiquitin-dependent degradation of cell-cycle regulatory proteins. Depending on the species, the active form of the APC/C consists of 14-15 different proteins that assemble into a 20-subunit complex with a mass of approximately 1.3 MDa. A hybrid approach of single-particle electron microscopy and protein crystallography of individual APC/C subunits has been applied to generate pseudo-atomic models of various functional states of the complex. Three approaches for assigning regions of the EM-derived APC/C density map to specific APC/C subunits are described. This information was used to dock atomic models of APC/C subunits, determined either by protein crystallography or homology modelling, to specific regions of the APC/C EM map, allowing the generation of a pseudo-atomic model corresponding to 80% of the entire complex.

iMODFIT: efficient and robust flexible fitting based on vibrational analysis in internal coordinates

Here, we employed the collective motions extracted from Normal Mode Analysis (NMA) in internal coordinates (torsional space) for the flexible fitting of atomic-resolution structures into electron microscopy (EM) density maps. The proposed methodology was validated using a benchmark of simulated cases, highlighting its robustness over the full range of EM resolutions and even over coarse-grained representations. A systematic comparison with other methods further showcased the advantages of this proposed methodology, especially at medium to lower resolutions. Using this method, computational costs and potential overfitting problems are naturally reduced by constraining the search in low-frequency NMA space, where covalent geometry is implicitly maintained. This method also effectively captures the macromolecular changes of a representative set of experimental test cases. We believe that this novel approach will extend the currently available EM hybrid methods to the atomic-level interpretation of large conformational changes and their functional implications.

Enhanced sampling and overfitting analyses in structural refinement of nucleic acids into electron microscopy maps

Flexible fitting computational algorithms are often useful to interpret low-resolution maps of many macromolecular complexes generated by electron microscopy (EM) imaging. One such atomistic simulation technique is molecular dynamics flexible fitting (MDFF), which has been widely applied to generate structural models of large ribonucleoprotein assemblies such as the ribosome. We have previously shown that MDFF simulations of globular proteins are sensitive to the resolution of the target EM map and the strength of restraints used to preserve the secondary structure elements during fitting (Vashisth, H.; et al. Structure 2012, 20, 1453-1462). In this work, we aim to systematically examine the quality of structural models of various nucleic acids obtained via MDFF by varying the map resolution and the strength of structural restraints. We also demonstrate how an enhanced conformational sampling technique for proteins, temperature-accelerated molecular dynamics (TAMD), can be combined with MDFF for the structural refinement of nucleic acids in EM maps. Finally, we also demonstrate application of TAMD-assisted MDFF (TAMDFF) on a RNA/protein complex and suggest that TAMDFF is a viable strategy for enhanced conformational fitting in target maps of ribonucleoprotein complexes.

Using enhanced sampling and structural restraints to refine atomic structures into low-resolution electron microscopy maps

For a variety of problems in structural biology, low-resolution maps generated by electron microscopy imaging are often interpreted with the help of various flexible-fitting computational algorithms. In this work, we systematically analyze the quality of final models of various proteins obtained via molecular dynamics flexible fitting (MDFF) by varying the map-resolution, strength of structural restraints, and the steering forces. We find that MDFF can be extended to understand conformational changes in lower-resolution maps if larger structural restraints and lower steering forces are used to prevent overfitting. We further show that the capabilities of MDFF can be extended by combining it with an enhanced conformational sampling method, temperature-accelerated molecular dynamics (TAMD). Specifically, either TAMD can be used to generate better starting configurations for MDFF fitting or TAMD-assisted MDFF (TAMDFF) can be performed to accelerate conformational search in atomistic simulations.

Outcome of the first electron microscopy validation task force meeting

This Meeting Review describes the proceedings and conclusions from the inaugural meeting of the Electron Microscopy Validation Task Force organized by the Unified Data Resource for 3DEM (http://www.emdatabank.org) and held at Rutgers University in New Brunswick, NJ on September 28 and 29, 2010. At the workshop, a group of scientists involved in collecting electron microscopy data, using the data to determine three-dimensional electron microscopy (3DEM) density maps, and building molecular models into the maps explored how to assess maps, models, and other data that are deposited into the Electron Microscopy Data Bank and Protein Data Bank public data archives. The specific recommendations resulting from the workshop aim to increase the impact of 3DEM in biology and medicine.

Conventions and workflows for using Situs

Situs is a modular program package for the multi-scale modeling of atomic resolution structures and low-resolution biophysical data from electron microscopy, tomography or small-angle X-ray scattering. This article provides an overview of recent developments in the Situs package, with an emphasis on workflows and conventions that are important for practical applications. The modular design of the programs facilitates scripting in the bash shell that allows specific programs to be combined in creative ways that go beyond the original intent of the developers. Several scripting-enabled functionalities, such as flexible transformations of data type, the use of symmetry constraints or the creation of two-dimensional projection images, are described. The processing of low-resolution biophysical maps in such workflows follows not only first principles but often relies on implicit conventions. Situs conventions related to map formats, resolution, correlation functions and feature detection are reviewed and summarized. The compatibility of the Situs workflow with CCP4 conventions and programs is discussed.

Macromolecular assembly structures by comparative modeling and electron microscopy

Advances in electron microscopy allow for structure determination of large biological machines at increasingly higher resolutions. A key step in this process is fitting component structures into the electron microscopy-derived density map of their assembly. Comparative modeling can contribute by providing atomic models of the components, via fold assignment, sequence-structure alignment, model building, and model assessment. All four stages of comparative modeling can also benefit from consideration of the density map. In this chapter, we describe numerous types of modeling problems restrained by a density map and available protocols for finding solutions. In particular, we provide detailed instructions for density map-guided modeling using the Integrative Modeling Platform (IMP), MODELLER, and UCSF Chimera.

Symmetry-restrained flexible fitting for symmetric EM maps

Many large biological macromolecules have inherent structural symmetry, being composed of a few distinct subunits, repeated in a symmetric array. These complexes are often not amenable to traditional high-resolution structural determination methods, but can be imaged in functionally relevant states using cryo-electron microscopy (cryo-EM). A number of methods for fitting atomic-scale structures into cryo-EM maps have been developed, including the molecular dynamics flexible fitting (MDFF) method. However, quality and resolution of the cryo-EM map are the major determinants of a method's success. In order to incorporate knowledge of structural symmetry into the fitting procedure, we developed the symmetry-restrained MDFF method. The new method adds to the cryo-EM map-derived potential further restraints on the allowed conformations of a complex during fitting, thereby improving the quality of the resultant structure. The benefit of using symmetry-based restraints during fitting, particularly for medium to low-resolution data, is demonstrated for three different systems.

UCSF Chimera, MODELLER, and IMP: an integrated modeling system

Structural modeling of macromolecular complexes greatly benefits from interactive visualization capabilities. Here we present the integration of several modeling tools into UCSF Chimera. These include comparative modeling by MODELLER, simultaneous fitting of multiple components into electron microscopy density maps by IMP MultiFit, computing of small-angle X-ray scattering profiles and fitting of the corresponding experimental profile by IMP FoXS, and assessment of amino acid sidechain conformations based on rotamer probabilities and local interactions by Chimera.

BCL::EM-Fit: rigid body fitting of atomic structures into density maps using geometric hashing and real space refinement

Cryo-electron microscopy (cryoEM) can visualize large macromolecular assemblies at resolutions often below 10Å and recently as good as 3.8-4.5 Å. These density maps provide important insights into the biological functioning of molecular machineries such as viruses or the ribosome, in particular if atomic-resolution crystal structures or models of individual components of the assembly can be placed into the density map. The present work introduces a novel algorithm termed BCL::EM-Fit that accurately fits atomic-detail structural models into medium resolution density maps. In an initial step, a "geometric hashing" algorithm provides a short list of likely placements. In a follow up Monte Carlo/Metropolis refinement step, the initial placements are optimized by their cross correlation coefficient. The resolution of density maps for a reliable fit was determined to be 10 Å or better using tests with simulated density maps. The algorithm was applied to fitting of capsid proteins into an experimental cryoEM density map of human adenovirus at a resolution of 6.8 and 9.0 Å, and fitting of the GroEL protein at 5.4 Å. In the process, the handedness of the cryoEM density map was unambiguously identified. The BCL::EM-Fit algorithm offers an alternative to the established Fourier/Real space fitting programs. BCL::EM-Fit is free for academic use and available from a web server or as downloadable binary file at http://www.meilerlab.org.

Exploring the spatial and temporal organization of a cell's proteome

To increase our current understanding of cellular processes, such as cell signaling and division, knowledge is needed about the spatial and temporal organization of the proteome at different organizational levels. These levels cover a wide range of length and time scales: from the atomic structures of macromolecules for inferring their molecular function, to the quantitative description of their abundance, and spatial distribution in the cell. Emerging new experimental technologies are greatly increasing the availability of such spatial information on the molecular organization in living cells. This review addresses three fields that have significantly contributed to our understanding of the proteome's spatial and temporal organization: first, methods for the structure determination of individual macromolecular assemblies, specifically the fitting of atomic structures into density maps generated from electron microscopy techniques; second, research that visualizes the spatial distributions of these complexes within the cellular context using cryo electron tomography techniques combined with computational image processing; and third, methods for the spatial modeling of the dynamic organization of the proteome, specifically those methods for simulating reaction and diffusion of proteins and complexes in crowded intracellular fluids. The long-term goal is to integrate the varied data about a proteome's organization into a spatially explicit, predictive model of cellular processes.

Scoring functions for cryoEM density fitting

In fitting atomic structures into cryoEM density maps of macromolecular assemblies, the cross-correlation function (CCF) is the most prevalent method of scoring the goodness-of-fit. However, there are still many possible, less studied ways of scoring fits. In this paper, we introduce four scores new to cryoEM fitting and compare their performance to three known scores. Our benchmark consists of (a) 4 protein assemblies with simulated maps at 5-20 Å resolution, including the heptameric ring of GroEL; and (b) 4 experimental maps of GroEL at ∼6-23 Å resolution with corresponding fitted atomic models. We perturb each fit 1000 times and assess each new fit with each score. The correlation between a score and the Cα RMSD of each fit from the "correctly" fitted structure shows that the CCF is one of the best scores, but in certain situations could be augmented or even replaced by other scores. For instance, our implementation of a score based on mutual information outperforms or is comparable to the CCF in almost all test cases, and our new "envelope score" works as well as the CCF at sub-nanometer resolution but is an order of magnitude faster to calculate. The results also suggest that the width of the Gaussian function used to blur the atomic structure into a density map can significantly affect the fitting process. Finally, we show that our score-testing method, when combined with the Laplacian CCF or the mutual information scores, can be used as a statistical tool for improving cryoEM density fitting.

Dimeric SecA couples the preprotein translocation in an asymmetric manner

The Sec translocase mediates the post-translational translocation of a number of preproteins through the inner membrane in bacteria. In the initiatory translocation step, SecB targets the preprotein to the translocase by specific interaction with its receptor SecA. The latter is the ATPase of Sec translocase which mediates the post-translational translocation of preprotein through the protein-conducting channel SecYEG in the bacterial inner membrane. We examined the structures of Escherichia coli Sec intermediates in solution as visualized by negatively stained electron microscopy in order to probe the oligomeric states of SecA during this process. The symmetric interaction pattern between the SecA dimer and SecB becomes asymmetric in the presence of proOmpA, and one of the SecA protomers predominantly binds to SecB/proOmpA. Our results suggest that during preprotein translocation, the two SecA protomers are different in structure and may play different roles.

Using Sculptor and Situs for simultaneous assembly of atomic components into low-resolution shapes

We describe an integrated software system called Sculptor that combines visualization capabilities with molecular modeling algorithms for the analysis of multi-scale data sets. Sculptor features extensive special purpose visualization techniques that are based on modern GPU programming and are capable of representing complex molecular assemblies in real-time. The integration of graphics and modeling offers several advantages. The user interface not only eases the usually steep learning curve of pure algorithmic techniques, but it also permits instant analysis and post-processing of results, as well as the integration of results from external software. Here, we implemented an interactive peak-selection strategy that enables the user to explore a preliminary score landscape generated by the colors tool of Situs. The interactive placement of components, one at a time, is advantageous for low-resolution or ambiguously shaped maps, which are sometimes difficult to interpret by the fully automatic peak selection of colors. For the subsequent refinement of the preliminary models resulting from both interactive and automatic peak selection, we have implemented a novel simultaneous multi-body docking in Sculptor and Situs that softly enforces shape complementarities between components using the normalization of the cross-correlation coefficient. The proposed techniques are freely available in Situs version 2.6 and Sculptor version 2.0.

Matching structural densities from different biophysical origins with gain and bias

The registration of volumetric structures in real space involves geometric and density transformations that align a target map and a probe map in the best way possible. Many computational docking strategies exist for finding the geometric transformations that superimpose maps, but the problem of finding an optimal density transformation, for the purposes of difference calculations or segmentation, has received little attention in the literature. We report results based on simulated and experimental electron microscopy maps, showing that a single scale factor (gain) may be insufficient when it comes to minimizing the density discrepancy between an aligned target and probe. We propose an affine transformation, with gain and bias, that is parameterized by known surface isovalues and by an interactive centering of the "cancellation peak" in the surface thresholded difference map histogram. The proposed approach minimizes discrepancies across a wide range of interior densities. Owing to having only two parameters, it avoids overfitting and requires only minimal knowledge of the probe and target maps. The linear transformation also preserves phases and relative amplitudes in Fourier space. The histogram matching strategy was implemented in the newly revised volhist tool of the Situs package, version 2.6.