To obtain a structural model of a macromolecular assembly by single-particle EM, a large number of particle images need to be collected, aligned, clustered, averaged, and finally assembled via reconstruction into a 3D density map. This process is limited by the number and quality of the particle images, the accuracy of the initial model, and the compositional and conformational heterogeneity. Here, we describe a structure determination method that avoids the reconstruction procedure. The atomic structures of the individual complex components are assembled by optimizing a match against 2D EM class-average images, an excluded volume criterion, geometric complementarity, and optional restraints from proteomics and chemical cross-linking experiments. The optimization relies on a simulated annealing Monte Carlo search and a divide-and-conquer message-passing algorithm. Using simulated and experimentally determined EM class averages for 12 and 4 protein assemblies, respectively, we show that a few class averages can indeed result in accurate models for complexes of as many as five subunits. Thus, integrative structural biology can now benefit from the relative ease with which the EM class averages are determined.integrative modeling | structural determination | computational biology E M is an increasingly useful approach for structural characterization of macromolecular assemblies (1-3). Different flavors of EM include electron crystallography, single-particle EM, and electron tomography (4), although tomography is generally limited to resolutions worse than 20 Å. Single-particle EM can be used with negative-stained or cryogenically frozen (cryo-EM) samples. Cryo-EM particularly presents some attractive advantages: It preserves a near-native conformation of the molecules, it can be applied to the study of large assemblies, and it can theoretically achieve atomic resolution (5). Currently, the standard way to analyze single-particle EM data is to align collected 2D single-particle images, cluster the images, calculate a 2D classaverage image for each cluster, and finally perform 3D reconstruction to obtain a 3D density map. Pseudoatomic models for many assemblies have been generated by fitting X-ray crystallographic structures and/or comparative models of the individual components into a density map of the whole assembly (6).Obtaining a high-resolution density map requires a large number of single-particle images and depends critically on determining an initial low-resolution density map as a template for reconstruction, as well as a high signal-to-noise ratio (SNR) in the particle images (7-10). Several methods exist to obtain the template map: randomconical reconstruction (11), common lines determination (12)(13)(14), and maximization of the posterior probability of observing the set of class averages (15). A challenge for template construction arises when the available class averages do not provide significant coverage of all the orientations of the complex; in such a case, an accurate template map cannot be compute...