Molding Atomic Structures into Intermediate- Resolution Cryo-EM Density Maps of Ribosomal Complexes Using Real-Space Refinement

Single-Particle Reconstruction of Biological Molecules-Story in a Sample (Nobel Lecture)

Pictures tell a thousand words: The development of single-particle cryo-electron microscopy set the stage for high-resolution structure determination of biological molecules. In his Nobel lecture, J. Frank describes the ground-breaking discoveries that have enabled the development of cryo-EM. The method has taken biochemistry into a new era.

The joys and perils of flexible fitting

While performing their functions, biological macromolecules often form large, dynamically changing macromolecular assemblies. Only a relatively small number of such assemblies have been accessible to the atomic-resolution techniques X-ray crystallography and NMR. Electron microscopy in conjunction with image reconstruction has become the preferred alternative for revealing the structures of such macromolecular complexes. However, for most assemblies the achievable resolution is too low to allow accurate atomic modeling directly from the data. Yet, useful models often can be obtained by fitting atomic models of individual components into a low-resolution reconstruction of the entire assembly. Several algorithms for achieving optimal fits in this context were developed recently, many allowing considerable degrees of flexibility to account for binding-induced conformational changes of the assembly components. This chapter describes the advantages and potential pitfalls of these methods and puts them into perspective with alternative approaches such as iterative modular fitting of rigid-body domains.

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.

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 Å.

Putting structure into context: fitting of atomic models into electron microscopic and electron tomographic reconstructions

A complete understanding of complex dynamic cellular processes such as cell migration or cell adhesion requires the integration of atomic level structural information into the larger cellular context. While direct atomic-level information at the cellular level remains inaccessible, electron microscopy, electron tomography and their associated computational image processing approaches have now matured to a point where sub-cellular structures can be imaged in three dimensions at the nanometer scale. Atomic-resolution information obtained by other means can be combined with this data to obtain three-dimensional models of large macromolecular assemblies in their cellular context. This article summarizes some recent advances in this field.

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.

Accurate flexible fitting of high-resolution protein structures into cryo-electron microscopy maps using coarse-grained pseudo-energy minimization

Cryo-electron microscopy (cryo-EM) has been widely used to explore conformational states of large biomolecular assemblies. The detailed interpretation of cryo-EM data requires the flexible fitting of a known high-resolution protein structure into a low-resolution cryo-EM map. To this end, we have developed what we believe is a new method based on a two-bead-per-residue protein representation, and a modified form of the elastic network model that allows large-scale conformational changes while maintaining pseudobonds and secondary structures. Our method minimizes a pseudo-energy which linearly combines various terms of the modified elastic network model energy with a cryo-EM-fitting score and a collision energy that penalizes steric collisions. Unlike previous flexible fitting efforts using the lowest few normal modes, our method effectively utilizes all normal modes so that both global and local structural changes can be fully modeled. We have validated our method for a diverse set of 10 pairs of protein structures using simulated cryo-EM maps with a range of resolutions and in the absence/presence of random noise. We have shown that our method is both accurate and efficient compared with alternative techniques, and its performance is robust to the addition of random noise. Our method is also shown to be useful for the flexible fitting of three experimental cryo-EM maps.

Flexible fitting of high-resolution x-ray structures into cryoelectron microscopy maps using biased molecular dynamics simulations

A methodology for flexible fitting of all-atom high-resolution structures into low-resolution cryoelectron microscopy (cryo-EM) maps is presented. Flexibility of the modeled structure is simulated by classical molecular dynamics and an additional effective potential is introduced to enhance the fitting process. The additional potential is proportional to the correlation coefficient between the experimental cryo-EM map and a synthetic map generated for an all-atom structure being fitted to the map. The additional forces are calculated as a gradient of the correlation coefficient. During the molecular dynamics simulations under the additional forces, the molecule undergoes a conformational transition that maximizes the correlation coefficient, which results in a high-accuracy fit of all-atom structure into a cryo-EM map. Using five test proteins that exhibit structural rearrangement during their biological activity, we demonstrate performance of our method. We also test our method on the experimental cryo-EM of elongation factor G and show that the model obtained is comparable to previous studies. In addition, we show that overfitting can be avoided by assessing the quality of the fitted model in terms of correlation coefficient and secondary structure preservation.

Biomolecular pleiomorphism probed by spatial interpolation of coarse models

In low resolution structures of biological assemblies one can often observe conformational deviations that require a flexible rearrangement of structural domains fitted at the atomic level. We are evaluating interpolation methods for the flexible alignment of atomic models based on coarse models. Spatial interpolation is well established in image-processing and visualization to describe the overall deformation or warping of an object or an image. Combined with a coarse representation of the biological system by feature vectors, such methods can provide a flexible approximation of the molecular structure. We have compared three well-known interpolation techniques and evaluated the results by comparing them with constrained molecular dynamics. One method, inverse distance weighting interpolation, consistently produced models that were nearly indistinguishable on the alpha carbon level from the molecular dynamics results. The method is simple to apply and enables flexing of structures by non-expert modelers. This is useful for the basic interpretation of volumetric data in biological applications such as electron microscopy. The method can be used as a general interpretation tool for sparsely sampled motions derived from coarse models.

Multiple subunit fitting into a low-resolution density map of a macromolecular complex using a gaussian mixture model

Recently, electron microscopy measurement of single particles has enabled us to reconstruct a low-resolution 3D density map of large biomolecular complexes. If structures of the complex subunits can be solved by x-ray crystallography at atomic resolution, fitting these models into the 3D density map can generate an atomic resolution model of the entire large complex. The fitting of multiple subunits, however, generally requires large computational costs; therefore, development of an efficient algorithm is required. We developed a fast fitting program, “gmfit”, which employs a Gaussian mixture model (GMM) to represent approximated shapes of the 3D density map and the atomic models. A GMM is a distribution function composed by adding together several 3D Gaussian density functions. Because our model analytically provides an integral of a product of two distribution functions, it enables us to quickly calculate the fitness of the density map and the atomic models. Using the integral, two types of potential energy function are introduced: the attraction potential energy between a 3D density map and each subunit, and the repulsion potential energy between subunits. The restraint energy for symmetry is also employed to build symmetrical origomeric complexes. To find the optimal configuration of subunits, we randomly generated initial configurations of subunit models, and performed a steepest-descent method using forces and torques of the three potential energies. Comparison between an original density map and its GMM showed that the required number of Gaussian distribution functions for a given accuracy depended on both resolution and molecular size. We then performed test fitting calculations for simulated low-resolution density maps of atomic models of homodimer, trimer, and hexamer, using different search parameters. The results indicated that our method was able to rebuild atomic models of a complex even for maps of 30 Å resolution if sufficient numbers (eight or more) of Gaussian distribution functions were employed for each subunit, and the symmetric restraints were assigned for complexes with more than three subunits. As a more realistic test, we tried to build an atomic model of the GroEL/ES complex by fitting 21-subunit atomic models into the 3D density map obtained by cryoelectron microscopy using the C7 symmetric restraints. A model with low root mean-square deviations (14.7 Å) was obtained as the lowest-energy model, showing that our fitting method was reasonably accurate. Inclusion of other restraints from biological and biochemical experiments could further enhance the accuracy.

What transmission electron microscopes can visualize now and in the future

Our review concentrates on the progress made in high-resolution transmission electron microscopy (TEM) in the past decade. This includes significant improvements in sample preparation by quick-freezing aimed at preserving the specimen in a close-to-native state in the high vacuum of the microscope. Following advances in cold stage and TEM vacuum technology systems, the observation of native, frozen hydrated specimens has become a widely used approach. It fostered the development of computer guided, fully automated low-dose data acquisition systems allowing matched pairs of images and diffraction patterns to be recorded for electron crystallography, and the collection of entire tilt-series for electron tomography. To achieve optimal information transfer to atomic resolution, field emission electron guns combined with acceleration voltages of 200-300 kV are now routinely used. The outcome of these advances is illustrated by the atomic structure of mammalian aquaporin-O and by the pore-forming bacterial cytotoxin ClyA resolved to 12 A. Further, the Yersinia injectisome needle, a bacterial pseudopilus and the binding of phalloidin to muscle actin filaments were chosen to document the advantage of the high contrast offered by dedicated scanning transmission electron microscopy (STEM) and/or the STEM's ability to measure the mass of protein complexes and directly link this to their shape. Continued progress emerging from leading research laboratories and microscope manufacturers will eventually enable us to determine the proteome of a single cell by electron tomography, and to more routinely solve the atomic structure of membrane proteins by electron crystallography.

Cryo-EM structure of the DNA-dependent protein kinase catalytic subunit at subnanometer resolution reveals alpha helices and insight into DNA binding

The DNA-dependent protein kinase catalytic subunit (DNA-PKcs) regulates the nonhomologous end joining pathway for repair of double-stranded DNA (dsDNA) breaks. Here, we present a 7A resolution structure of DNA-PKcs determined by cryo-electron microscopy single-particle reconstruction. This structure is composed of density rods throughout the molecule that are indicative of alpha helices and reveals structural features not observed in lower resolution EM structures. Docking of homology models into the DNA-PKcs structure demonstrates that up to eight helical HEAT repeat motifs fit well within the density. Surprisingly, models for the kinase domain can be docked into either the crown or base of the molecule at this resolution, although real space refinement suggests that the base location is the best fit. We propose a model for the interaction of DNA with DNA-PKcs in which one turn of dsDNA enters the central channel and interacts with a resolved alpha-helical protrusion.

Combining efficient conformational sampling with a deformable elastic network model facilitates structure refinement at low resolution

Structural studies of large proteins and protein assemblies are a difficult and pressing challenge in molecular biology. Experiments often yield only low-resolution or sparse data that are not sufficient to fully determine atomistic structures. We have developed a general geometry-based algorithm that efficiently samples conformational space under constraints imposed by low-resolution density maps obtained from electron microscopy or X-ray crystallography experiments. A deformable elastic network (DEN) is used to restrain the sampling to prior knowledge of an approximate structure. The DEN restraints dramatically reduce over-fitting, especially at low resolution. Cross-validation is used to optimally weight the structural information and experimental data. Our algorithm is robust even for noise-added density maps and has a large radius of convergence for our test case. The DEN restraints can also be used to enhance reciprocal space simulated annealing refinement.

Three-dimensional visualization of FKBP12.6 binding to an open conformation of cardiac ryanodine receptor

The cardiac isoform of the ryanodine receptor (RyR2) from dog binds predominantly a 12.6-kDa isoform of the FK506-binding protein (FKBP12.6), whereas RyR2 from other species binds both FKBP12.6 and the closely related isoform FKBP12. The role played by FKBP12.6 in modulating calcium release by RyR2 is unclear at present. We have used cryoelectron microscopy and three-dimensional (3D) reconstruction techniques to determine the binding position of FKBP12.6 on the surface of canine RyR2. Buffer conditions that should favor the "open" state of RyR2 were used. Quantitative comparison of 3D reconstructions of RyR2 in the presence and absence of FKBP12.6 reveals that FKBP12.6 binds along the sides of the square-shaped cytoplasmic region of the receptor, adjacent to domain 9, which forms part of the four clamp (corner-forming) structures. The location of the FKBP12.6 binding site on "open" RyR2 appears similar, but slightly displaced (by 1-2 nm) from that found previously for FKBP12 binding to the skeletal muscle ryanodine receptor that was in the buffer that favors the "closed" state. The conformation of RyR2 containing bound FKBP12.6 differs considerably from that depleted of FKBP12.6, particularly in the transmembrane region and in the clamp structures. The x-ray structure of FKBP12.6 was docked into the region of the 3D reconstruction that is attributable to bound FKBP12.6, to show the relative orientations of amino acid residues (Gln-31, Asn-32, Phe-59) that have been implicated as being critical in interactions with RyR2. A thorough understanding of the structural basis of RyR2-FKBP12.6 interaction should aid in understanding the roles that have been proposed for FKBP12.6 in heart failure and in certain forms of sudden cardiac death.

Combining electron microscopy and comparative protein structure modeling

Recently, advances have been made in methods and applications that integrate electron microscopy density maps and comparative modeling to produce atomic structures of macromolecular assemblies. Electron microscopy can benefit from comparative modeling through the fitting of comparative models into electron microscopy density maps. Also, comparative modeling can benefit from electron microscopy through the use of intermediate-resolution density maps in fold recognition, template selection and sequence-structure alignment.