Mathematical Methods for Images and SurfacesView this Special Issue
Research Article | Open Access
Songjun Zeng, Hongrong Liu, Qibin Yang, "Application of Symmetry Adapted Function Method for Three-Dimensional Reconstruction of Octahedral Biological Macromolecules", International Journal of Biomedical Imaging, vol. 2010, Article ID 195274, 11 pages, 2010. https://doi.org/10.1155/2010/195274
Application of Symmetry Adapted Function Method for Three-Dimensional Reconstruction of Octahedral Biological Macromolecules
A method for three-dimensional (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 Red-cell 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 signal-to-noise 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.
The determination of three-dimensional (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 Cryo-electron microscopy (referred to as “CryoEM”) has been successively used to solve 3D structures up to subnanometer resolution [1–6], even near-atomic resolution, such as the 3.8 Å resolution structures of the CPV  and rotavirus , the 4.2 Å GroEL structure , and the 4.5 Å Epsilon 15 bacteriophage structure . Such large assemblies often are either too large or too heterogeneous to be able to study by the conventionalX-ray 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 , FREALIGN  using the direct Fourier inversion method, MRC  using the Fourier-Bessel synthesis method and spherical symmetry-adapted functions (SAFs) method [16, 17]. The SAF method was first realized to be a more efficient method indicated by Crowther in their pioneer paper  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 X-ray scattering data . Navaza  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  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 Degp24 , Hepatitis B small surface antigen particles (HbsAg)  tetrahedral (small heat shock protein ACR1 ), dihedral (auxilin-bound clathrin coat  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 symmetry-adapted 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 Degp24 (3Cs0.pdb)  and the Red-Cell L Ferritin (1rcc.pdb)  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 Symmetry-Adapted Functions (OSAFs)
Due to SAFs being the linear combinations of the spherical harmonics , 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 4-fold 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 Legendre-associated 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 high-order OSAFs, such as the general methods , Algebraic method , and the recursive approach . 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  where is normalized constant.
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 :
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 Fourier-Bessel method , 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 .
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 Degp24 and the Red-cell 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 real-space 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 DegP24 and the Red-cell L Ferritin proteins with and without noise, respectively. Figures 4 and 5 show the two-dimensional section maps perpendicular to axis with viewing along the 3-fold axis of DegP24 and 4-fold axis of Red-cell L Ferritin proteins, respectively. From the two-dimensional 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)  listed in Figures 6(e) and 7(e) were calculated, respectively. According to the criterion, the nominal resolution of reconstructed results of the DegP24 protein is approximately equal to 5.8 Å, 4.6 Å, and 4.3 Å and those of Red-cell 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 Red-cell 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 : 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 , but other SAFs have certain ability to suppress the influence of noise.
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 , FREALIGN , 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 Z-X-Z convention by “clockwise rotation’’ which is identical to the EMAN’s program and is different from Z-Y-Z convention such as FREALIGN’s. The relationship of orientation definition between OSAF method and Z-Y-Z convention can be described as follows which is in the same way as EMAN:
where is the Z-X-Z convention adopted in our OSAF method and EMAN, is the Z-Y-Z 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.
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).
- 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.
- J. F. Conway, N. Cheng, A. Zlotnick, P. T. Wingfield, S. J. Stahl, and A. C. Steven, “Visualization of a 4-helix bundle in the hepatitis B virus capsid by cryo-electron microscopy,” Nature, vol. 386, no. 6620, pp. 91–94, 1997.
- 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.
- 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.
- 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.
- 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.
- X. Yu, L. Jin, and Z. H. Zhou, “3.88 Å structure of cytoplasmic polyhedrosis virus by cryo-electron microscopy,” Nature, vol. 453, pp. 415–419, 2008.
- X. Zhang, E. Settembre, C. Xu et al., “Near-atomic resolution using electron cryomicroscopy and single-particle reconstruction,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 6, pp. 1867–1872, 2008.
- 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.
- W. Jiang, M. L. Baker, J. Jakana, P. R. Weigele, J. King, and W. Chiu, “Backbone structure of the infectious 15 virus capsid revealed by electron cryomicroscopy,” Nature, vol. 451, no. 7182, pp. 1130–1134, 2008.
- 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.
- 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.
- S. J. Ludtke, P. R. Baldwin, and W. Chiu, “EMAN: semiautomated software for high-resolution single-particle reconstructions,” Journal of Structural Biology, vol. 128, no. 1, pp. 82–97, 1999.
- N. Grigorieff, “FREALIGN: high-resolution refinement of single particle structures,” Journal of Structural Biology, vol. 157, no. 1, pp. 117–125, 2007.
- 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.
- J. Navaza, “On the three-dimensional reconstruction of icosahedral particles,” Journal of Structural Biology, vol. 144, pp. 13–23, 2003.
- H. Liu, L. Cheng, S. Zeng, C. Cai, Z. H. Zhou, and Q. Yang, “Symmetry-adapted spherical harmonics method for high-resolution 3D single-particle reconstructions,” Journal of Structural Biology, vol. 161, no. 1, pp. 64–73, 2008.
- R. A. Crowther, “Procedures for three-dimensional 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.
- S. W. Provencher and R. H. Vogel, “Three-dimensional reconstruction from electron micrographs of disordered specimens. I. Method,” Ultramicroscopy, vol. 25, no. 3, pp. 209–221, 1988.
- R. H. Vogel and S. W. Provencher, “Three-dimensional reconstruction from electron micrographs of disordered specimens. II. Implementation and results,” Ultramicroscopy, vol. 25, no. 3, pp. 223–239, 1988.
- Y. Zheng, P. C. Doerschuk, and J. E. Johnson, “Determination of three-dimensional low-resolution viral structure from solution X-ray scattering data,” Biophysical Journal, vol. 69, no. 2, pp. 619–639, 1995.
- K. K. Kim, R. Kim, and S.-H. Kim, “Crystal structure of a small heat-shock protein,” Nature, vol. 394, pp. 595–599, 1998.
- D. A. Haley, M. P. Bova, Q.-L. Huang, H. S. Mchaourab, and P. L. Stewart, “Small heat-shock protein structures reveal a continuum from symmetric to variable assemblies,” Journal of Molecular Biology, vol. 298, no. 2, pp. 261–272, 2000.
- 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.
- 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.
- 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.
- A. Fotin, Y. Cheng, N. Grigorieff, T. Walz, S. C. Harrison, and T. Kirchhausen, “Structure of an auxilin-bound clathrin coat and its implications for the mechanism of uncoating,” Nature, vol. 432, pp. 649–653, 2004.
- J. Trikha, E. C. Theil, and N. M. Allewell, “High resolution crystal structures of amphibian red-cell L ferritin: potential roles for structural plasticity and solvation in function,” Journal of Molecular Biology, vol. 248, no. 5, pp. 949–967, 1995.
- 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.
- P.-D. Fan and J.-Q. Chen, “Algebraic solutions for the octahedral group: group chain ,” International Journal of Quantum Chemistry, vol. 74, no. 1, pp. 7–22, 1999.
- B. Schmidt and P. Zdanska, “Solution of the time-dependent Schrödinger equation for highly symmetric potentials,” Computer Physics Communications, vol. 127, pp. 290–308, 2000.
- R. A. Crowther, D. J. DeRosier, and A. Klug, “The reconstruction of a three-dimensional 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.
- 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.
- M. G. Rossmann, “Fitting atomic models into electron-microscopy maps,” Acta Crystallographica Section D, vol. 56, no. 10, pp. 1341–1349, 2000.
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.