research papers
Ionic scattering factors of atoms that compose biological molecules
aBiostructural Mechanism Laboratory, RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan, and bComputational Life Science Laboratory, Graduate School of Medical Life Science, Yokohama City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
*Correspondence e-mail: yone@spring8.or.jp
Ionic scattering factors of atoms that compose biological molecules have been computed by the multi-configuration Dirac–Fock method. These ions are chemically unstable and their scattering factors had not been reported except for O−. Yet these factors are required for the estimation of partial charges in protein molecules and The electron scattering factors of these ions are particularly important as the electron scattering curves vary considerably between neutral and charged atoms in the spatial-resolution range explored in structural biology. The calculated X-ray and electron scattering factors have then been parameterized for the major scattering curve models used in X-ray and electron protein crystallography and single-particle cryo-EM. The X-ray and electron scattering factors and the fitting parameters are presented for future reference.
Keywords: form factors; electron crystallography; single-particle cryo-EM; structure refinement; X-ray crystallography; imaging; structure determination.
1. Introduction
Scattering factors for X-rays and electrons are essential for the analysis of experimental data probed with an X-ray and an electron beam. X-rays are scattered by electrons around the atom, yielding an electron-density map, whereas electrons are scattered by Coulomb potential, yielding a Coulomb-potential map. Hence the scattering factors are directly related to the charged state, particularly the electron scattering factors, which vary considerably between the neutral and charged atoms in the regime when sinθ/λ is smaller than 0.1 Å−1 and do not match until sinθ/λ reaches 0.2–0.3 Å−1 (Fig. 1). Here θ represents half the scattering angle, λ is the wavelength of the incident X-ray or electrons and the latter values of sinθ/λ correspond to spatial resolutions of 1/(0.2 × 2) – 1/(0.3 × 2) = 2.5 – 1.67 Å. Thus, the charge state essentially affects the structure analysis of proteins and protein complexes including by electron crystallography and single-particle cryo-EM.
The charge is actually delocalized over several atoms in the protein molecules and for theoretical partial charge distributions in amino acids and a nucleic acid). Although this approach is more qualitative than using form factors representing chemical bonding calculated for typical small molecules (Chang et al., 1999; Yamashita & Kidera, 2001; Zhong et al., 2002), it allows for the information related to charges to be extracted by electron two-dimensional (Mitsuoka et al., 1999) and three-dimensional (3D) crystallography (Yonekura et al., 2015; Yonekura & Maki-Yonekura, 2016) and single-particle cryo-EM (Yonekura & Maki-Yonekura, 2016). Also, ionic scattering factors could lead to a more accurate structure against cryo-EM data (Yonekura & Maki-Yonekura, 2016). However, those studies assumed partial charges for a limited number of atoms including O and H of the side chains in titratable residues such as aspartate, glutamate, lysine, arginine and histidine.
yielding partial charges which can be regarded as effective scattering factors. Partial charges should be assigned to all of the atoms in proteins and protein complexes (see Fig. 2The scattering factors of atoms and ions for X-rays and electrons were calculated (e.g. Doyle & Turner, 1968; Schmidt & Weiss, 1979; Rez et al., 1994; Wang et al., 1996; Su & Coppens, 1997; Macchi & Coppens, 2001) and parameterized previously (e.g. Doyle & Turner, 1968; Rez et al., 1994; Waasmaier & Kirfel, 1995; Peng et al., 1996; Su & Coppens, 1997; Peng, 1998; Macchi & Coppens, 2001; Colliex et al., 2006), but there have been no reported scattering factors for chemically unstable ions involved in proteins and protein complexes except for O−. To fill the gap, we calculated the scattering factors of those ions by the multi-configuration Dirac–Fock (MCDF) method in this study. We only selected ions for the atoms, C, N, O, P and S that compose amino acids and which are the major building blocks of biological molecules. The electron scattering factor of H+ is also presented. We then provide coefficients for the purpose of fitting the calculated scattering factors to the major curve models used in X-ray and electron crystallography and single-particle cryo-EM.
2. Calculations
2.1. Calculation of radial wavefunctions
Radial wavefunctions in the ground state were calculated for C, C+, C−, N, N+, N−, O, O+, O−, O2−, P, P+, P−, P2+, P3+, S, S+, S−, S2+, S2−, S3+ and S4+ by the MCDF method. The GRASP2K package (version 1.1, Jönsson et al., 2013) was used to perform the extended optimal level calculation for self-consistent fields of wavefunctions. Isolated negative ions are known to be unstable with a short lifetime in vacuum and tend to suffer from the inherent lack of convergence by MCDF calculations. A practical approach to this problem is placing positive charges outside a given radius around the ion (Watson, 1958; Schmidt & Weiss, 1979; Rez et al., 1994), but the resultant scattering curves differ significantly in the range sinθ/λ < 0.4 Å−1, depending on the size of the radius (Schmidt & Weiss, 1979). Previous studies (Wang et al., 1996; Macchi & Coppens, 2001) showed that the Dirac–Fock calculation yielded good solutions for O− and halides with no surrounding positive charges. Based on these studies, we carried out MCDF calculations in the absence of surrounding positive charges; this treatment was straightforward for O−, P− and S− as well as for the neutral atoms and positive ions with the exception of P3+ and S4+. For C− and N−, starting from initial estimates of radial functions for the neutral atom and/or increasing array sizes in the GRASP2K programs yielded reasonable solutions. No stable solutions were obtained using this approach for the closed-shell ions: O2−, S2−, P3+ and S4+.
2.2. Conversion to scattering factors
The scattering factors were calculated from relativistic wavefunctions as described in Su & Coppens (1997). Briefly, the radial charge density ρ(r) at a radius of r can be calculated as
where Pi(r) and Qi(r) represent the major and minor components, respectively, of the radial wavefunction for the ith orbital, and Nqi is the generalized occupation number for the corresponding orbital. The X-ray scattering factor Fx(s) at s = sinθ/λ can then be converted from ρ(r) by Fourier transform in polar coordinates as
The integral in equation (2) was obtained by the composite-Simpson summation at a fine step of r given by cubic spline interpolation. The scattering factors were computed up to s = 12 Å−1. A python script called scsumrhofft.py (provided in the supporting information) was written for these calculations in equations (1) and (2).
2.3. Parameterization of scattering factors
The ScatCurve package (Yonekura & Maki-Yonekura, 2016) was used to parameterize the X-ray scattering factors for the major curve models used in X-ray protein crystallography expressed in the form
where n = 4 for the four Gaussians plus a constant model (Colliex et al., 2006) and n = 5 for the five Gaussians plus a constant model (Waasmaier & Kirfel, 1995). The program scatcurvefit was modified to achieve a more robust minimization of the difference between the scattering factor and the model.
The X-ray scattering factor Fx was converted to the electron scattering factor Fel for an atom with Z0 electrons and Z nuclear charges using the Mott formula:
where m0 represents the electron mass, e is the 0 is the permittivity of free space, h is Planck's constant, and (m0e2/8π∊0h2) = 0.02393366. ΔZ is the ionic charge and is defined as ΔZ = Z – Z0. The scattering factor of H+ was also obtained from equation (4) with Fx of H+ = 0 and ΔZ = 1. By using scatcurvefit (Yonekura & Maki-Yonekura, 2016) again, the electron scattering factors were parameterized for the five Gaussians plus a charge term model (Peng, 1998) expressed as
The calculated scattering factors and fitting to the curve model were evaluated from an Rscat factor (Peng, 1998; Yonekura & Maki-Yonekura, 2016) by using scatcurvediff and scatcurvefit. Rscat is defined as
where F is the scattering amplitude computed in this work and Fc is the reference value or fitted value for the curve model.
3. Results and discussion
3.1. Validation of scattering factors
For validation of the X-ray scattering factors obtained from the MCDF calculations, we compared the calculated values for C, N, O, P, S and O− with the reference values, which are tabulated in the International Tables for Crystallography (Colliex et al., 2006). The tables contain these data up to s = 6 Å−1 for neutral atoms and 1.5 Å−1 for O−. Rscat values defined in equation (6) were calculated in these ranges (Table 1). The slight differences of 0.02–0.01% indicate that the values calculated in this study are consistent with the reference data. Plots of the calculated scattering factors and the reference curves appeared almost identical.
|
Relativistic effects are known to be small for light atoms (e.g. Wang et al., 1996). Despite this, we included relativistic effects as most of the scattering factors in Colliex et al. (2006) adopted the relativistic calculation, even for the atoms concerned in this study. The MCDF calculation for those atoms does not require much computing cost on the current PC system. The reference scattering factor for O− was non-relativistically computed (Colliex et al., 2006) and the difference from our relativistic calculation is very small (Table 1), confirming a very small contribution from the relativistic term.
The scattering factors of the same atom with and without a charge match well when s = 0.2–0.3 Å−1 (Fig. 1). The curves of the calculated X-ray and electron scattering factors show these expected appearances, see Figs. 3 and 4, respectively. Figs. 3 and 4 also indicate that the point, where the curves of different charge states become overlapped, varies for each atom. In contrast, the electron scattering factors of H given in the International Tables for Crystallography (Colliex et al., 2006) and of H+ derived in Hirai et al. (2007) do not match even beyond s = 0.5 Å−1 (see Fig. 1 of Yonekura & Maki-Yonekura, 2016). We calculated the electron scattering factor of H+ from equation (4) by simply setting Fx of H+ to 0 and ΔZ = 1. The new curve matches that of neutral H as sinθ/λ increases (Fig. 4). The calculated X-ray and electron scattering factors, including the new values for H+, are tabulated up to s = 6 Å−1 in Tables S1 and S2, respectively (see supporting information).
We also tried the standard quantum chemistry programs Gaussian 09 (Frisch et al., 2016) and GAMESS (Schmidt et al., 1993) for the calculation of the scattering factors, with the hope of obtaining the scattering factors of O2−, S2−, P3+ and S4+, as the MCDF calculation gave no stable solution for these ions. calculations with ROHF, UHF, MP2,3,4 and CCSD(T) were performed using the various Gaussian-type basis sets such as cc-pVTZ, cc-pV5Z, aug-cc-pVTZ, d-aug-cc-pVTZ and so on. Charge density ρ was then calculated on a sufficiently fine grid and converted to the radial scattering factor. The Rscat values, however, were around 0.5–1.0% against the references for neutral atoms and O−, and were much larger than the errors obtained from the MCDF calculations (Table 1), although the basis sets above were supposed to be available for quantitative calculation. Plots of the scattering factors appear to vary from the reference, particularly in the range of low sinθ/λ. This inconsistency may reflect that the Gaussian-type basis set cannot describe the cusps near the nucleus and inter-electron coalescent points (e.g. Pachucki & Komasa, 2004). Therefore, we did not adopt these calculations.
3.2. Parameterization of scattering factors
We then parameterized the X-ray scattering factors in Table S1 for the four Gaussians plus a constant (Table 2) and five Gaussians plus a constant models (Table 3), and the electron scattering factors in Table S2 for the five Gaussians plus a charge term model (Table 4). Data up to s = 6 Å−1 were used for all X-ray scattering factors and data from s = 0.04 to 6 Å−1 for the electron scattering factors of the positive ions. Only data from s = 0.04 to 1.5 Å−1 were used for the electron scattering factors of the negative ions, as Rscat values between the scattering factors and the fitted curves were getting worse when including data beyond s = 1.5 Å−1. Calculated coefficients and Rscat values are summarized in Tables 2, 3 and 4. All the Rscat values are less than 0.09% and typically around 0.05–0.005%. The electron scattering factors of the negative ions give slightly worse Rscat values, but these values are comparable with those reported for other ions (Peng, 1998). Equations (4) and (5) are identical for H+ when all a1–5 = 0, and do not need to be parameterized for the five Gaussians plus a charge term model. Scattering factors of multivalent cations of phosphorus and sulfur may be necessary for some chemicals such as nicotinamide adenine dinucleotide phosphates [see the crystallographic information files (Brown & McMahon, 2002) in the PHENIX library (Adams et al., 2010)]. Tables S1, S2 and 2–4 are provided in the supporting information for future reference.
|
|
|
3.3. Radiation damage
Density losses in the negatively charged side chains of aspartate and glutamate were often observed in recent single-particle reconstructions (e.g. Bartesaghi et al., 2014; Fromm et al., 2015; Hryc et al., 2017) and are likely to be interpreted as radiation damage (Fioravanti et al., 2007) even when reconstructed from data sets collected with a low accumulated electron dose. Part of the density loss should be attributed to negatively charged atoms, and this can be checked if a map calculated without lower-resolution data recovers the corresponding density (Yonekura et al., 2015). In contrast, an intense beam such as an X-ray-free electron laser causes ionization damage to samples and this may affect the diffraction patterns (Hau-Riege, 2007). The scattering factors of these radicals losing some electrons may help us to treat the data properly (Hau-Riege, 2007).
4. Conclusions
We have provided the X-ray and electron scattering factors of C+, C−, N+, N−, O+, O−, P+, P−, P2+, S+, S−, S2+ and S3+, the electron scattering factor of H+ and the coefficients for the four Gaussians plus a constant, the five Gaussians plus a constant and the five Gaussians plus a charge term models. The scattering factors of partially charged atoms and the fitting parameters for a curve model can be obtained by a linear combination of those for neutral and fully ionized atoms using scatcurvecomb (as in Yonekura & Maki-Yonekura, 2016). We are now testing the calculated values in this study for structure against data obtained by electron 3D crystallography and single-particle cryo-EM.
Supporting information
Supplementary Tables. DOI: https://doi.org/10.1107/S2052252518005237/fq5001sup1.pdf
Supplementary code. DOI: https://doi.org/10.1107/S2052252518005237/fq5001sup2.zip
Acknowledgements
We thank Toshiyuki Azuma for comments on features of divalent anions. Calculations with Gaussian 09 were performed in the HOKUSAI GreatWave HPC system, RIKEN.
Funding information
This work was supported by the Japan Science and Technology Agency SENTAN program (to KY), Japan Society for the Promotion of Science Grant-in-Aid for Challenging Exploratory Research Grant 24657111 (to KY), and Japan Society for the Promotion of Science Grant-in-Aid for Scientific Research Grant 16 H04757 (to KY).
References
Adams, P. D., Afonine, P. V., Bunkóczi, G., Chen, V. B., Davis, I. W., Echols, N., Headd, J. J., Hung, L.-W., Kapral, G. J., Grosse-Kunstleve, R. W., McCoy, A. J., Moriarty, N. W., Oeffner, R., Read, R. J., Richardson, D. C., Richardson, J. S., Terwilliger, T. C. & Zwart, P. H. (2010). Acta Cryst. D66, 213–221. Web of Science CrossRef CAS IUCr Journals Google Scholar
Bartesaghi, A., Matthies, D., Banerjee, S., Merk, A. & Subramaniam, S. (2014). Proc. Natl Acad. Sci. USA, 111, 11709–11714. Web of Science CrossRef CAS PubMed Google Scholar
Brown, I. D. & McMahon, B. (2002). Acta Cryst. B58, 317–324. Web of Science CrossRef CAS IUCr Journals Google Scholar
Chang, S., Head-Gordon, T., Glaeser, R. M. & Downing, K. H. (1999). Acta Cryst. A55, 305–313. Web of Science CrossRef CAS IUCr Journals Google Scholar
Colliex, C., Cowley, J. M., Dudarev, S. L., Fink, M., Gjønnes, J., Hilderbrandt, R., Howie, A., Lynch, D. F., Peng, L. M., Ren, G., Ross, A. W., Smith, V. H., Spence, J. C. H., Steeds, J. W., Wang, J., Whelan, M. J. & Zvyagin, B. B. (2006). Electron diffraction. In International Tables for Crystallography, Vol. C, edited by E. Prince, pp. 263–281. Dordrecht: Kluwer. Google Scholar
Doyle, P. A. & Turner, P. S. (1968). Acta Cryst. A24, 390–397. CrossRef IUCr Journals Web of Science Google Scholar
Fioravanti, E., Vellieux, F. M. D., Amara, P., Madern, D. & Weik, M. (2007). J. Synchrotron Rad. 14, 84–91. Web of Science CrossRef CAS IUCr Journals Google Scholar
Frisch, M. J. et al. (2016). Gaussian 09, Revision A.02, Gaussian Inc., Wallingford, CT, USA. Google Scholar
Fromm, S. A., Bharat, T. A., Jakobi, A. J., Hagen, W. J. & Sachse, C. (2015). J. Struct. Biol. 189, 87–97. Web of Science CrossRef PubMed Google Scholar
Hau-Riege, S. P. (2007). Phys. Rev. A, 76, 042511. Google Scholar
Hirai, T., Mitsuoka, K., Kidera, A. & Fujiyoshi, Y. (2007). J. Electron Microsc. 56, 131–140. Web of Science CrossRef Google Scholar
Hryc, C. F., Chen, D. H., Afonine, P. V., Jakana, J., Wang, Z., Haase-Pettingell, C., Jiang, W., Adams, P. D., King, J. A., Schmid, M. F. & Chiu, W. (2017). Proc. Natl Acad. Sci. USA, 114, 3103–3108. CrossRef CAS Google Scholar
Jönsson, P., Gaigalas, G., Bieroń, J., Froese Fischer, C. & Grant, I. P. (2013). Comput. Phys. Commun. 184, 2197–2203. Google Scholar
Macchi, P. & Coppens, P. (2001). Acta Cryst. A57, 656–662. Web of Science CrossRef CAS IUCr Journals Google Scholar
Mitsuoka, K., Hirai, T., Murata, K., Miyazawa, A., Kidera, A., Kimura, Y. & Fujiyoshi, Y. (1999). J. Mol. Biol. 286, 861–882. Web of Science CrossRef PubMed CAS Google Scholar
Pachucki, K. & Komasa, J. (2004). Chem. Phys. Lett. 389, 209–211. CrossRef CAS Google Scholar
Peng, L.-M. (1998). Acta Cryst. A54, 481–485. Web of Science CrossRef CAS IUCr Journals Google Scholar
Peng, L.-M., Ren, G., Dudarev, S. L. & Whelan, M. J. (1996). Acta Cryst. A52, 257–276. CrossRef CAS Web of Science IUCr Journals Google Scholar
Rez, D., Rez, P. & Grant, I. (1994). Acta Cryst. A50, 481–497. CrossRef CAS Web of Science IUCr Journals Google Scholar
Schmidt, M. W., Baldridge, K. K., Boatz, J. A., Elbert, S. T., Gordon, M. S., Jensen, J. H., Koseki, S., Matsunaga, N., Nguyen, K. A., Su, S., Windus, T. L., Dupuis, M. & Montgomery, J. A. (1993). J. Comput. Chem. 14, 1347–1363. CrossRef CAS Web of Science Google Scholar
Schmidt, P. C. & Weiss, A. (1979). Z. Naturforsch. 34, 1471–1481. CrossRef Google Scholar
Su, Z. & Coppens, P. (1997). Acta Cryst. A53, 749–762. CrossRef CAS Web of Science IUCr Journals Google Scholar
Waasmaier, D. & Kirfel, A. (1995). Acta Cryst. A51, 416–431. CrossRef CAS Web of Science IUCr Journals Google Scholar
Wang, J., Smith, V. H., Bunge, C. F. & Jáuregui, R. (1996). Acta Cryst. A52, 649–658. CrossRef CAS Web of Science IUCr Journals Google Scholar
Watson, R. E. (1958). Phys. Rev. 111, 1108–1110. CrossRef CAS Web of Science Google Scholar
Yamashita, H. & Kidera, A. (2001). Acta Cryst. A57, 518–525. Web of Science CrossRef CAS IUCr Journals Google Scholar
Yonekura, K. & Maki-Yonekura, S. (2016). J. Appl. Cryst. 49, 1517–1523. CrossRef CAS IUCr Journals Google Scholar
Yonekura, K., Kato, K., Ogasawara, M., Tomita, M. & Toyoshima, C. (2015). Proc. Natl Acad. Sci. USA, 112, 3368–3373. Web of Science CrossRef CAS PubMed Google Scholar
Zhong, S., Dadarlat, V. M., Glaeser, R. M., Head-Gordon, T. & Downing, K. H. (2002). Acta Cryst. A58, 162–170. Web of Science CrossRef CAS IUCr Journals Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.