Research Article  Open Access
Application of Symmetry Adapted Function Method for ThreeDimensional Reconstruction of Octahedral Biological Macromolecules
Abstract
A method for threedimensional (3D) reconstruction of macromolecule assembles, that is, octahedral symmetrical adapted functions (OSAFs) method, was introduced in this paper and a series of formulations for reconstruction by OSAF method were derived. To verify the feasibility and advantages of the method, two octahedral symmetrical macromolecules, that is, heat shock protein and the Redcell L Ferritin, were utilized as examples to implement reconstruction by the OSAF method. The schedule for simulation was designed as follows: 2000 random orientated projections of single particles with predefined Euler angles and centers of origins were generated, then different levels of noises that is signaltonoise ratio (S/N) , and 0.8 were added. The structures reconstructed by the OSAF method were in good agreement with the standard models and the relative errors of the structures reconstructed by the OSAF method to standard structures were very little even for high level noise. The facts mentioned above account for that the OSAF method is feasible and efficient approach to reconstruct structures of macromolecules and have ability to suppress the influence of noise.
1. Introduction
The determination of threedimensional (3D) structures of macromolecular assemblies plays a key role in understanding their functions and properties. In the course of structure reconstruction of single particles during the last several decades, the Cryoelectron microscopy (referred to as “CryoEM”) has been successively used to solve 3D structures up to subnanometer resolution [1–6], even nearatomic resolution, such as the 3.8 Å resolution structures of the CPV [7] and rotavirus [8], the 4.2 Å GroEL structure [9], and the 4.5 Å Epsilon 15 bacteriophage structure [10]. Such large assemblies often are either too large or too heterogeneous to be able to study by the conventionalXray crystallography and nuclear magnetic resonance (NMR) [11, 12]. It is considered that the CryoEM is an indispensable approach for determining the 3D structures of macromolecular complexes. And many different software packages have been developed in the world wide laboratories for 3D reconstructions, such as EMAN [13], FREALIGN [14] using the direct Fourier inversion method, MRC [15] using the FourierBessel synthesis method and spherical symmetryadapted functions (SAFs) method [16, 17]. The SAF method was first realized to be a more efficient method indicated by Crowther in their pioneer paper [18] three decades ago. Provencher and Vogel implemented 3D reconstruction by the SAF method taking both simulated structures and biological objects as samples [19, 20]. Zheng et al. used an icosahedral SAF (ISAF) method to determine the structures of viruses from solution Xray scattering data [21]. Navaza [16] systematically developed formulations for 3D reconstruction of icosahedral viruses including ab initio determination of origins and orientations of particles and interpolation of data in the reciprocal space by the ISAF method. Recently our group [17] reconstructed icosahedral symmetry biological objects (HBV, etc.) by icosahedral SAF approach which showed that SAF method is an efficient approach to reduce the influence of noise and achieve high resolution because of its ability of completely utilizing symmetry operation of the object being studied. Up to date, all the SAF method is only used in the reconstruction of icosahedral symmetry macromolecular. Due to a variety of symmetry of the macromolecular existing in nature except icosahedral symmetry, such as octahedral (small heat shock protein hsp16.5) [22, 23], heat shock protein Degp_{24} [24], Hepatitis B small surface antigen particles (HbsAg) [25] tetrahedral (small heat shock protein ACR1 [26]), dihedral (auxilinbound clathrin coat [27] catalase, ribulose bisphosphate, glutamine synthetase, carboxylase/oxygenase) and so forth, we would like to extend SAF method to reconstruct the objects having any other symmetry, for example, octahedral, tetrahedral, dihedral symmetries.
In this paper, we would like to concentrate our attention to the octahedral symmetryadapted function (OSAF) method. And a series of formulations for 3D reconstruction of octahedral symmetry macromolecules had been derived. To verify the feasibility and the advantages of this approach, two octahedral symmetrical simulated data such as heat shock protein Degp_{24} (3Cs0.pdb) [24] and the RedCell L Ferritin (1rcc.pdb) [28] downloaded from protein data bank (PDB) have been reconstructed by the OSAF method at high resolution. The results demonstrate that the OSAF method can retrieve the 3D structures of the octahedral symmetrical objects at high resolution and effectively suppress the influence of noise.
2. Principle of the OSAF Method
For readily understanding the OSAF method for 3D reconstruction, it is necessary to describe the principle of the OSAF method briefly.
2.1. Octahedral SymmetryAdapted Functions (OSAFs)
Due to SAFs being the linear combinations of the spherical harmonics [29], the major problem is how to find the coefficients of the OSAF (see expression (1)). According to the conventional definition, we choose the axes of the Cartesian coordinates system along three 4 fold axes of an octahedral symmetry group, respectively, and the relationship between the Cartesian and spherical pole coordinates is as convention. Consequently, the OSAF can be written as follows: where the represents the OSAFs, the denotes the normalized spherical harmonics, and the is the combination coefficient with required by the axis being along a 4fold axis, is multiplicity of a given order of an OSAF. All the for are listed in the Table 1, when , the values of the can be obtained by where is the integer part of the and the reminder of divided by 12.

To calculate OSAFs, at first, one should calculate the normalized spherical harmonics by the formula
where the is the Legendreassociated polynomial. Then, it is essential to solve the combination coefficients in (1) for getting the OSAFs. Up to now, there are many methods to obtain the highorder OSAFs, such as the general methods [30], Algebraic method [31], and the recursive approach [32]. But we prefer the recursive method since it is less sensitive to the computation errors and more stable for achieving the higher order OSAF as pointed out by Schmidt [32] where is normalized constant.
According to (4), any highorder OSAF can be generated only by using three lower order seed OSAFs with , 6, 9. The coefficients of the seed OSAFs are listed in the Table 2.

For example,
Due to the properties of the normalized spherical harmonic functions, the and satisfy the following relationships:
So that where the asterisk* denotes complex conjugation, and denote the real and the imaginary parts of the OSAF, respectively.
According to (8) one may use the to fit the real part of structure factor and the imaginary part of structure factor of biological objects with octahedral symmetry in reciprocal space. All the and of OSAFs with are listed in Table 3. The first column gives the orders of OSAF, the second gives the multiplicity of OSAF, and in the third denotes the and , respectively, the final column presents the combination coefficients of the OSAF, and the , and mean , and 12, respectively. For example, according to (8), we can gain the OSAFs of :

Figures 1(a) and 1(b) show the density contour lines of the OSAFs looking down along a fourfold axis with , and 2, respectively. We can get all the higher order of OSAFs up to 1000 via (4).
(a)
(b)
It should be pointed out that the OSAFs () are the orthogonal complete basis if the multiplicity is taken all the values, then any function with octahedral symmetry can be represented in terms of the linear combination of .
2.2. Reconstruction Principle
It is well known that the structure of macromolecular complexes can be described as its potential functions which are determined by Fourier inversion transformation of the structure factors , and its expression in spherical coordinates is where and denote the vectors in real and Fourier spaces, respectively.
In contrast to Crowther’s FourierBessel method [33], we utilize OSAF to express where is the fitting coefficient of the OSAF, and its value depends on the Fourier radius of a spherical shell and ; is the maximum multiplicity of a given order .
According to (8), (11) can be further expressed as the following two parts, that is, real and imaginary parts, where and denote the real and imaginary parts of , respectively.
Substituting (12) into (10), one can finally obtain where labels the spherical Bessel functions, and its recurrence relationship can be seen in reference [34].
The reconstruction by OSAF can be carried out in the following procedure.
(1)Calculate the OSAFs by (1) up to the required order.(2)Construct two linear equation groups with experimental determined structure factors and in the reciprocal space according to (12).(3)Find the fitting coefficients and by solving the above two linear equation groups by means of the least square method.(4)Determine 3D structures of octahedral symmetrical objects according to (13).3. Implementation of Reconstruction by the OSAF Method
To verify the feasibility and advantages of the OSAF method for reconstruction of macromolecules with octahedral symmetry, two biological objects with octahedral symmetry, heat shock protein Degp_{24} and the Redcell L Ferritin, were taken as examples. The atomic structures were downloaded from PDB (3cs0.pdb and 1rcc.pdb). The both 3D structures with 4.0 Å resolution were generated as standard structure models (SSMs) by the EMAN’s pdb2mrc procedure. Then two thousand random projections of these two proteins with predefined orientations and centers were created using realspace projection. Then random noise was added to each projection at 3 different signal to noise ratios (), that is, 0.1, 0.5., and 0.8 for 3D reconstruction according to the definition of which is described as below
where is the average value of the signal, and the average value of the noise.
Finally, 3D structures of two models by making use of all the 2000 projections (, 0.5, and 0.8) with predefined random Euler angles and center parameters were reconstructed at high resolutions by the OSAF approach. Figures 2 and 3 show typical projections of the DegP_{24} and the Redcell L Ferritin proteins with and without noise, respectively. Figures 4 and 5 show the twodimensional section maps perpendicular to axis with viewing along the 3fold axis of DegP_{24} and 4fold axis of Redcell L Ferritin proteins, respectively. From the twodimensional section maps, one may find that some noises still present in the final structures. Usually a perfect reconstruction is impossible in the case of added noise. However it can be seen intuitively from Figures 4 and 5 that the structures reconstructed by the OSAF method are in good agreement with the SSMs even with heavy noise . Figures 6 and 7 show the comparison of 3D reconstructed structures of above two proteins in high resolution, respectively, by the OSAF method. Although there is no obvious discrepancy between the 3D density map reconstructed by the OSAF method for , a few differences still can be identified. For quantitative comparison, the Fourier shell correlations (FSCs) [35] listed in Figures 6(e) and 7(e) were calculated, respectively. According to the criterion, the nominal resolution of reconstructed results of the DegP_{24} protein is approximately equal to 5.8 Å, 4.6 Å, and 4.3 Å and those of Redcell L Ferritin protein, are approximately equal to 5.6 Å, 4.7 Å, and 4.2 Å with , and 0.8, respectively. It can be seen from Figures 6 and 7 that the reconstructed structures with have the approximately nominal resolutions as the SSMs and the attainable resolutions decrease with the increase of the added noises ( and ). Naturally, more particles are needed to achieve high resolution for very low . To show the quality of the reconstructed structures by the OSAF method for different s in the real space quantitatively, the relative errors (REs) of the reconstructed 3D structures of Redcell L Ferritin model with added different noise levels such as , and 0.8 deviated to the SSMs have been calculated and presented in Figure 8. The formula for calculating RE is described as follows [36]: where RE denotes the Relative Errors between the SSMs and reconstructed structures, represents potential of SSM and is that of the normalized structure reconstructed by OSAF method. From Figure 8, one may find that with the increase of the Fourier frequency, the RE increases gradually. Furthermore it is apparent that in the case of the low Fourier frequency with nominal resolution of 14.4 Å, the relative errors listed in Table 4 show that REs keep almost constant for the different s, which means that the reconstructed structures are hardly influenced by added noise even for . As the structure reconstructed by the OSAF method is of very low RE, that is to say, the structure reconstructed by the OSAF method is very close to the real structure. As the added noises increase, the RE increases. The fact mentioned above implies that the OSAF method is feasible and efficient approach to reconstruct structures of macromolecules and can suppress the influence of noise since the OSAF method can completely utilize 24 symmetry operations of the octahedral symmetry. To achieve advantage of symmetry operation, one should use symmetry adapted function (SAF), for example, icosahedral, octahedral, tetrahedral, dihedral SAF, all these functions have ability to suppress influence of noise in different extent depending on the number of symmetry operation. The icosahedron have 60 symmetry operations which have the strongest ability which is verified by our former paper [17], but other SAFs have certain ability to suppress the influence of noise.

(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(e)
(a)
(b)
(c)
(d)
(e)
Since the of raw CryoEM data is very low, one may need a large number of particles to reconstruct a 3D structure to achieve high resolution. Therefore the program for dealing with these huge particles for 3D reconstruction is very time consumption. It is essential to reduce the computation time of the 3D reconstruction from huge particles at high resolution. To achieve fast computation, we managed to carry out the reconstruction in an asymmetrical unit of an octahedral symmetry, and therefore the calculation was speeded up 24 times so that the reconstruction can be performed with a PC computer which will be described in another paper in detail. Table 5 shows the algorithm time by the OSAF method. All the above tests were carried out at a general PC flat with the Pentium D 3.2 GHz CPU and the 2 G RAM. From the Table 5 one can see that the OSAF method is very fast even for high resolution reconstruction.

4. Conclusions
(1)A set of formulations for 3D reconstruction of macromolecular assemblies with octahedral symmetry by the OSAF method has been established.(2)The OSAF method is feasible and efficiently suppresses the influence of the noise because of its sufficiently utilizing the symmetry of the objects.(3)The calculation can be greatly speed up by dealing with the reconstruction in an asymmetrical unit of the octahedral symmetry group.It should be pointed out that in the simulation, one may use projections with predetermined centers and orientations to reconstruct structures; however in practice, one should reconstruct based on experimented measured data with unknown centers and orientations. In this case, one should first determine the center and orientation of a projection. At this stage, we did not write a program to determine the center and orientation by the OSAF method itself yet, So far we should use the other program such as EMAN [13], FREALIGN [14], and so forth, to determine the orientation and center parameters of a particle, then adopt the OSAF method for reconstruction. The orientation definition in our OSAF method is the ZXZ convention by “clockwise rotation’’ which is identical to the EMAN’s program and is different from ZYZ convention such as FREALIGN’s. The relationship of orientation definition between OSAF method and ZYZ convention can be described as follows which is in the same way as EMAN:
where is the ZXZ convention adopted in our OSAF method and EMAN, is the ZYZ convention adopted in SPIDER, IMAGIC, MRC, and FREALIGN.
The OSAF method for reconstruction is just at beginning stage, there is a plenty of space for optimizing the program. We believe that this method has a prospective future. A reconstruction with experimental data is proceeding based on the principle mentioned above and will be reported later on.
Acknowledgments
This work is supported by the National Natural Scientific Foundation of China (no. 10874144), National Basic Research Program of China (Grant no. 2010CB912400) and the Scientific Foundation of Education Department of Hunan Province (no. 08C885).
References
 B. Böttcher, S. A. Wynne, and R. A. Crowther, “Determination of the fold of the core protein of hepatitis B virus by electron cryomicroscopy,” Nature, vol. 386, no. 6620, pp. 88–91, 1997. View at: Publisher Site  Google Scholar
 J. F. Conway, N. Cheng, A. Zlotnick, P. T. Wingfield, S. J. Stahl, and A. C. Steven, “Visualization of a 4helix bundle in the hepatitis B virus capsid by cryoelectron microscopy,” Nature, vol. 386, no. 6620, pp. 91–94, 1997. View at: Publisher Site  Google Scholar
 R. Matadeen, A. Patwardhan, B. Gowen et al., “The Escherichia coli large ribosomal subunit at 7.5?Å resolution,” Structure, vol. 7, no. 12, pp. 1575–1583, 1999. View at: Publisher Site  Google Scholar
 Z. H. Zhou, M. Dougherty, J. Jakana, J. He, F. J. Rixon, and W. Chiu, “Seeing the herpesvirus capsid at 8.5 Å,” Science, vol. 288, no. 5467, pp. 877–880, 2000. View at: Publisher Site  Google Scholar
 Y. Modis, B. L. Trus, and S. C. Harrison, “Atomic model of the papillomavirus capsid,” EMBO Journal, vol. 21, no. 18, pp. 4754–4762, 2002. View at: Publisher Site  Google Scholar
 S. J. Ludtke, D.H. Chen, J.L. Song, D. T. Chuang, and W. Chiu, “Seeing GroEL at 6 Å resolution by single particle electron cryomicroscopy,” Structure, vol. 12, no. 7, pp. 1129–1136, 2004. View at: Publisher Site  Google Scholar
 X. Yu, L. Jin, and Z. H. Zhou, “3.88 Å structure of cytoplasmic polyhedrosis virus by cryoelectron microscopy,” Nature, vol. 453, pp. 415–419, 2008. View at: Publisher Site  Google Scholar
 X. Zhang, E. Settembre, C. Xu et al., “Nearatomic resolution using electron cryomicroscopy and singleparticle reconstruction,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 6, pp. 1867–1872, 2008. View at: Publisher Site  Google Scholar
 S. J. Ludtke, M. L. Baker, D.H. Chen, J.L. Song, D. T. Chuang, and W. Chiu, “De novo backbone trace of GroEL from single particle electron cryomicroscopy,” Structure, vol. 16, no. 3, pp. 441–448, 2008. View at: Publisher Site  Google Scholar
 W. Jiang, M. L. Baker, J. Jakana, P. R. Weigele, J. King, and W. Chiu, “Backbone structure of the infectious $\epsilon $15 virus capsid revealed by electron cryomicroscopy,” Nature, vol. 451, no. 7182, pp. 1130–1134, 2008. View at: Publisher Site  Google Scholar
 M. M. Golas, B. Sander, C. L. Will, R. Lührmann, and H. Stark, “Molecular architecture of the multiprotein splicing factor SF3b,” Science, vol. 300, no. 5621, pp. 980–984, 2003. View at: Publisher Site  Google Scholar
 Z. H. Zhou, H. Zhang, J. Jakana, X.Y. Lu, and J.Q. Zhang, “Cytoplasmic polyhedrosis virus structure at 8 Å by electron cryomicroscopy: structural basis of capsid stability and mRNA processing regulation,” Structure, vol. 11, no. 6, pp. 651–663, 2003. View at: Publisher Site  Google Scholar
 S. J. Ludtke, P. R. Baldwin, and W. Chiu, “EMAN: semiautomated software for highresolution singleparticle reconstructions,” Journal of Structural Biology, vol. 128, no. 1, pp. 82–97, 1999. View at: Publisher Site  Google Scholar
 N. Grigorieff, “FREALIGN: highresolution refinement of single particle structures,” Journal of Structural Biology, vol. 157, no. 1, pp. 117–125, 2007. View at: Publisher Site  Google Scholar
 R. A. Crowther, R. Henderson, and J. M. Smith, “The MRC image processing programs,” Journal of Structural Biology, vol. 116, no. 1, pp. 9–16, 1996. View at: Publisher Site  Google Scholar
 J. Navaza, “On the threedimensional reconstruction of icosahedral particles,” Journal of Structural Biology, vol. 144, pp. 13–23, 2003. View at: Google Scholar
 H. Liu, L. Cheng, S. Zeng, C. Cai, Z. H. Zhou, and Q. Yang, “Symmetryadapted spherical harmonics method for highresolution 3D singleparticle reconstructions,” Journal of Structural Biology, vol. 161, no. 1, pp. 64–73, 2008. View at: Publisher Site  Google Scholar
 R. A. Crowther, “Procedures for threedimensional reconstruction of spherical viruses by Fourier synthesis from electron micrographs,” Philosophical Transactions of the Royal Society of London. Series B, vol. 261, no. 837, pp. 221–230, 1971. View at: Google Scholar
 S. W. Provencher and R. H. Vogel, “Threedimensional reconstruction from electron micrographs of disordered specimens. I. Method,” Ultramicroscopy, vol. 25, no. 3, pp. 209–221, 1988. View at: Google Scholar
 R. H. Vogel and S. W. Provencher, “Threedimensional reconstruction from electron micrographs of disordered specimens. II. Implementation and results,” Ultramicroscopy, vol. 25, no. 3, pp. 223–239, 1988. View at: Google Scholar
 Y. Zheng, P. C. Doerschuk, and J. E. Johnson, “Determination of threedimensional lowresolution viral structure from solution Xray scattering data,” Biophysical Journal, vol. 69, no. 2, pp. 619–639, 1995. View at: Google Scholar
 K. K. Kim, R. Kim, and S.H. Kim, “Crystal structure of a small heatshock protein,” Nature, vol. 394, pp. 595–599, 1998. View at: Publisher Site  Google Scholar
 D. A. Haley, M. P. Bova, Q.L. Huang, H. S. Mchaourab, and P. L. Stewart, “Small heatshock protein structures reveal a continuum from symmetric to variable assemblies,” Journal of Molecular Biology, vol. 298, no. 2, pp. 261–272, 2000. View at: Publisher Site  Google Scholar
 T. Krojer, J. Sawa, E. Schäfer, H. R. Saibil, M. Ehrmann, and T. Clausen, “Structural basis for the regulated protease and chaperone function of DegP,” Nature, vol. 453, pp. 885–890, 2008. View at: Publisher Site  Google Scholar
 R. J. C. Gilbert, L. Beales, D. Blond et al., “Hepatitis B small surface antigen particles are octahedral,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 41, pp. 14783–14788, 2005. View at: Publisher Site  Google Scholar
 C. K. Kennaway, J. L. P. Benesch, U. Gohlke et al., “Dodecameric structure of the small shock protein ACR1 from Mycobacterium tuberculosis,” The Journal of Biological Chemistry, vol. 280, pp. 33419–33425, 2005. View at: Publisher Site  Google Scholar
 A. Fotin, Y. Cheng, N. Grigorieff, T. Walz, S. C. Harrison, and T. Kirchhausen, “Structure of an auxilinbound clathrin coat and its implications for the mechanism of uncoating,” Nature, vol. 432, pp. 649–653, 2004. View at: Publisher Site  Google Scholar
 J. Trikha, E. C. Theil, and N. M. Allewell, “High resolution crystal structures of amphibian redcell L ferritin: potential roles for structural plasticity and solvation in function,” Journal of Molecular Biology, vol. 248, no. 5, pp. 949–967, 1995. View at: Publisher Site  Google Scholar
 C. J. Bradley and A. P. Cracknell, The Mathematical Theory of Symmetry in Solid, Clarendon, Oxford, UK, 1972.
 S. L. Altnann, “On the symmetries of spherical harmonics,” Proceedings of the Cambridge Philosophical Society, vol. 53, no. 2, pp. 343–367, 1957. View at: Google Scholar
 P.D. Fan and J.Q. Chen, “Algebraic solutions for the octahedral group: group chain $O\supset {C}_{4}$,” International Journal of Quantum Chemistry, vol. 74, no. 1, pp. 7–22, 1999. View at: Google Scholar
 B. Schmidt and P. Zdanska, “Solution of the timedependent Schrödinger equation for highly symmetric potentials,” Computer Physics Communications, vol. 127, pp. 290–308, 2000. View at: Publisher Site  Google Scholar
 R. A. Crowther, D. J. DeRosier, and A. Klug, “The reconstruction of a threedimensional structure from projections and its application to electron microscopy,” Proceedings of the Royal Society of London. Series A, vol. 317, no. 1530, pp. 319–340, 1970. View at: Publisher Site  Google Scholar
 M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover Press, New York, NY, USA, 1965.
 W. O. Saxton and W. Baumeister, “The correlation averaging of a regularly arranged bacterial cell envelope protein,” Journal of Microscopy, vol. 127, no. 2, pp. 127–138, 1982. View at: Google Scholar
 M. G. Rossmann, “Fitting atomic models into electronmicroscopy maps,” Acta Crystallographica Section D, vol. 56, no. 10, pp. 1341–1349, 2000. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2010 Songjun Zeng et al. 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.