CC BY
21
RESEARCH ARTICLES
Three-Dimensional Structure оf Cytochrome c Nitrite Reductase As Determined by Cryo-Electron Microscopy
T. N. Baymukhametov1, Y. M. Chesnokov1, E. B. Pichkur1, K. M. Boyko12, T. V. Tikhonova2, A. G. Myasnikov3 4 5, A. L. Vasiliev16, A. V. Lipkin1, V. O. Popov12*, M. V. Kovalchuk16 National Research Center «Kurchatov Institute», Akademika Kurchatova Sqr., 1, Moscow, 123182, Russia
2Bach Institute of Biochemistry, Research Center of Biotechnology of the Russian Academy of Sciences, Leninsky Ave., 33, bldg. 2, Moscow, 119071, Russia
3Petersburg Nuclear Physics Institute named by B.P. Konstantinov of NRC «Kurchatov Institute», 188300, Leningradskaya Oblast, Gatchina, mkr. Orlova roshcha, 1, Russia
4Centre for Integrative Biology, Department of Integrated Structural Biology, Institute of Genetics and of Molecular and Cellular Biology, 67404, 1 rue Laurent Fries, Illkirch, France 5University of California San Francisco Mission Bay, Genentech Hall, San Francisco, CA, 941582517, USA
6Shubnikov Institute of Crystallography of Federal Scientific Research Centre «Crystallography and Photonics»
Russian Academy of Sciences, Leninsky Ave., 59, Moscow, 1 19333, Russia *E-mail: [email protected]
Received February 05, 2018; in final form June 02, 2018
Copyright © 2018 Park-media, Ltd. This is an open access article distributed under the Creative Commons Attribution License,which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
ABSTRACT The structure of cytochrome c nitrite reductase from the bacterium Thioalkalivibrio nitratireducens was determined by cryo-electron microscopy (cryo-EM) at a 2.56 A resolution. Possible structural heterogeneity of the enzyme was assessed. The backbone and side-chain orientations in the cryo-EM-based model are, in general, similar to those in the high-resolution X-ray diffraction structure of this enzyme.
KEYWORDS single particle analysis, high-resolution cryo-electron microscopy, structural biology, cytochrome c nitrite reductase.
ABBREVIATIONS TvNiR - cytochrome c nitrite reductase from Thioalkalivibrio nitratireducens, cryo-EM -cryo-electron microscopy, CMOS - complementary metal-oxide semiconductor, CTF - contrast transfer function, FSC - Fourier shell correlation.
INTRODUCTION
With the recent advances in cryo-electron microscopy (cryo-EM) associated with improvements in spatial resolution and a decrease in the lower limit of the molecular weight accessible through this technique, the method has begun to rival X-ray crystallography [1, 2]. Owing to the unique opportunity of extracting structural information on heterogeneous objects in a nearly native state and a rather simple sample preparation [3], cryo-EM has become a powerful tool in modern structural biology [4, 5]. Cryo-EM is currently the only method capable of addressing such challenges as the structure determination of difficult-to-crystallize membrane proteins and searching for different states
of molecular complexes, the method being complementary to classical X-ray crystallography [6].
In this work, the three-dimensional structure of cytochrome c nitrite reductase from the bacterium Thioalkalivibrio nitratireducens (TvNiR) [7] was studied for the first time by single-particle analysis [8] - a cryo-EM method for protein structure determination. This enzyme catalyzes the reduction of nitrite to ammonia without the release of reaction intermediates from the active site [7]. In solution and in the crystalline state, the enzyme exists as a stable hexamer that possess a bipyramidal shape with a characteristic height of ~150 A and a base of ~120 A [9]. The hexamer boasts D3 symmetry and a molecular weight of ~360 kDa, due to
Defocus estimation
Distribution of mean defocus value
C
<D
m m £
<u -Q £
s -i
S -
8 -
DefocusV, pm Astigmatism distribution
>80 nm ~4% of data
M 40 60 SO 100 120 140 160 160 200
Astigmatism, nm
<D
g a
<U _Q
<-1.5 pm ~1% of data
D
<u
g a
<u
_Q
i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—t—i—i—i—r—1
0.0 -0.2 -0.4 -0.6 -0.6 "1.0 -1.2 -1.4 -16 -1.6 -2.0
Defocus, pm Distribution of resolution estimation
>3.8 Â ~3.5% of data
I-1--1-1-1-r
11 13
-n
15
Resolution, Ä
B
Fig. 1. Principal characteristics of the initial data based on CTF parameters estimated with the Gctf program. A - the defocus distribution along two orthogonal axes; the color corresponds to the density of values at a specified coordinate. B—D - the distributions of the average defocus, astigmatism, and resolution assessment, respectively, along with the threshold values. The data with parameters in the ranges indicated in red were excluded from further processing
which TvNiR is a very convenient object for cryo-EM. The goal of this work was to determine the structure of TvNiR by cryo-EM and compare it with the high-resolution X-ray structure determined earlier [9-11].
EXPERIMENTAL
Protein isolation and purification
Native TvNiR was isolated and purified in two steps, according to a procedure described earlier [7]. Ani-on-exchange chromatography was performed on a 35-mL column prepacked with DEAE Sepharose Fast Flow at 4°C using a BioLogic LP system (BioRad, USA). The column was pre-equilibrated with 25 mM potassium phosphate buffer, pH 7.0. After loading of the extract and washing of the column with the same buffer, the protein was eluted with a linear gradient of 0-1.0 M NaCl. Gel filtration chromatography was performed
on an AKTA FPLC system (Amersham Biosciences, USA) equipped with a SuperdexTM200 10/300 column equilibrated with 50 mM potassium phosphate buffer, pH 7.0, supplemented with 0.15 M NaCl. For further structural studies, the protein was concentrated to 10 mg/mL.
Cryo-EM sample preparation
In order to determine the optimal protein concentration in solution, concentrations in a range of 0.1 to 6.0 mg/mL were tested. The protein solutions were applied to Lacey Carbon 300 mesh copper grids (Ted Pella, USA). The experimental data were collected using Quantifoil R1.2/1.3 300 mesh grids (Quantifoil, Germany) coated with a carbon support film containing regular arrays of 1.2-^m circular holes spaced by 2.5 ^m. The grids were glow-discharged for 30 s using a PELCO easyGlow system designed for hydrophilization
Summary of data collection parameters, cryo-EM map reconstruction, and structure refinement statistics for TvNiR
Data collection
Accelerating voltage, kV 300
Magnification 75000x
Pixel size, Â 0.86
Exposure time, s 1.5
Number of image stacks 3055
Total dose per stack, e-/Â2 ~100
Number of images per stack 30
Defocus range, ^m [-1.5; -0.5]
Defocus step size, ^m 0.1
Cryo-EM map reconstruction
Final number of particle images 33891
Symmetry group D3
FSC0 5 (with masking / without masking) 2.86 / 3.19
FSC0143 (with masking / without masking) 2.56 / 2.82
Resolution (average), A 2.56
Structure refinement
FSC 0.8679
Rf (weighted, overall), % 28.70
Average B-factor, A2 80.08
Rmsd bond lengths, A2 0.018
Rmsd bond angles, deg 1.989
MolProbity score, % 2.55
All-atom clashscore, % 13.31
Good rotamers, % 91.68
Ramachandran outliers, % 0.19
EMDB code EMD-0020
(Ted Pella, USA) at a chamber pressure of 0.26 mbar and a current of 25 mA. A sample (3 ^L, 6 mg/mL) was applied to the grids, and the grids were plunge-frozen in liquid ethane using a Vitrobot Mark IV vitrification device (Thermo Fisher Scientific, USA); the following parameters were set: blot force, 6; drain time, 1 s; temperature of the climate chamber, 4 °C; and humidity, 98 ± 2%.
Cryo-electron microscopy
The grids were transferred at liquid nitrogen temperature to a Titan Krios cryo-electron microscope (Thermo Fisher Scientific, USA) equipped with a Schottky-type field emission electron gun (FEI XFEG, the Netherlands), a spherical aberration corrector (CEOS GmbH, Germany), and a CMOS-based Falcon II direct electron detector (Thermo Fisher Scientific, USA). A total of 3,055 image stacks were recorded in an automatic mode using the EPU package (version 1.9.1.16REL) (Thermo Fisher Scientific, USA) with a total exposure time of 1.5 s. The microscope was operated at an accelerating voltage of 300 kV and 75000Xmagnification
Initial data: 3055 stacks of images
Global and local drift correction (Motioncor2, GPU)
2.
3055 images
CTF estimation (Gctf, GPU); data sorting
3.
2851 images
Template generation Particle picking
(Gautomatch+ Relion, GPU) (Gautomatch, GPU)
4.
44
435 558 particles
2D classification (Relion, GPU)
4.1
Initial model generation (EMAN2.1)
5.
6.
7.
171 130 particles
Density map refinement without symmetry and mask (Relion, GPU)
171 130 particles
Density map refinement with symmetry and mask (Relion, GPU)
2.9Â
2.65Â
171 130 particles
3D classification and examination of possible structural heterogeneity (Relion, GPU)
[39058; 25212; 33891; 22196; 22324; 28449] particles
Density map refinement for each class (Relion, GPU) Local resolution analysis (ResMap]
9.
2.56 Â
Final map 2.56 Â
Refinement of atomic model (REFMAC5)
Fig. 2. Main data processing steps
corresponding to a pixel size of 0.86 A at the specimen level, with objective lens defocused between -1.5 ^m and -0.5 ^m, with a step of 0.1 ^m, using a total dose of ~100e-/A2 evenly distributed across the image stack. The main data collection parameters are summarized in the Table. The principal data characteristics are shown in Fig. 1A-D.
RESULTS
Cryo-EM map reconstruction
The images were processed in several consecutive steps, presented in Fig. 2, using the computing resources of the Federal Collective Usage Center Complex for Simulation and Data Processing for Mega-Science
Fig. 3. Preliminary data processing after correction for drift and estimation of the CTF parameters. A - an image with projections of the particles picked with Gautomatch; B - the Fourier spectrum of this image is shown in the upper-right section of the panel and the rotational average of the power spectrum is in the lower-right section of the panel; Thon ring fitting is presented in the left section. The circle, the radius of which in direct space corresponds to the assessment of the maximum resolution (3.32 A), is highlighted
Facilities at the National Research Centre Kurchatov Institute, equipped with Nvidia Tesla K80 graphics accelerators. The following image processing packages were employed: Motioncor2 (version 1.0.5) [12], Gctf (version 1.18) [13], Gautomatch (version 0.56) (K. Zhang, MRC Laboratory of Molecular Biology, Cambridge, UK, http://www.mrc-lmb.cam.ac.uk/kzhang/ Gautomatch), and Relion (version 2.1) [14]. All these packages were optimized for computations on graphics processors.
In the first step (Fig. 2), 3,055 initial image stacks were individually corrected for beam-induced motion using Motioncor2. The following two averaged and corrected image sets were obtained: (1) images filtered depending on the electron dose exposed to the sample (Dose Weighting) [15], which were used for classification and refinement processes; and (2) non-filtered images that were used to estimate contrast transfer function (CTF) parameters. In the second step (Fig. 2), the CTF parameters were estimated with the Gctf program. For each image, the information limit (resolution assessment), defocus, and astigmatism were estimated based on Thon ring fitting. The distribution's tail points below the thresholds shown in Fig. 1B-D were excluded from further processing.
Therefore, 2,851 selected images with defocus, astigmatism, and resolution parameters not higher than -1.5 ^m (in magnitude), 80 nm, and 3.8 A, respectively, were used in all subsequent steps. A typical image after correction for drift is shown in Fig. 3.
In the third step (Fig. 2), particles were picked with Gautomatch. Initially, the procedure was applied for a subset of images recorded with a high degree of de-focus. The 2D Gaussian function with a half-width corresponding to the characteristic size of an object was used as a template. The resulting set of particles was subjected to 2D classification in Relion. The classes containing projections of the object were utilized as templates to pick particles from the total data set. The resulting set contained 435,558 coordinates of possible particle positions.
The fourth step (Fig. 2) involved two sequential rounds of classification. In the first round, the particle images were divided into 40 classes. Then, the images that were combined into classes not containing projections of the object or those depicting artifacts, such as ice crystals, surface contamination, and carbon edges, were excluded from the data set. In the second round, the remaining images were clustered into 50 classes, followed by the exclusion of the particles belonging to
6 2D class population
- 171130
tn O -o images of particles
g O " o
o © - © * •
s- •
"1-1-1-1-1-■-r
0 10 20 30 40 50
Class number
Fig. 4. Fourth data processing step. A - images merged in each class after the second classification round. Particle images belonging to classes 1 to 44 (numbered in blue) were used for further analysis and density map construction. Particle images belonging to classes 45 to 50 (numbered in red), which led to overtraining of reconstruction algorithms, were excluded from further data processing; B - class population. The colors of the points correspond to the colors of the numbers in panel A
classes in which the boundaries of the object were not well-defined (Fig. 4 A, B).
After classification, 171,130 particle images containing structural information on the object were selected. A low-resolution initial model was built with the EMAN2 image processing package [16] by the Monte-Carlo method, taking into account the known symmetry of the object based on 44 projections obtained by averaging images present in the most populated classes, after the second classification round (Fig. 4B).
In the fifth step (Fig. 2), the low-resolution model was refined by reassigning the Euler angles and fitting projections of the object to the cryo-EM map in each sequential iteration of the EM algorithm implemented in Relion [17-19]. The refinement was carried out using the 3D auto-refine procedure without resort-
ing to a-priori data on the symmetry of the object after postprocessing and applying the binary mask that defined the boundary conditions for the calculation of cross-correlation coefficients between two independently refined maps [20]. The resolution of the final map was 2.9 A, as determined by the FSC = 0.143 criterion [21].
The symmetry-imposed refinement in the sixth step (Fig. 2) resulted in resolution improvement to 2.65 A. In the seventh step, the same particle set (Fig. 2) was subjected to 3D classification [22, 23] into six classes (Fig. 5) without angular or translational searches, using the same mask according to a procedure described in [24]. In the eighth step (Fig. 2), the maps were repeatedly refined for each class. This approach allowed us to select a data subset composed of 33,891 particle images (median defocus was -0.86 ^m, with the values varying from -1.48 to -0.18 ^m, as estimated with Gctf), which corresponded to the third 3D class and yielded a map with the best resolution of 2.56 A (Fig. 5). In order to increase the resolution, the steps 6-8 (Fig. 2) were carried out with imposition of D3 symmetry and using a binary mask (Fig. 6A) that was created by applying 5-pixel isotropic extension, smoothing the boundaries by 5 pixels, and using an isosurface threshold of 0.02.
Structure refinement with REFMAC5
The crystal structure of TvNiR at the best resolution of 1.4 A (RCSB ID code is 3FO3) served as a starting model for the refinement. Solvent molecules, except for those corresponding to experimental density peaks, were removed from the dimer located in the asymmetric unit of the crystal. The TvNiR hexamer used for the refinement was generated by application of the appropriate symmetry operations. The refinement was performed with REFMAC5 [25] implemented in the CCP-EM suite [26].
A map of the third 3D class (Fig. 5), which had the best resolution, was used as the experimental cryo-EM map. To prevent the model from overfitting, the following approach was employed [27]. The images, from which the entire map was produced, were randomly divided into two subsets, and these subsets were used for the calculation of two independent 'half maps' and cross-validation. The high resolution of the experimental data allowed us to perform the refinement with no restraints, except for restraints on deviations of the bond lengths from the average value (jelly-body refinement [25]). To more correctly estimate steric constraints, hydrogen atoms were included in the refinement in fixed positions. After 30 refinement cycles, the fit of the refined model to the experimental density was visually inspected with Coot [28]. The influence of map sharpening/blurring on the refinement was assessed
FJ^B «s, : n 5.01 Â < t " > 3.02 Â 3 r.r. » ■ *■ it vf 2.56 Â
4 * * "V* t 2.86 Â 5 < • * - e . ""> vf 3.10 Â 6 ft?, ; * h Vf 3.06 Â
3
a
o
Q-
Class population with an increase of iteration step
Class 1 39058 Class 2 25212 Class 3 33891 Class 4 22196 Class 5 22324 Class 6 28449
10 15
Iteration step
B
Fig. 5. Results of 3D classification. A - central sections of the classes. Images from 1 to 44 (blue) were used for further analysis and density map reconstruction. Images from 45 to 50 (red) were excluded from further data processing, B - changes in the class populations depending on the step of iteration. Dots are colored in accordance with a numbering color on the A panel
by performing a series of refinements with different degrees of blurring, from -150 to +150. The appropriate parameters for further refinement were selected based on the best Rf and FSC, which corresponded to a blur parameter of -80. In addition to the visual inspection, the quality of the refined model was validated with Molprobity [29]. The refined model was compared with the crystal structure of the protein using PDBe-FOLD [30]. The Ramachandran plot analysis showed that the residues Gly285 and His361 of all subunits were in disallowed regions; however, these residues in all subunits had a well-defined electron density.
DISCUSSION
The 3D classification revealed no significant structural heterogeneity in the sample at the resolution level achieved, but it made it possible to select a subset of particles containing information on high spatial frequencies and providing maximum resolution for the final cryo-EM map (Fig. 6 B,C).
The distribution of projections for the angle classification, which was used in the refinement of the map without imposition of symmetry and the map for the third 3D class, is shown in Fig. 6D. It can be seen that the object has no preferred orientations in the amorphous ice layer. The assessment of the average resolution by the FSC = 0.143 and FSC = 0.5 criteria based on the results of postprocessing is shown in Fig. 7.
High symmetry of the object significantly simplifies the cryo-EM reconstruction, thereby compensating for a small number of projections with a rather high in-
A
90°
60°
C
D
Fig. 6. Final density map with the best resolution of 2.56 Â. A - the mask used for reconstruction and resolution assessment is shown by a gray semi-transparent isosurface; B—C - density maps in different projections, individual subunits of the protein hexamer are highlighted in colors; D - distribution of projections for the angle classification used
B
Fig. 7. Resolution assessment. A - local resolution for the density map of the third 3D class analyzed with MonoRes [31]; B - Fourier shell correlation (FSC) curves for the density maps without imposition of symmetry (gray curve), with symmetry constraints for the complete set of particle images (black curve), and for a particle set belonging to the third class (green curve)
Leucine (Leu389)
Threinine (Thr71)
Alanine (Ala67)
Valine (Val394)
Glutamine (Gln298)
Cysteine (Cys379)
Glycine (Gly285)
Methionine Phenylalanine Tryptophan (Met501) (Phe490) (Trp444)
Proline (Pro307)
Isoleucine (Ile246)
Tyrosine (Tyr425)
Asparagine (Asn432)
Serine (Ser399)
Fig. 8. Quality of the experimental density for the third 3D class. A - for the side chains of selected residues of the refined model; B - for the active-site heme and its environment
Lysine (Lys90)
Arginine (Arg239)
Histidine (His349)
Aspartic acid Glutamic acid (Aps130) (Glu417)
A
formation limit for this experiment and a low signal-to-noise ratio due to the low molecular weight of the object. The results of the refinement in Relion 2.1 after 3D classification clearly show the influence of the data quality on the final resolution. When applying symmetry of the object, the structural model was constructed using a total of 34,000 particle images from the 3D class that had the highest resolution of 2.56 A (Fig. 5).
The fitting of the enzyme X-ray structure (RCSB ID code 3FO3) to the experimental 2.56 A resolution cryo-EM map and the subsequent refinement yielded a final structure with a Molprobity score of 2.55 and the following parameters: R = 28.70, FSC = 0.8679
01 I ' average
(Table). The high-quality map allowed us not only to trace the polypeptide chain, but also to identify the side chains of residues (Fig. 8A), including the unique covalent bond between the active-site tyrosine and cysteine (Fig. 8B) and the side chains of residues Arg52, Arg316, and Lys456 on the surface of the enzyme molecule, which were invisible in the crystal structure. The density observed in the active site of the enzyme was identified as phosphate based on the composition of the buffer used for crystallization (Fig. 8B). The binding of inorganic anions is characteristic of the active site of TvNiR, which is well-known from X-ray crystallography [9, 10].
The superposition of the refined cryo-EM structure on the initial crystal structure shows that their backbone structures are, in general, similar. The side-chain orientations are also similar, except for some residues on the surface of the enzyme, e.g., Asp40, Arg51, Glu337, Glu341, etc., which may be attributed to their relatively high flexibility. The RMSD of the Ca atoms of these structures are not larger than 0.36 A. Therefore, the cryo-EM structure of TvNiR is in good agreement with the crystal structure of the enzyme determined
earlier by X-ray crystallography. The cryo-EM map was deposited in the Electron Microscopy Data Bank (EMDB) under the accession code EMD-0020.
CONCLUSION
The structure of cytochrome c nitrite reductase from the bacterium T. nitratireducens was studied by cryo-EM single-particle analysis; a cryo-EM map with a 2.56 A resolution was obtained, and the appropriate three-dimensional model was constructed. The optimal algorithm was found for data collection and processing in order to achieve high resolution. A comparison of the three-dimensional TvNiR structures determined by X-ray crystallography (1.40 A) and cryo-EM (2.56 A) revealed no significant differences. At the resolution level achieved by cryo-EM, TvNiR does not exhibit structural heterogeneity. •
This work was carried out using the computing resources of the Federal Collective Usage Center "Complex for Simulation and Data Processing for Mega-Science Facilities" at the National Research Centre "Kurchatov Institute," http://ckp.nrcki.ru (financial support from the Ministry of Education and Science of the Russian Federation, the unique project identifier RFMEFI62117X0016) and the facilities of the Resource Center of Probe and Electron Microscopy at the Centre for Nano, Bio, Info, Cognitive, and Social
Sciences and Technologies (NBICS Centre) of the National Research Centre "Kurchatov Institute" and
was supported by the Russian Science Foundation (grant No. 18-41-06001 supporting the data collection and processing and grant No. 14-24-00172 supporting
the protein isolation and purification) and the Government Task No. 01201351365 (cultivation of the producer strain).
REFERENCES
1. Callaway E. // Nature. 2015. V. 525. P. 172-174.
2. Merk A., Bartesaghi A., Banerjee S., Falconieri V., Rao P., Davis M.I., Pragani R., Boxer M.B., Earl L.A., Milne J.L.S., et al. // Cell. 2016. V. 165. № 7. P. 1698-1707.
3. Iancu C. V., Tivol W.F., Schooler J.B., Dias D.P., Henderson G.P., Murphy G.E., Wright E.R., Li Z., Yu Z., Briegel A., et al. // Nat. Protoc. 2007. V. 1. № 6. P. 2813-2819.
4. Cheng Y. // Cell. 2015. V. 161. № 3. P. 450-457.
5. Gelfand A. // Biomedical Computation Review. 2016. V.1 . P. 13-21.
6. Bai X., McMullan G., Scheres S.H.W. // Trends Biochem. Sci. 2015. V. 40. № 1. P. 49-57.
7. Tikhonova T.V., Slutsky A., Antipov A.N., Boyko K.M., Pol-yakov K.M., Sorokin D.Y., Zvyagilskaya R.A., Popov V.O. // Biochim. Biophys. Acta Proteins Proteomics. 2006. V. 1764. № 4. P. 715-723.
8. Frank J., Zhu J., Penczek P., Li Y., Srivastava S., Verschoor A., Radermacher M., Grassucci R., Lata R.K., Agrawal R.K.
// Nature. 1995. V. 376. № 6539. P. 441-444.
9. Polyakov K.M., Boyko K.M., Tikhonova T.V., Slutsky A., Antipov A.N., Zvyagilskaya R.A., Popov A.N., Bourenkov G.P., Lamzin V.S., Popov V.O. // J. Mol. Biol. 2009. V. 389. № 5. P. 846-862.
10. Trofimov A.A., Polyakov K.M., Boyko K.M., Tikhonova T.V., Safonova T.N., Tikhonov A.V., Popov A.N., Popov V.O. // Acta Crystallogr. Sect. D Biol. Crystallogr. 2010. V. 66. № 10. P. 1043-1047.
11. Trofimov A.A., Polyakov K.M., Boiko K.M., Filimonen-kov A.A., Dorovatovskii P.V., Tikhonova T.V., Popov V.O., Koval'chuk M.V. // Crystallogr. Reports. 2010. V. 55. № 1. P. 61-67.
12. Zheng S., Palovcak E., Armache J.-P., Cheng Y., Agard D. // bioRxiv. 2016. P. 1-30.
13. Zhang K. // J. Struct. Biol. 2016. V. 193. № 1. P. 1-12.
14. Kimanius D., Forsberg B.O., Scheres S.H.W., Lindahl E. // Elife. 2016. V. 5.
15. Grant T., Grigorieff N. // Elife. 2015. V. 4. P. e06980.
16. Tang G., Peng L., Baldwin P.R., Mann D.S., Jiang W., Rees I., Ludtke S.J. // J. Struct. Biol. 2007. V. 157. № 1. P. 38-46.
17. Dempster A.P.A., Laird N.M.N., Rubin D.D.B. Maximum likelihood from incomplete data via the EM algorithm. // Journal of the royal statistical society. Series B (methodological). 1977. P. 1-38
18. Scheres S.H.W. // J. Struct. Biol. 2012. V. 180. № 3. P. 519-530.
19. Scheres S.H.W. // J. Mol. Biol. 2012. V. 415. № 2. P. 406418.
20. Chen S., McMullan G., Faruqi A.R., Murshudov G.N., Short J.M., Scheres S.H.W., Henderson R. // Ultramicrosco-py. 2013. V. 135. P. 24-35.
21. Rosenthal P.B., Henderson R. // J. Mol. Biol. 2003. V. 333. № 4. P. 721-745.
22. Scheres S.H.W. // Methods in Enzymology. 2016. V. 579. Academic Press. P. 125-157.
23. Ludtke S.J. // Methods in Enzymology. 2016. V. 579. Academic Press. P. 159-189.
24. Hirschi M., Herzik M.A., Wie J., Suo Y., Borschel W.F., Ren D., Lander G.C., Lee S.Y. // Nature. 2017. V. 550. № 7676. P. 411-414.
25. Murshudov G.N., Skubák P., Lebedev A.A., Pannu N.S., Steiner R.A., Nicholls R.A., Winn M.D., Long F., Vagin A.A. // Acta Crystallogr. Sect. D Biol. Crystallogr. 2011. V. 67. № 4. P. 355-367.
26. Wood C., Burnley T., Patwardhan A., Scheres S., Topf M., Roseman A., Winn M. // Acta Crystallogr. Sect. D Biol. Crystallogr. 2015. V. 71. P. 123-126.
27. Brown A., Long F., Nicholls R.A., Toots J., Emsley P., Murshudov G. // Acta Crystallogr. Sect. D Biol. Crystallogr. 2015. V. 71. P. 136-153.
2B. Emsley P., Lohkamp B., Scott W.G., Cowtan K. // Acta Crystallogr. Sect. D Biol. Crystallogr. 2010. V. 66. № 4. P. 4B6-501.
29. Chen V.B., Arendall W.B., Headd J. J., Keedy D.A., Im-mormino R.M., Kapral G.J., Murray L.W., Richardson J.S., Richardson D.C. // Acta Crystallogr. Sect. D Biol. Crystal-logr. 2010. V. 66. № 1. P. 12-21.
30. Krissinel E., Henrick K. // Acta Crystallogr. Sect. D. 2004. V. 60. № 12-1. P. 2256-226B.
31. Vilas J.L., Gómez-Blanco J., Conesa P., Melero R., Miguel de la Rosa-Trevín J., Otón J., Cuenca J., Marabini R., Carazo J.M., Vargas J., et al. // Structure. 201B. V. 26. № 2. P. 337-344. e4.
