Application of Symmetry Adapted Function Method for Three-Dimensional Reconstruction of Octahedral Biological Macromolecules

Abstract

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 Degp₂₄ 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) =0.1,0.5, 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 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 [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 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 [13], FREALIGN [14] using the direct Fourier inversion method, MRC [15] 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 [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 X-ray 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], Hepatitis B small surface antigen particles (HbsAg) [25] tetrahedral (small heat shock protein ACR1 [26]), dihedral (auxilin-bound 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 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 Degp₂₄ (3Cs0.pdb) [24] and the Red-Cell 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 Symmetry-Adapted 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: 𝑂_𝑙,𝜇(𝜃,𝜑)=^𝑙_{𝑚=−𝑙,𝑚=0(mod4)}𝐵^𝜇_𝑙,𝑚𝑌_𝑙,𝑚(𝜃,𝜑), (1) where the 𝑂_𝑙,𝜇(𝜃,𝜑) represents the OSAFs, the 𝑌_𝑙,𝑚(𝜃,𝜑) denotes the normalized spherical harmonics, and the 𝐵^𝜇_𝑙,𝑚 is the combination coefficient with 𝑚=0(mod4) required by the 𝑍 axis being along a 4-fold axis, 𝜇(𝑙) is multiplicity of a given order 𝑙 of an OSAF. All the 𝜇(𝑙) for 𝑙≤11 are listed in the Table 1, when 𝑙≥12, the values of the 𝜇(𝑙) can be obtained by 𝑙𝜇(𝑙)=𝜇[mod(𝑙12)]+int(12), (2) where int(𝑙/12) is the integer part of the (𝑙/12) and mod(𝑙/12) the reminder of 𝑙 divided by 12.

Table 1 

The multiplicity 𝜇(𝑙) of the OSAF with 𝑙≤11.

To calculate OSAFs, at first, one should calculate the normalized spherical harmonics by the formula 𝑌^𝑚_𝑙(𝜃,𝜑)={(2𝑙+1)(𝑙−|𝑚|)!4𝜋(𝑙+|𝑚|)!}𝑃^𝑚_𝑙(cos𝜃)exp(𝑖𝑚𝜑), (3)

where the 𝑃^𝑚_𝑙(cos𝜃) 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 [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] 𝐵^𝜇_{𝑙1+𝑙2,𝑀}=𝛾(𝑙₁,𝑙₂,𝑙₁+𝑙₂)×^𝑙_𝑚=−𝑙(2𝑙₁)!(2𝑙₂)!(𝑙₁+𝑙₂+𝑀)!(𝑙₁+𝑙₂−𝑀)![2(𝑙₁+𝑙₂)+1]!(𝑙₁+𝑚)!(𝑙₁−𝑚)!(𝑙₂+𝑚−𝑀)!(𝑙₂−𝑚+𝑀)!×𝐵^𝜇1_𝑙1,𝑚𝐵^𝜇2_{𝑙2,𝑀−𝑚}, (4)where 𝛾(𝑙₁,𝑙₂,𝑙₁+𝑙₂) is normalized constant.

According to (4), any high-order OSAF can be generated only by using three lower order seed OSAFs with 𝑙=4, 6, 9. The coefficients 𝐵^𝜇_𝑙,𝑚 of the seed OSAFs are listed in the Table 2.

Table 2 

The coefficients of the OSAF for 𝑙=4,6,9.

For example, 𝑂_4,1(√(𝜃,𝜑)=21⋅𝑌⁰₄(𝜃,𝜑)+15/2⋅(𝑌₄⁻⁴(𝜃,𝜑)+𝑌⁴₄(𝜃,𝜑)))6,𝑂_9,1(√(𝜃,𝜑)=51⋅(𝑌⁴₉(𝜃,𝜑)−𝑌₉⁻⁴(𝜃,𝜑))−√21⋅(𝑌⁸₉(𝜃,𝜑)−𝑌₉⁻⁸(𝜃,𝜑))).12 (5)

Due to the properties of the normalized spherical harmonic functions, the 𝑌_𝑙,𝑚(𝜃,𝜑) and 𝐵^𝜇_𝑙,𝑚 satisfy the following relationships: 𝐵^𝜇_𝑙,𝑚=(−1)^𝑙−𝑚𝐵^𝜇_𝑙,−𝑚,𝑌^∗_𝑙,𝑚(𝜃,𝜑)=(−1)^𝑚𝑌_𝑙,−𝑚(𝜃,𝜑). (6)

So that 𝐵^𝜇_𝑙,𝑚𝑌_𝑙,𝑚(𝜃,𝜑)+𝐵^𝜇_𝑙,−𝑚𝑌_𝑙,−𝑚(𝜃,𝜑),=𝐵^𝜇_𝑙,𝑚[𝑌_𝑙,𝑚(𝜃,𝜑)+(−1)^𝑙𝑌^∗_𝑙,𝑚(𝜃,𝜑)].When𝑙iseven,𝑚≠0→𝑂^𝑐_𝑙,𝜇(𝜃,𝜑)=2𝐵^𝜇_𝑙,𝑚𝑃^𝑚_𝑙(cos𝜃)cos𝑚𝜑,𝑚=0→𝑂^𝑐_𝑙,𝜇(𝜃,𝜑)=𝐵^𝜇_𝑙,𝑚𝑃^𝑚_𝑙(cos𝜃)cos𝑚𝜑.When𝑙isodd,𝑚≠0→𝑂^𝑠_𝑙,𝜇(𝜃,𝜑)=2𝐵^𝜇_𝑙,𝑚𝑃^𝑚_𝑙(cos𝜃)sin𝑚𝜑,𝑚=0→𝑂^𝑠_𝑙,𝜇(𝜃,𝜑)=0, (7)(8) 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 𝑙≤12 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 (0),(4),(8), and (12) mean 𝑚=0,4,8, and 12, respectively. For example, according to (8), we can gain the OSAFs of 𝑙=6,9: 𝑂^𝑐_6,1(𝜃,𝜑)=0.35355339059×𝑃⁰₆(cos𝜃)cos(0)−0.66143782777×(2.0×𝑃⁴₆𝑂(cos𝜃)cos(4𝜑)),_𝑠9,1(𝜃,𝜑)=0.59511903571×(2.0×𝑃⁴₉(cos𝜃)sin(4𝜑))−0.38188130791×(2.0×𝑃⁸₉(cos𝜃)sin(8𝜑)). (9)

Table 3 

The OSAFs of all 𝑙≤12.

Figures 1(a) and 1(b) show the density contour lines of the OSAFs looking down along a fourfold axis with 𝑙=40, 𝜇=1 and 2, respectively. We can get all the higher order of OSAFs up to 1000 via (4).

(a)

(b)

Figure 1 

Density contour maps of the OSAFs looking down along a fourfold axis. (a) 𝑙=40, 𝜇=1, (b) 𝑙=40, 𝜇=2.

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 𝜌(𝐫)=𝐹(𝐑)exp(−2𝜋𝑖𝐑⋅𝐫)𝑅²sinΘ𝑑𝑅𝑑Θ𝑑Φ, (10) where 𝐫 and 𝐑 denote the vectors in real and Fourier spaces, respectively.

In contrast to Crowther’s Fourier-Bessel method [33], we utilize OSAF to express 𝐹(𝑅)𝐹(𝑅)=𝐹(𝑅,Θ,Φ)=^∞𝑛_𝑙=0^𝑙_𝜇=1𝑓_𝑙,𝜇(𝑅)𝑂_𝑙,𝜇(Θ,Φ), (11) 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, 𝐹_𝑟(𝑅,Θ,Φ)=^∞𝑛_𝑙=0(mod)2^𝐿_𝜇=1𝑓_{𝑙even,𝜇}(𝑅)𝑂^𝑐_𝑙,𝜇(Θ,Φ),𝑖𝐹_𝑖(𝑅,Θ,Φ)=𝑖^∞𝑛_𝑙=1(mod)2^𝑙_𝜇=1𝑓_𝑙odd,𝜇(𝑅)𝑂^𝑠_𝑙,𝜇(Θ,Φ), (12) where 𝐹_𝑟(𝑅,Θ,Φ) and 𝐹_𝑖(𝑅,Θ,Φ) denote the real and imaginary parts of 𝐹(𝑅,Θ,Φ), respectively.

Substituting (12) into (10), one can finally obtain 𝜌(𝑟,𝜃,𝜑)=4𝜋[^∞_𝑙=0mod2(−𝑖)^𝑙𝑛𝐿_𝜇=1(^∞₀𝑓_{𝑙even,𝜇}(𝑅)𝑗_𝑙(2𝜋𝑅𝑟)𝑅²𝑑𝑅)×𝑂_𝑙,𝜇(𝜃,𝜑)+^∞_𝑙=1mod2(−𝑖)^𝑙𝑖×^𝑛𝐿_𝜇=1(^∞₀𝑓_𝑙odd,𝜇(𝑅)𝑗_𝑙(2𝜋𝑅𝑟)𝑅²𝑑𝑅)𝑂_𝑙,𝜇(𝜃,𝜑)], (13) where 𝑗_𝑙(2𝜋𝑅𝑟) 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 𝑓_{𝑙even,𝜇}(𝑅) and 𝑓_𝑙odd,𝜇(𝑅) 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₂₄ 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 𝑆𝑁=signalnoise, (14)

where signal is the average value of the signal, and noise the average value of the noise.

Finally, 3D structures of two models by making use of all the 2000 projections (𝑆/𝑁=0.1, 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₂₄ and the Red-cell L Ferritin proteins with 𝑆/𝑁=0.1,0.5,0.8 and without noise, respectively. Figures 4 and 5 show the two-dimensional section maps perpendicular to 𝑍 axis with 𝑧=zero viewing along the 3-fold axis of DegP₂₄ 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 𝑆/𝑁=0.1. 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 𝑆/𝑁=0.1,0.5,0.8, 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 FSC=0.5 criterion, the nominal resolution of reconstructed results of the DegP₂₄ 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 𝑆/𝑁=0.1,0.5, and 0.8, respectively. It can be seen from Figures 6 and 7 that the reconstructed structures with 𝑆/𝑁=0.8 have the approximately nominal resolutions as the SSMs and the attainable resolutions decrease with the increase of the added noises (𝑆/𝑁=0.5 and 𝑆/𝑁=0.1). Naturally, more particles are needed to achieve high resolution for very low 𝑆/𝑁=0.1. 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 𝑆/𝑁=0.1,0.5, 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]: RE=∑|𝜌_𝑜−𝜌_𝑟|∑|𝜌_𝑜|×100%, (15) 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 𝑆/𝑁=0.1. 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.

Table 4 

The Relative Errors with different 𝑆/𝑁 at low Fourier frequency (𝑅=1/14.4 Å^-1).

(a)

(b)

(c)

(d)

Figure 2 

Three representative projection images of the DegP24 model generated from PDB data. (a) Without noise, (b) with S/N = 0.8, (c) with S/N = 0.5, (d) with S/N = 0.1.

(a)

(b)

(c)

(d)

Figure 3 

Three representative projection images of the erythrocyte L Ferritin model generated from PDB data. (a) Without noise, (b) with 𝑆/𝑁=0.8, (c) with 𝑆/𝑁=0.5, (d) with 𝑆/𝑁=0.1.

(a)

(b)

(c)

(d)

Figure 4 

Maps of two-dimensional sections perpendicular to 𝑍 axis with 𝑍=0 of the DegP24 model and reconstruction results viewing along the 3-fold axis by the OSAF method: (a) the standard model, (b) the reconstructed results with 𝑆/𝑁=0.1, (c) 𝑆/𝑁=0.5, (d) 𝑆/𝑁=0.8.

(a)

(b)

(c)

(d)

Figure 5 

Maps of two-dimensional sections perpendicular to 𝑍 axis with 𝑍=0 of the erythrocyte L Ferritin model and reconstruction results viewing along the 3-fold axis by the DSAF method: (a) the standard model, (b) the reconstructed results with 𝑆/𝑁=0.1, (c) 𝑆/𝑁=0.5, (d) 𝑆/𝑁=0.8.

(a)

(b)

(c)

(d)

(e)

Figure 6 

3D structures reconstructed by the OSAF method viewing along the 4-fold axis of DegP24 protein with 𝑆/𝑁=0.1, 0.5, and 0.8. (a) 3D structures of SSM, (b) the result for 𝑆/𝑁=0.1, (c) 𝑆/𝑁=0.5, (d) 𝑆/𝑁=0.8, (e) the FSC curves for different noise levels.

(a)

(b)

(c)

(d)

(e)

Figure 7 

3D structure viewing along the 4-fold axis of erythrocyte L Ferritin model with 𝑆/𝑁=0.1, 0.5, and 0.8. (a) The standard 3D structures, (b) the result for 𝑆/𝑁=0.1, (c) 𝑆/𝑁=0.5, (d) 𝑆/𝑁=0.8, (e) the FSC curves for different noise levels.

Figure 8 

The relative error between the erythrocyte L Ferritin model and the results reconstructed by the OSAF method.

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.

Table 5 

The algorithm time of the OSAF method for reconstructing DegP₂₄ model at low and high resolution (S means second).

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 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: 𝜙_𝑧1=𝜙_𝑧+𝜋2,𝜃_𝑥=𝜃_𝑦,𝜓_𝑧1=𝜓_𝑧−𝜋2, (16)

where 𝜙_𝑧1,𝜃_𝑥,𝜓_𝑧1 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.

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