Combining X-ray crystallography and electron microscopy

Docking strategies for predicting protein-ligand interactions and their application to structure-based drug design

Molecular docking stands as a pivotal element in the realm of computer-aided drug design (CADD), consistently contributing to advancements in pharmaceutical research. In essence, it employs computer algorithms to identify the "best" match between two molecules, akin to solving intricate three-dimensional jigsaw puzzles. At a more stringent level, the molecular docking challenge entails predicting the accurate bound association state based on the atomic coordinates of two molecules. This process assumes particular significance in unraveling the mechanistic intricacies of physicochemical interactions at the atomic scale. Notably, the application of docking, especially in the context of protein-small molecule interactions, holds wide-ranging implications for structure-based drug design, given the prevalent use of small compounds as drug candidates. This study provides an overview of docking methodologies, delves into recent key developments, elucidates the physicochemical underpinnings of molecular recognition in protein-ligand interactions, and concludes by addressing the applications of docking in virtual screening, alongside current challenges within existing docking methods.

Conventional Electron Microscopy, Cryogenic Electron Microscopy, and Cryogenic Electron Tomography of Viruses

Electron microscopy (EM) techniques have been crucial for understanding the structure of biological specimens such as cells, tissues and macromolecular assemblies. Viruses and related viral assemblies are ideal targets for structural studies that help to define essential biological functions. Whereas conventional EM methods use chemical fixation, dehydration, and staining of the specimens, cryogenic electron microscopy (cryo-EM) preserves the native hydrated state. Combined with image processing and three-dimensional reconstruction techniques, cryo-EM provides three-dimensional maps of these macromolecular complexes from projection images, at atomic or near-atomic resolutions. Cryo-EM is also a major technique in structural biology for dynamic studies of functional complexes, which are often unstable, flexible, scarce, or transient in their native environments. State-of-the-art techniques in structural virology now extend beyond purified symmetric capsids and focus on the asymmetric elements such as the packaged genome and minor structural proteins that were previously missed. As a tool, cryo-EM also complements high-resolution techniques such as X-ray diffraction and NMR spectroscopy; these synergistic hybrid approaches provide important new information. Three-dimensional cryogenic electron tomography (cryo-ET), a variation of cryo-EM, goes further, and allows the study of pleomorphic and complex viruses not only in their physiological state but also in their natural environment in the cell, thereby bridging structural studies at the molecular and cellular levels. Cryo-EM and cryo-ET have been applied successfully in basic research, shedding light on fundamental aspects of virus biology and providing insights into threatening viruses, including SARS-CoV-2, responsible for the COVID-19 pandemic.

Temperature-induced DNA density transition in phage λ capsid revealed with contrast-matching SANS

Structural details of a genome packaged in a viral capsid are essential for understanding how the structural arrangement of a viral genome in a capsid controls its release dynamics during infection, which critically affects viral replication. We previously found a temperature-induced, solid-like to fluid-like mechanical transition of packaged λ-genome that leads to rapid DNA ejection. However, an understanding of the structural origin of this transition was lacking. Here, we use small-angle neutron scattering (SANS) to reveal the scattering form factor of dsDNA packaged in phage λ capsid by contrast matching the scattering signal from the viral capsid with deuterated buffer. We used small-angle X-ray scattering and cryoelectron microscopy reconstructions to determine the initial structural input parameters for intracapsid DNA, which allows accurate modeling of our SANS data. As result, we show a temperature-dependent density transition of intracapsid DNA occurring between two coexisting phases-a hexagonally ordered high-density DNA phase in the capsid periphery and a low-density, less-ordered DNA phase in the core. As the temperature is increased from 20 °C to 40 °C, we found that the core-DNA phase undergoes a density and volume transition close to the physiological temperature of infection (~37 °C). The transition yields a lower energy state of DNA in the capsid core due to lower density and reduced packing defects. This increases DNA mobility, which is required to initiate rapid genome ejection from the virus capsid into a host cell, causing infection. These data reconcile our earlier findings of mechanical DNA transition in phage.

Novel Artificial Intelligence-Based Approaches for Ab Initio Structure Determination and Atomic Model Building for Cryo-Electron Microscopy

Single particle cryo-electron microscopy (cryo-EM) has emerged as the prevailing method for near-atomic structure determination, shedding light on the important molecular mechanisms of biological macromolecules. However, the inherent dynamics and structural variability of biological complexes coupled with the large number of experimental images generated by a cryo-EM experiment make data processing nontrivial. In particular, ab initio reconstruction and atomic model building remain major bottlenecks that demand substantial computational resources and manual intervention. Approaches utilizing recent innovations in artificial intelligence (AI) technology, particularly deep learning, have the potential to overcome the limitations that cannot be adequately addressed by traditional image processing approaches. Here, we review newly proposed AI-based methods for ab initio volume generation, heterogeneous 3D reconstruction, and atomic model building. We highlight the advancements made by the implementation of AI methods, as well as discuss remaining limitations and areas for future development.

Imaging Infection Across Scales of Size: From Whole Animals to Single Molecules

Infectious diseases are a leading cause of global morbidity and mortality, and the threat of infectious diseases to human health is steadily increasing as new diseases emerge, existing diseases reemerge, and antimicrobial resistance expands. The application of imaging technology to the study of infection biology has the potential to uncover new factors that are critical to the outcome of host-pathogen interactions and to lead to innovations in diagnosis and treatment of infectious diseases. This article reviews current and future opportunities for the application of imaging to the study of infectious diseases, with a particular focus on the power of imaging objects across a broad range of sizes to expand the utility of these approaches. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

A viral genome packaging motor transitions between cyclic and helical symmetry to translocate dsDNA

Molecular segregation and biopolymer manipulation require the action of molecular motors to do work by applying directional forces to macromolecules. The additional strand conserved E (ASCE) ring motors are an ancient family of molecular motors responsible for diverse biological polymer manipulation tasks. Viruses use ASCE segregation motors to package their genomes into their protein capsids and provide accessible experimental systems due to their relative simplicity. We show by cryo-EM-focused image reconstruction that ASCE ATPases in viral double-stranded DNA (dsDNA) packaging motors adopt helical symmetry complementary to their dsDNA substrates. Together with previous data, our results suggest that these motors cycle between helical and planar configurations, providing a possible mechanism for directional translocation of DNA. Similar changes in quaternary structure have been observed for proteasome and helicase motors, suggesting an ancient and common mechanism of force generation that has been adapted for specific tasks over the course of evolution.

Analysis of the Dynamic Proteasome Structure by Cross-Linking Mass Spectrometry

The 26S proteasome is a macromolecular complex that degrades proteins maintaining cell homeostasis; thus, determining its structure is a priority to understand its function. Although the 20S proteasome's structure has been known for some years, the highly dynamic nature of the 19S regulatory particle has presented a challenge to structural biologists. Advances in cryo-electron microscopy (cryo-EM) made it possible to determine the structure of the 19S regulatory particle and showed at least seven different conformational states of the proteasome. However, there are still many questions to be answered. Cross-linking mass spectrometry (CLMS) is now routinely used in integrative structural biology studies, and it promises to take integrative structural biology to the next level, answering some of these questions.

Bursting out: linking changes in nanotopography and biomechanical properties of biofilm-forming Escherichia coli to the T4 lytic cycle

The bacteriophage infection cycle has been extensively studied, yet little is known about the nanostructure and mechanical changes that lead to bacterial lysis. Here, atomic force microscopy was used to study in real time and in situ the impact of the canonical phage T4 on the nanotopography and biomechanics of irreversibly attached, biofilm-forming E. coli cells. The results show that in contrast to the lytic cycle in planktonic cells, which ends explosively, anchored cells that are in the process of forming a biofilm undergo a more gradual lysis, developing distinct nanoscale lesions (~300 nm in diameter) within the cell envelope. Furthermore, it is shown that the envelope rigidity and cell elasticity decrease (>50% and >40%, respectively) following T4 infection, a process likely linked to changes in the nanostructure of infected cells. These insights show that the well-established lytic pathway of planktonic cells may be significantly different from that of biofilm-forming cells. Elucidating the lysis paradigm of these cells may advance biofilm removal and phage therapeutics.

Haruspex: A Neural Network for the Automatic Identification of Oligonucleotides and Protein Secondary Structure in Cryo-Electron Microscopy Maps

In recent years, three-dimensional density maps reconstructed from single particle images obtained by electron cryo-microscopy (cryo-EM) have reached unprecedented resolution. However, map interpretation can be challenging, in particular if the constituting structures require de-novo model building or are very mobile. Herein, we demonstrate the potential of convolutional neural networks for the annotation of cryo-EM maps: our network Haruspex has been trained on a carefully curated set of 293 experimentally derived reconstruction maps to automatically annotate RNA/DNA as well as protein secondary structure elements. It can be straightforwardly applied to newly reconstructed maps in order to support domain placement or as a starting point for main-chain placement. Due to its high recall and precision rates of 95.1 % and 80.3 %, respectively, on an independent test set of 122 maps, it can also be used for validation during model building. The trained network will be available as part of the CCP-EM suite.

Anticipating innovations in structural biology

In this review, we describe how the interplay among science, technology and community interests contributed to the evolution of four structural biology data resources. We present the method by which data deposited by scientists are prepared for worldwide distribution, and argue that data archiving in a trusted repository must be an integral part of any scientific investigation.

Cryo-EM reconstruction of the Cafeteria roenbergensis virus capsid suggests novel assembly pathway for giant viruses

Whereas the protein composition and overall shape of several giant virus capsids have been described, the mechanism by which these large capsids assemble remains enigmatic. Here, we present a reconstruction of the capsid of Cafeteria roenbergensis virus (CroV), one of the largest viruses analyzed by cryo-electron microscopy (cryo-EM) to date. The CroV capsid has a diameter of 3,000 Å and a Triangulation number of 499. Unlike related mimiviruses, the CroV capsid is not decorated with glycosylated surface fibers, but features 30 Å-long surface protrusions that are formed by loops of the major capsid protein. Based on the orientation of capsomers in the cryo-EM reconstruction, we propose that the capsids of CroV and related giant viruses are assembled by a newly conceived assembly pathway that initiates at a five-fold vertex and continuously proceeds outwards in a spiraling fashion.

Structural and mechanistic insights into human splicing factor SF3b complex derived using an integrated approach guided by the cryo-EM density maps

Pre-mRNA splicing in eukaryotes is performed by the spliceosome, a highly complex macromolecular machine. SF3b is a multi-protein complex which recognizes the branch point adenosine of pre-mRNA as part of a larger U2 snRNP or U11/U12 di-snRNP in the dynamic spliceosome machinery. Although a cryo-EM map is available for human SF3b complex, the structure and relative spatial arrangement of all components in the complex are not yet known. We have recognized folds of domains in various proteins in the assembly and generated comparative models. Using an integrative approach involving structural and other experimental data, guided by the available cryo-EM density map, we deciphered a pseudo-atomic model of the closed form of SF3b which is found to be a "fuzzy complex" with highly flexible components and multiplicity of folds. Further, the model provides structural information for 5 proteins (SF3b10, SF3b155, SF3b145, SF3b130 and SF3b14b) and localization information for 4 proteins (SF3b10, SF3b145, SF3b130 and SF3b14b) in the assembly for the first time. Integration of this model with the available U11/U12 di-snRNP cryo-EM map enabled elucidation of an open form. This now provides new insights on the mechanistic features involved in the transition between closed and open forms pivoted by a hinge region in the SF3b155 protein that also harbors cancer causing mutations. Moreover, the open form guided model of the 5' end of U12 snRNA, which includes the branch point duplex, shows that the architecture of SF3b acts as a scaffold for U12 snRNA: pre-mRNA branch point duplex formation with potential implications for branch point adenosine recognition fidelity.

De Novo modeling in cryo-EM density maps with Pathwalking

As electron cryo-microscopy (cryo-EM) can now frequently achieve near atomic resolution, accurate interpretation of these density maps in terms of atomistic detail has become paramount in deciphering macromolecular structure and function. However, there are few software tools for modeling protein structure from cryo-EM density maps in this resolution range. Here, we present an extension of our original Pathwalking protocol, which can automatically trace a protein backbone directly from a near-atomic resolution (3-6Å) density map. The original Pathwalking approach utilized a Traveling Salesman Problem solver for backbone tracing, but manual adjustment was still required during modeling. In the new version, human intervention is minimized and we provide a more robust approach for backbone modeling. This includes iterative secondary structure identification, termini detection and the ability to model multiple subunits without prior segmentation. Overall, the new Pathwalking procedure provides a more complete and robust tool for annotating protein structure function in near-atomic resolution density maps.

Refinement of atomic models in high resolution EM reconstructions using Flex-EM and local assessment

As the resolutions of Three Dimensional Electron Microscopic reconstructions of biological macromolecules are being improved, there is a need for better fitting and refinement methods at high resolutions and robust approaches for model assessment. Flex-EM/MODELLER has been used for flexible fitting of atomic models in intermediate-to-low resolution density maps of different biological systems. Here, we demonstrate the suitability of the method to successfully refine structures at higher resolutions (2.5-4.5Å) using both simulated and experimental data, including a newly processed map of Apo-GroEL. A hierarchical refinement protocol was adopted where the rigid body definitions are relaxed and atom displacement steps are reduced progressively at successive stages of refinement. For the assessment of local fit, we used the SMOC (segment-based Manders' overlap coefficient) score, while the model quality was checked using the Qmean score. Comparison of SMOC profiles at different stages of refinement helped in detecting regions that are poorly fitted. We also show how initial model errors can have significant impact on the goodness-of-fit. Finally, we discuss the implementation of Flex-EM in the CCP-EM software suite.

X-ray structure determination using low-resolution electron microscopy maps for molecular replacement

Structures of multisubunit macromolecular machines are primarily determined either by electron microscopy (EM) or by X-ray crystallography. In many cases, a structure for a complex can be obtained at low resolution (at a coarse level of detail) with EM and at a higher resolution (with finer detail) by X-ray crystallography. The integration of these two structural techniques is becoming increasingly important for the generation of atomic models of macromolecular complexes. A low-resolution EM image can be a powerful tool for obtaining the 'phase' information that is missing from an X-ray crystallography experiment; however, integration of EM and X-ray diffraction data has been technically challenging. Here we present a step-by-step protocol that explains how low-resolution EM maps can be placed in the crystallographic unit cell by molecular replacement, and how initial phases computed from the placed EM density are extended to high resolution by averaging maps over noncrystallographic symmetry. As the resolution gap between EM and X-ray crystallography continues to narrow, the use of EM maps to help with X-ray crystal structure determination, as described in this protocol, will become increasingly effective.

Transmission electron microscopy and the molecular structure of icosahedral viruses

The field of structural virology developed in parallel with methodological advances in X-ray crystallography and cryo-electron microscopy. At the end of the 1970s, crystallography yielded the first high resolution structure of an icosahedral virus, the T=3 tomato bushy stunt virus at 2.9Å. It took longer to reach near-atomic resolution in three-dimensional virus maps derived from electron microscopy data, but this was finally achieved, with the solution of complex icosahedral capsids such as the T=25 human adenovirus at ∼3.5Å. Both techniques now work hand-in-hand to determine those aspects of virus assembly and biology that remain unclear. This review examines the trajectory followed by EM imaging techniques in showing the molecular structure of icosahedral viruses, from the first two-dimensional negative staining images of capsids to the latest sophisticated techniques that provide high resolution three-dimensional data, or snapshots of the conformational changes necessary to complete the infectious cycle.

Structural comparison of four different antibodies interacting with human papillomavirus 16 and mechanisms of neutralization

Cryo-electron microscopy (cryo-EM) was used to solve the structures of human papillomavirus type 16 (HPV16) complexed with fragments of antibody (Fab) from three different neutralizing monoclonals (mAbs): H16.1A, H16.14J, and H263.A2. The structure-function analysis revealed predominantly monovalent binding of each Fab with capsid interactions that involved multiple loops from symmetry related copies of the major capsid protein. The residues identified in each Fab-virus interface map to a conformational groove on the surface of the capsomer. In addition to the known involvement of the FG and HI loops, the DE loop was also found to constitute the core of each epitope. Surprisingly, the epitope mapping also identified minor contributions by EF and BC loops. Complementary immunological assays included mAb and Fab neutralization. The specific binding characteristics of mAbs correlated with different neutralizing behaviors in pre- and post-attachment neutralization assays.

PCalign: a method to quantify physicochemical similarity of protein-protein interfaces

Background

Structural comparison of protein-protein interfaces provides valuable insights into the functional relationship between proteins, which may not solely arise from shared evolutionary origin. A few methods that exist for such comparative studies have focused on structural models determined at atomic resolution, and may miss out interesting patterns present in large macromolecular complexes that are typically solved by low-resolution techniques.

Results

We developed a coarse-grained method, PCalign, to quantitatively evaluate physicochemical similarities between a given pair of protein-protein interfaces. This method uses an order-independent algorithm, geometric hashing, to superimpose the backbone atoms of a given pair of interfaces, and provides a normalized scoring function, PC-score, to account for the extent of overlap in terms of both geometric and chemical characteristics. We demonstrate that PCalign outperforms existing methods, and additionally facilitates comparative studies across models of different resolutions, which are not accommodated by existing methods. Furthermore, we illustrate potential application of our method to recognize interesting biological relationships masked by apparent lack of structural similarity.

Conclusions

PCalign is a useful method in recognizing shared chemical and spatial patterns among protein-protein interfaces. It outperforms existing methods for high-quality data, and additionally facilitates comparison across structural models with different levels of details with proven robustness against noise.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0471-x) contains supplementary material, which is available to authorized users.

Locking and blocking the viral landscape of an alphavirus with neutralizing antibodies

Alphaviruses are serious, sometimes lethal human pathogens that belong to the family Togaviridae. The structures of human Venezuelan equine encephalitis virus (VEEV), an alphavirus, in complex with two strongly neutralizing antibody Fab fragments (F5 and 3B4C-4) have been determined using a combination of cryo-electron microscopy and homology modeling. We characterize these monoclonal antibody Fab fragments, which are known to abrogate VEEV infectivity by binding to the E2 (envelope) surface glycoprotein. Both of these antibody Fab fragments cross-link the surface E2 glycoproteins and therefore probably inhibit infectivity by blocking the conformational changes that are required for making the virus fusogenic. The F5 Fab fragment cross-links E2 proteins within one trimeric spike, whereas the 3B4C-4 Fab fragment cross-links E2 proteins from neighboring spikes. Furthermore, F5 probably blocks the receptor-binding site, whereas 3B4C-4 sterically hinders the exposure of the fusion loop at the end of the E2 B-domain.

IMPORTANCE

Alphaviral infections are transmitted mainly by mosquitoes. Venezuelan equine encephalitis virus (VEEV) is an alphavirus with a wide distribution across the globe. No effective vaccines exist for alphaviral infections. Therefore, a better understanding of VEEV and its associated neutralizing antibodies will help with the development of effective drugs and vaccines.

Finding the right fit: chiseling structures out of cryo-electron microscopy maps

Cryo-electron microscopy is a central tool for studying the architecture of macromolecular complexes at subnanometer resolution. Interpretation of an electron microscopy map requires its computational integration with data about the structure's components from all available sources, notably atomic models. Selecting a protocol for EM density-guided integrative structural modeling depends on the resolution and quality of the EM map as well as the available complimentary datasets. Here, we review rigid, flexible, and de novo integrative fitting into EM maps and provide guidelines and considerations for the design of modeling experiments. Finally, we discuss efforts towards establishing unified criteria for map and model assessment and validation.

Probing the structure of human protein disulfide isomerase by chemical cross-linking combined with mass spectrometry

Protein disulfide-isomerase (PDI) is a four-domain flexible protein that catalyzes the formation of disulfide bonds in the endoplasmic reticulum. Here we have analyzed native PDI purified from human placenta by chemical cross-linking followed by mass spectrometry (CXMS). In addition to PDI the sample contained soluble calnexin and ERp72. Extensive cross-linking was observed within the PDI molecule, both intra- and inter-domain, as well as between the different components in the mixture. The high sensitivity of the analysis in the current experiments, combined with a likely promiscuous interaction pattern of the involved proteins, revealed relatively densely populated cross-link heat maps. The established X-ray structure of the monomeric PDI could be confirmed; however, the dimer as presented in the existing models does not seem to be prevalent in solution as modeling on the observed cross-links revealed new models of dimeric PDI. The observed inter-protein cross-links confirmed the existence of a peptide binding area on calnexin that binds strongly both PDI and ERp72. On the other hand, interaction sites on PDI and ERp72 could not be uniquely identified, indicating a more non-specific interaction pattern.

BIOLOGICAL SIGNIFICANCE

The present work demonstrates the use of chemical cross-linking and mass spectrometry (CXMS) for the determination of a solution structure of natural human PDI and its interaction with the chaperones ERp72 and calnexin. The data shows that the dimeric structure of PDI may be more diverse than indicated by present models. We further observe that the temperature influences the cross-linking pattern of PDI, but this does not influence the overall folding pattern of the molecule.

Asymmetric perturbations of signalling oligomers

This review focuses on rapid and reversible noncovalent interactions for symmetric oligomers of signalling proteins. Symmetry mismatch, transient symmetry breaking and asymmetric perturbations via chemical (ligand binding) and physical (electric or mechanic) effects can initiate the signalling events. Advanced biophysical methods can reveal not only structural symmetries of stable membrane-bound signalling proteins but also asymmetric functional transition states. Relevant techniques amenable to distinguish between symmetric and asymmetric architectures are discussed including those with the capability of capturing low-populated transient conformational states. Typical examples of signalling proteins are overviewed for symmetry breaking in dimers (GPCRs, growth factor receptors, transcription factors); trimers (acid-sensing ion channels); tetramers (voltage-gated cation channels, ionotropic glutamate receptor, CNG and CHN channels); pentameric ligand-gated and mechanosensitive channels; higher order oligomers (gap junction channel, chaperonins, proteasome, virus capsid); as well as primary and secondary transporters. In conclusion, asymmetric perturbations seem to play important functional roles in a broad range of communicating networks.

Combined approaches to flexible fitting and assessment in virus capsids undergoing conformational change

Fitting of atomic components into electron cryo-microscopy (cryoEM) density maps is routinely used to understand the structure and function of macromolecular machines. Many fitting methods have been developed, but a standard protocol for successful fitting and assessment of fitted models has yet to be agreed upon among the experts in the field. Here, we created and tested a protocol that highlights important issues related to homology modelling, density map segmentation, rigid and flexible fitting, as well as the assessment of fits. As part of it, we use two different flexible fitting methods (Flex-EM and iMODfit) and demonstrate how combining the analysis of multiple fits and model assessment could result in an improved model. The protocol is applied to the case of the mature and empty capsids of Coxsackievirus A7 (CAV7) by flexibly fitting homology models into the corresponding cryoEM density maps at 8.2 and 6.1 Å resolution. As a result, and due to the improved homology models (derived from recently solved crystal structures of a close homolog – EV71 capsid – in mature and empty forms), the final models present an improvement over previously published models. In close agreement with the capsid expansion observed in the EV71 structures, the new CAV7 models reveal that the expansion is accompanied by ∼5° counterclockwise rotation of the asymmetric unit, predominantly contributed by the capsid protein VP1. The protocol could be applied not only to viral capsids but also to many other complexes characterised by a combination of atomic structure modelling and cryoEM density fitting.

Molecular replacement then and now

The `phase problem' in crystallography results from the inability to directly measure the phases of individual diffracted X-ray waves. While intensities are directly measured during data collection, phases must be obtained by other means. Several phasing methods are available (MIR, SAR, MAD, SAD and MR) and they all rely on the premise that phase information can be obtained if the positions of marker atoms in the unknown crystal structure are known. This paper is dedicated to the most popular phasing method, molecular replacement (MR), and represents a personal overview of the development, use and requirements of the methodology. The first description of noncrystallographic symmetry as a tool for structure determination was explained by Rossmann and Blow [Rossmann & Blow (1962), Acta Cryst. 15, 24-31]. The term `molecular replacement' was introduced as the name of a book in which the early papers were collected and briefly reviewed [Rossmann (1972), The Molecular Replacement Method. New York: Gordon & Breach]. Several programs have evolved from the original concept to allow faster and more sophisticated searches, including six-dimensional searches and brute-force approaches. While careful selection of the resolution range for the search and the quality of the data will greatly influence the outcome, the correct choice of the search model is probably still the main criterion to guarantee success in solving a structure using MR. Two of the main parameters used to define the `best' search model are sequence identity (25% or more) and structural similarity. Another parameter that may often be undervalued is the quality of the probe: there is clearly a relationship between the quality and the correctness of the chosen probe and its usefulness as a search model. Efforts should be made by all structural biologists to ensure that their deposited structures, which are potential search probes for future systems, are of the best possible quality.

From lows to highs: using low-resolution models to phase X-ray data

The study of virus structures has contributed to methodological advances in structural biology that are generally applicable (molecular replacement and noncrystallographic symmetry are just two of the best known examples). Moreover, structural virology has been instrumental in forging the more general concept of exploiting phase information derived from multiple structural techniques. This hybridization of structural methods, primarily electron microscopy (EM) and X-ray crystallography, but also small-angle X-ray scattering (SAXS) and nuclear magnetic resonance (NMR) spectroscopy, is central to integrative structural biology. Here, the interplay of X-ray crystallography and EM is illustrated through the example of the structural determination of the marine lipid-containing bacteriophage PM2. Molecular replacement starting from an ~13 Å cryo-EM reconstruction, followed by cycling density averaging, phase extension and solvent flattening, gave the X-ray structure of the intact virus at 7 Å resolution This in turn served as a bridge to phase, to 2.5 Å resolution, data from twinned crystals of the major coat protein (P2), ultimately yielding a quasi-atomic model of the particle, which provided significant insights into virus evolution and viral membrane biogenesis.

Computational methods for constructing protein structure models from 3D electron microscopy maps

Protein structure determination by cryo-electron microscopy (EM) has made significant progress in the past decades. Resolutions of EM maps have been improving as evidenced by recently reported structures that are solved at high resolutions close to 3Å. Computational methods play a key role in interpreting EM data. Among many computational procedures applied to an EM map to obtain protein structure information, in this article we focus on reviewing computational methods that model protein three-dimensional (3D) structures from a 3D EM density map that is constructed from two-dimensional (2D) maps. The computational methods we discuss range from de novo methods, which identify structural elements in an EM map, to structure fitting methods, where known high resolution structures are fit into a low-resolution EM map. A list of available computational tools is also provided.

Cross-validation in cryo-EM-based structural modeling

Single-particle cryo-EM is a powerful approach to determine the structure of large macromolecules and assemblies thereof in many cases at subnanometer resolution. It has become popular to refine or flexibly fit atomic models into density maps derived from cryo-EM experiments. These density maps are typically significantly lower in resolution than electron density maps obtained from X-ray diffraction experiments, such that the number of parameters that need to be determined is much larger than the number of experimental observables. Overfitting and misinterpretation of the density, thus, become a serious problem. For diffraction data, a cross-validation approach was introduced almost 20 y ago; however, no such approach has been described yet for structure refinement against cryo-EM density maps, although the overfitting problem is, because of the lower resolution, significantly larger. We present a cross-validation approach for real-space refinement against cryo-EM density maps in analogy to cross-validation typically used in crystallography. Our approach is able to detect overfitting and allows for optimizing the choice of restraints used in the refinement. The approach is shown on three protein structures with simulated data and experimental data of the rotavirus double-layer particle. Because cross-validation requires splitting the dataset into at least two independent sets, we further present an approach to quantify correlations between the structure factor sets. This analysis is also helpful for other cross-validation applications, such as refinements against diffraction data or 3D reconstructions of cryo-EM density maps.

Consensus among multiple approaches as a reliability measure for flexible fitting into cryo-EM data

Cryo-electron microscopy (cryo-EM) can provide low-resolution density maps of large macromolecular assemblies. As the number of structures deposited in the Protein Data Bank by fitting a high-resolution structure into a low-resolution cryo-EM map is increasing, there is a need to revise the protocols and improve the measures for fitting. A recent study suggested using a combination of multiple automated flexible fitting approaches to improve the interpretation of cryo-EM data. The current work further explores the use of multiple approaches by validating this "consensus" fitting approach and deriving a local reliability measure. Here four different flexible fitting approaches are applied for fitting an initial structure into a simulated density map of known target structure from a dataset of proteins. It is found that the models produced from different approaches often have a consensus in conformation and are also near to the target structure, whereas cases not showing consensus are away from the target. A high correlation is also observed between the RMSF profiles calculated with respect to the average and the target structures, which indicates that the relation between consensus and accuracy can also be extended to a per-residue level. Therefore, the RMSF among the fitted models is proposed as a local reliability measure, which can be used to assess the reliability of the fit at specific regions. Hence, we encourage the community to use consensus flexible fitting with different methods to report on local reliability of the resulting models and improve the interpretation of cryo-EM data.

Conventional electron microscopy, cryo-electron microscopy and cryo-electron tomography of viruses

Electron microscopy (EM) techniques have been crucial for understanding the structure of biological specimens such as cells, tissues and macromolecular assemblies. Viruses and related viral assemblies are ideal targets for structural studies that help to define essential biological functions. Whereas conventional EM methods use chemical fixation, dehydration, and staining of the specimens, cryo-electron microscopy (cryo-EM) preserves the native hydrated state. Combined with image processing and three-dimensional reconstruction techniques, cryo-EM provides 3D maps of these macromolecular complexes from projection images, at subnanometer to near-atomic resolutions. Cryo-EM is also a major technique in structural biology for dynamic studies of functional complexes, which are often unstable, flexible, scarce or transient in their native environments. As a tool, cryo-EM complements high-resolution techniques such as X-ray diffraction and NMR spectroscopy; these synergistic hybrid approaches provide important new information. Three-dimensional cryo-electron tomography goes further, and allows the study of viruses not only in their physiological state, but also in their natural environment in the cell, thereby bridging structural studies at the molecular and cellular levels.

CIM(®) monolithic anion-exchange chromatography as a useful alternative to CsCl gradient purification of bacteriophage particles

The use of anion-exchange chromatography was investigated as an alternative method to concentrate and purify bacterial viruses, and parameters for different bacteriophages were compared. Chromatography was performed with Convective Interactive Media(®) monoliths, with three different volumes and two matrix chemistries. Eleven morphologically distinct phages were tested, infecting five different bacterial species. For each of the phages tested, a protocol was optimized, including the choice of column chemistry, loading, buffer and elution conditions. The capacity and recovery of the phages on the columns varied considerably between phages. We conclude that anion-exchange chromatography with monoliths is a valid alternative to the more traditional CsCl purification, has upscaling advantages, but it requires more extensive optimization.

Gorgon and pathwalking: macromolecular modeling tools for subnanometer resolution density maps

The complex interplay of proteins and other molecules, often in the form of large transitory assemblies, are critical to cellular function. Today, X-ray crystallography and electron cryo-microscopy (cryo-EM) are routinely used to image these macromolecular complexes, though often at limited resolutions. Despite the rapidly growing number of macromolecular structures, few tools exist for modeling and annotating structures in the range of 3-10 Å resolution. To address this need, we have developed a number of utilities specifically targeting subnanometer resolution density maps. As part of the 2010 Cryo-EM Modeling Challenge, we demonstrated two of our latest de novo modeling tools, Pathwalking and Gorgon, as well as a tool for secondary structure identification (SSEHunter) and a new rigid-body/flexible fitting tool in Gorgon. In total, we submitted 30 structural models from ten different subnanometer resolution data sets in four of the six challenge categories. Each of our utlities produced accurate structural models and annotations across the various density maps. In the end, the utilities that we present here offer users a robust toolkit for analyzing and modeling protein structure in macromolecular assemblies at non-atomic resolutions.

Near-atomic resolution structural model of the yeast 26S proteasome

The 26S proteasome operates at the executive end of the ubiquitin-proteasome pathway. Here, we present a cryo-EM structure of the Saccharomyces cerevisiae 26S proteasome at a resolution of 7.4 Å or 6.7 Å (Fourier-Shell Correlation of 0.5 or 0.3, respectively). We used this map in conjunction with molecular dynamics-based flexible fitting to build a near-atomic resolution model of the holocomplex. The quality of the map allowed us to assign α-helices, the predominant secondary structure element of the regulatory particle subunits, throughout the entire map. We were able to determine the architecture of the Rpn8/Rpn11 heterodimer, which had hitherto remained elusive. The MPN domain of Rpn11 is positioned directly above the AAA-ATPase N-ring suggesting that Rpn11 deubiquitylates substrates immediately following commitment and prior to their unfolding by the AAA-ATPase module. The MPN domain of Rpn11 dimerizes with that of Rpn8 and the C-termini of both subunits form long helices, which are integral parts of a coiled-coil module. Together with the C-terminal helices of the six PCI-domain subunits they form a very large coiled-coil bundle, which appears to serve as a flexible anchoring device for all the lid subunits.

Real-space refinement with DireX: from global fitting to side-chain improvements

Single-particle cryo-electron microscopy (cryo-EM) has become an important tool to determine the structure of large biomolecules and assemblies thereof. However, the achievable resolution varies considerably over a wide range of about 3.5-20 Å. The interpretation of these intermediate- to low-resolution density maps in terms of atomic models is a big challenge and an area of active research. Here, we present our real-space structure refinement program DireX, which was developed primarily for cryo-EM-derived density maps. The basic principle and its main features are described. DireX employs Deformable Elastic Network (DEN) restraints to reduce overfitting by decreasing the effective number of degrees of freedom used in the refinement. Missing or reduced density due to flexible parts of the protein can lead to artifacts in the structure refinement, which is addressed through the concept of restrained grouped occupancy refinement. Furthermore, we describe the performance of DireX in the 2010 Cryo-EM Modeling Challenge, where we chose six density maps of four different proteins provided by the Modeling Challenge exemplifying typical refinement results at a large resolution range from 3 to 23 Å.

Fitting multimeric protein complexes into electron microscopy maps using 3D Zernike descriptors

A novel computational method for fitting high-resolution structures of multiple proteins into a cryoelectron microscopy map is presented. The method named EMLZerD generates a pool of candidate multiple protein docking conformations of component proteins, which are later compared with a provided electron microscopy (EM) density map to select the ones that fit well into the EM map. The comparison of docking conformations and the EM map is performed using the 3D Zernike descriptor (3DZD), a mathematical series expansion of three-dimensional functions. The 3DZD provides a unified representation of the surface shape of multimeric protein complex models and EM maps, which allows a convenient, fast quantitative comparison of the three-dimensional structural data. Out of 19 multimeric complexes tested, near native complex structures with a root-mean-square deviation of less than 2.5 Å were obtained for 14 cases while medium range resolution structures with correct topology were computed for the additional 5 cases.

Atomic model of the type III secretion system needle

Pathogenic bacteria using a type III secretion system (T3SS) to manipulate host cells cause many different infections including Shigella dysentery, typhoid fever, enterohaemorrhagic colitis and bubonic plague. An essential part of the T3SS is a hollow needle-like protein filament through which effector proteins are injected into eukaryotic host cells. Currently, the three-dimensional structure of the needle is unknown because it is not amenable to X-ray crystallography and solution NMR, as a result of its inherent non-crystallinity and insolubility. Cryo-electron microscopy combined with crystal or solution NMR subunit structures has recently provided a powerful hybrid approach for studying supramolecular assemblies, resulting in low-resolution and medium-resolution models. However, such approaches cannot deliver atomic details, especially of the crucial subunit-subunit interfaces, because of the limited cryo-electron microscopic resolution obtained in these studies. Here we report an alternative approach combining recombinant wild-type needle production, solid-state NMR, electron microscopy and Rosetta modelling to reveal the supramolecular interfaces and ultimately the complete atomic structure of the Salmonella typhimurium T3SS needle. We show that the 80-residue subunits form a right-handed helical assembly with roughly 11 subunits per two turns, similar to that of the flagellar filament of S. typhimurium. In contrast to established models of the needle in which the amino terminus of the protein subunit was assumed to be α-helical and positioned inside the needle, our model reveals an extended amino-terminal domain that is positioned on the surface of the needle, while the highly conserved carboxy terminus points towards the lumen.

Crystallography, evolution, and the structure of viruses

My undergraduate education in mathematics and physics was a good grounding for graduate studies in crystallographic studies of small organic molecules. As a postdoctoral fellow in Minnesota, I learned how to program an early electronic computer for crystallographic calculations. I then joined Max Perutz, excited to use my skills in the determination of the first protein structures. The results were even more fascinating than the development of techniques and provided inspiration for starting my own laboratory at Purdue University. My first studies on dehydrogenases established the conservation of nucleotide-binding structures. Having thus established myself as an independent scientist, I could start on my most cherished ambition of studying the structure of viruses. About a decade later, my laboratory had produced the structure of a small RNA plant virus and then, in another six years, the first structure of a human common cold virus. Many more virus structures followed, but soon it became essential to supplement crystallography with electron microscopy to investigate viral assembly, viral infection of cells, and neutralization of viruses by antibodies. A major guide in all these studies was the discovery of evolution at the molecular level. The conservation of three-dimensional structure has been a recurring theme, from my experiences with Max Perutz in the study of hemoglobin to the recognition of the conserved nucleotide-binding fold and to the recognition of the jelly roll fold in the capsid protein of a large variety of viruses.

Visualization of cell microtubules in their native state

BACKGROUND INFORMATION

Over the past decades, cryo-electron microscopy of vitrified specimens has yielded a detailed understanding of the tubulin and microtubule structures of samples reassembled in vitro from purified components. However, our knowledge of microtubule structure in vivo remains limited by the chemical treatments commonly used to observe cellular architecture using electron microscopy.

RESULTS

We used cryo-electron microscopy and cryo-electron tomography of vitreous sections to investigate the ultrastructure of microtubules in their cellular context. Vitreous sections were obtained from organotypic slices of rat hippocampus and from Chinese-hamster ovary cells in culture. Microtubules revealed their protofilament ultrastructure, polarity and, in the most favourable cases, molecular details comparable with those visualized in three-dimensional reconstructions of microtubules reassembled in vitro from purified tubulin. The resolution of the tomograms was estimated to be approx. 4 nm, which enabled the detection of luminal particles of approx. 6 nm in diameter inside microtubules.

CONCLUSIONS

The present study provides a first step towards a description of microtubules, in addition to other macromolecular assemblies, in an unperturbed cellular context at the molecular level. As the resolution appears to be similar to that obtainable with plunge-frozen samples, it should allow for the in vivo identification of larger macromolecular assemblies in vitreous sections of whole cells and tissues.

Reconstructing virus structures from nanometer to near-atomic resolutions with cryo-electron microscopy and tomography

The past few decades have seen tremendous advances in single-particle electron -cryo-microscopy (cryo-EM). The field has matured to the point that near-atomic resolution density maps can be generated for icosahedral viruses without the need for crystallization. In parallel, substantial progress has been made in determining the structures of nonicosahedrally arranged proteins in viruses by employing either single-particle cryo-EM or cryo-electron tomography (cryo-ET). Implicit in this course have been the availability of a new generation of electron cryo-microscopes and the development of the computational tools that are essential for generating these maps and models. This methodology has enabled structural biologists to analyze structures in increasing detail for virus particles that are in different morphogenetic states. Furthermore, electron imaging of frozen, hydrated cells, in the process of being infected by viruses, has also opened up a new avenue for studying virus structures "in situ". Here we present the common techniques used to acquire and process cryo-EM and cryo-ET data and discuss their implications for structural virology both now and in the future.

Consensus among flexible fitting approaches improves the interpretation of cryo-EM data

Cryo-elecron microscopy (cryo-EM) can provide important structural information of large macromolecular assemblies in different conformational states. Recent years have seen an increase in structures deposited in the Protein Data Bank (PDB) by fitting a high-resolution structure into its low-resolution cryo-EM map. A commonly used protocol for accommodating the conformational changes between the X-ray structure and the cryo-EM map is rigid body fitting of individual domains. With the emergence of different flexible fitting approaches, there is a need to compare and revise these different protocols for the fitting. We have applied three diverse automated flexible fitting approaches on a protein dataset for which rigid domain fitting (RDF) models have been deposited in the PDB. In general, a consensus is observed in the conformations, which indicates a convergence from these theoretically different approaches to the most probable solution corresponding to the cryo-EM map. However, the result shows that the convergence might not be observed for proteins with complex conformational changes or with missing densities in cryo-EM map. In contrast, RDF structures deposited in the PDB can represent conformations that not only differ from the consensus obtained by flexible fitting but also from X-ray crystallography. Thus, this study emphasizes that a "consensus" achieved by the use of several automated flexible fitting approaches can provide a higher level of confidence in the modeled configurations. Following this protocol not only increases the confidence level of fitting, but also highlights protein regions with uncertain fitting. Hence, this protocol can lead to better interpretation of cryo-EM data.

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.

Determining macromolecular assembly structures by molecular docking and fitting into an electron density map

Structural models of macromolecular assemblies are instrumental for gaining a mechanistic understanding of cellular processes. Determining these structures is a major challenge for experimental techniques, such as X-ray crystallography, NMR spectroscopy and electron microscopy (EM). Thus, computational modeling techniques, including molecular docking, are required. The development of most molecular docking methods has so far been focused on modeling of binary complexes. We have recently introduced the MultiFit method for modeling the structure of a multisubunit complex by simultaneously optimizing the fit of the model into an EM density map of the entire complex and the shape complementarity between interacting subunits. Here, we report algorithmic advances of the MultiFit method that result in an efficient and accurate assembly of the input subunits into their density map. The successful predictions and the increasing number of complexes being characterized by EM suggests that the CAPRI challenge could be extended to include docking-based modeling of macromolecular assemblies guided by EM.

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.

Modeling protein structure at near atomic resolutions with Gorgon

Electron cryo-microscopy (cryo-EM) has played an increasingly important role in elucidating the structure and function of macromolecular assemblies in near native solution conditions. Typically, however, only non-atomic resolution reconstructions have been obtained for these large complexes, necessitating computational tools for integrating and extracting structural details. With recent advances in cryo-EM, maps at near-atomic resolutions have been achieved for several macromolecular assemblies from which models have been manually constructed. In this work, we describe a new interactive modeling toolkit called Gorgon targeted at intermediate to near-atomic resolution density maps (10-3.5 Å), particularly from cryo-EM. Gorgon's de novo modeling procedure couples sequence-based secondary structure prediction with feature detection and geometric modeling techniques to generate initial protein backbone models. Beyond model building, Gorgon is an extensible interactive visualization platform with a variety of computational tools for annotating a wide variety of 3D volumes. Examples from cryo-EM maps of Rotavirus and Rice Dwarf Virus are used to demonstrate its applicability to modeling protein structure.

Modeling pilus structures from sparse data

Bacterial Type II secretion systems (T2SS) and type IV pili (T4P) biogenesis machineries share the ability to assemble thin filaments from pilin protein subunits in the plasma membrane. Here we describe in detail the calculation strategy that served to determine a detailed atomic model of the T2SS pilus from Klebsiella oxytoca (Campos et al., PNAS 2010). The strategy is based on molecular modeling with generalized distance restraints and experimental validation (salt bridge charge inversion; double cysteine substitution and crosslinking). It does not require directly fitting structures into an envelope obtained from electron microscopy, but relies on lower resolution information, in particular the symmetry parameters of the helix forming the pilus. We validate the strategy with T4P where either a higher resolution structure is available (for the gonococcal (GC) pilus from Neisseria gonorrhoeae), or where we can compare our results to additional experimental data (for Vibrio cholerae TCP). The models are of sufficient precision to compare the architecture of the different pili in detail.