BioMed Research International

Volume 2014 (2014), Article ID 904381, 8 pages

http://dx.doi.org/10.1155/2014/904381

Research Article

Crystallographic Analysis of Ground and Space Thermostable T1 Lipase Crystal Obtained via Counter Diffusion Method Approach

Sayangku Nor Ariati Mohamad Aris,1,2 Adam Leow Thean Chor,1,3 Mohd Shukuri Mohamad Ali,1,4 Mahiran Basri,1,5 Abu Bakar Salleh,1,4 and Raja Noor Zaliha Raja Abd. Rahman1,2

¹Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

²Department of Microbiology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

³Department of Cell and Molecular Biology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

⁴Department of Biochemistry, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

⁵Department of Chemistry, Faculty of Science, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

Received 1 April 2013; Revised 28 September 2013; Accepted 10 October 2013; Published 2 January 2014

Academic Editor: Michael Kalafatis

Copyright © 2014 Sayangku Nor Ariati Mohamad Aris 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.

Abstract.

Three-dimensional structure of thermostable lipase is much sought after nowadays as it is important for industrial application mainly found in the food, detergent, and pharmaceutical sectors. Crystallization utilizing the counter diffusion method in space was performed with the aim to obtain high resolution diffracting crystals with better internal order to improve the accuracy of the structure. Thermostable T1 lipase enzyme has been crystallized in laboratory on earth and also under microgravity condition aboard Progress spacecraft to the ISS in collaboration with JAXA (Japanese Aerospace Exploration Agency). This study is conducted with the aims of improving crystal packing and structure resolution. The diffraction data set for ground grown crystal was collected to 1.3 Å resolution and belonged to monoclinic C2 space group with unit cell parameters  Å,  Å, and  Å, whereas the diffraction data set for space grown crystal was collected to 1.1 Å resolution and belonged to monoclinic C2 space group with unit cell parameters  Å,  Å, and  Å. The major difference between the two crystal growth systems is the lack of convection and sedimentation in microgravity environment resulted in the growth of much higher quality crystals of T1 lipase.

1. Introduction

Lipases (triacylglycerol acylhydrolase) are enzymes which are present in many different organisms. It catalyze both the hydrolysis of triglycerides and the synthesis of esters formed from alcohol and long chain fatty acids. Enzymes from thermophiles often show more stability towards organic solvents and exhibit higher activity at elevated temperatures. Thus, this thermostable enzyme became an important requirement as a biocatalyst in industry including detergent, food, pharmaceuticals, cosmetic, textiles, fine chemicals, and biodiesel [1]. To better understand the mechanism of the enzyme, it is important to elucidate the three-dimensional structure of the enzyme. Performing precise X-ray data collection from high quality single crystal is the suitable way to solve the 3D structure accurately. The bottleneck of X-ray crystallography is obtaining high quality crystal which is influenced by the protein’s purity and how the crystals are grown [2]. An excellent crystallization method called the counter diffusion method was introduced by García-Ruiz and Moreno [3]. This technique requires a convection-free environment, which can be achieved using either gelled solutions, very thin capillaries, or microgravity conditions. The time for crystal growth may be extended by using longer gel tube. Furthermore, it is also possible to obtain larger crystals suitable for X-ray diffraction if capillaries with larger diameter are used [4]. There are two main benefits of growing crystals in counter diffusion method rather than vapour diffusion. The first benefit is the convection disturbances such as the effect of solution density or the fluctuating temperature which can be better controlled. The second benefit is the absence of crystal sedimentations. Sedimentation effect will result in the formation of crystals which usually grew at the bottom of the growth cell. By performing the crystallization inside a capillary, it can minimize or even suppress convection in solution on earth and create diffusive environment where capillary forces counterbalance gravity forces [5]. One of the challenges to get higher quality crystals is the natural convection that exists under normal earth gravity. On earth, gravity often has a negative impact on growing protein crystals. In microgravity, however, gravitational disturbances are at minimum level, thus allowing some crystals to grow in a more regular and perfect manner. Approach of crystallization in space using counter diffusion method has been conducted in order to reduce buoyancy-driven convection and avoid sedimentation of crystals. It will improve the yield of high quality crystals with better internal order which is suitable for diffraction analysis and thus give more accurate and reliable protein structure. An important objective for determining the three-dimensional structure of protein is to understand the structural mechanism and their biological activity. The ultimate goal of this study is to improve the crystal quality which will result in more accurate three-dimensional protein structures. Thus it will lead to a more precise understanding of biological function.

2. Materials and Methods

2.1. Materials

Resins for purification were purchased from GE Healthcare (Sweden). Capillary C-Tube Crystal-Tube Kit was purchased from Confocal Science Inc.

2.2. Purification of T1 Lipase

Geobacillus zalihae strain T1 lipase was overexpressed in pGEX vector in the prokaryotic system. Purification of T1 lipase was conducted according to the procedures of Leow et al. [6]. Mature T1 lipase was then purified via additional step, anion exchange chromatography using Q-sepharose HP column. The protein was eluted with 25 mM Tris-HCl (pH 9) containing linear gradient of 0 to 0.1 M sodium chloride at a flow rate of 0.5 ml/min. This step was applied to remove the remaining trace of impurities from earlier step. The purity of the fractions with high lipase activity was analysed by SDS-PAGE [7] and Native-PAGE. The desired protein was then desalted and kept in the storage buffer (25 mM Tris-HCl, pH 8).

2.3. Crystallization of T1 Lipase

Previously, crystallization of T1 lipase was achieved by sitting drop vapour diffusion method using reservoir solution composed of 1 M NaCl, 0.1 M NaH₂PO₄, 0.1 M KH₂PO₄, and 0.1 M MES pH 6.5 [8]. In the experiment described here, gel-tube method [4] which was modified from the original capillary counter diffusion method of García-Ruiz and Moreno [3] was applied for crystallization in order to control and improve the crystal quality. The glass capillaries of 0.5 mm in diameter and 60 mm in length from Crystal-Tube Kit (Confocal Science Inc.) were loaded with approximately 8 μL mixture of protein solution (5 mg/mL) and precipitant of 1 : 1 ratio. The upper end of the capillary was sealed with clay. The other end was inserted into the presoaked gel. The gel tube was cut to a length of 8 mm after being equilibrated by soaking it into the formulation of reservoir. The capillary was then placed in the cylinders of the syringe cases and the bottom of the cylinder cases was covered with its caps using a forceps. The caps of the cylinders of the syringe cases were carefully sealed. The area around the capillaries that is attached to the cylinder cases and the upper end of the capillaries were sealed with Araldite adhesive (Huntsman Advance Material). The crystallization of T1 lipase on the ISS was carried out in the Japanese Experiment Module Kibo of the JAXA (Japanese Aerospace Exploration Agency). The space sample unit aboard the Progress spacecraft to ISS. The canister containing the space unit samples was mounted to the experiment facility, the PCRF (Protein Crystallization Research Facilities) by a crew member. The temperature for both space and ground control was set to 20°C.

2.4. X-Ray Diffraction Analysis

Crystals of T1 lipase were mounted in cryoloop soaked in cryoprotectant solution consisting of 40% glycerol and then flash-cooled in a liquid nitrogen stream at 100 K. High resolution X-ray diffraction data for both space grown crystal and ground control of T1 lipase crystals were collected using Rayonix MAR225HE CCD detector at station BL41XU of SPring-8 (Harima, Japan). The detector was set at a distance of 130 mm with 0.5 seconds per frame exposure and a 0.5˚ per frame oscillation was used. All resulting data from this experiment were integrated and scaled using the programs DENZO and SCALEPACK from the HKL-2000 suite [9].

2.5. Structure Determination

Structure refinement and model building were performed by using CCP4 program suite (Collaborative Computational Project, Number 4 and COOT (Crystallographic Object-Oriented Toolkit). Meanwhile, structural analysis was carried out using the YASARA (Yet Another Scientific Artificial Reality Application) software version 10.2.1. [10].

3. Results and Discussion

T1 lipase was successfully purified to homogeneity with a molecular weight of 43 kDa. The purity was confirmed by the presence of a single band on SDS-PAGE and Native-PAGE (Figure 1). T1 lipase was extensively purified via ion exchange chromatography to gain higher purity and to improve its uniformity. The purity and uniformity of the protein sample are paramount important in crystallization. Incorporation of the impurities or nonuniformity in the crystal will cause lattice strain and reduce the crystalline order, which will eventually decrease the quality of the diffraction pattern and thus become unsuitable for X-ray structural analysis [11]. To date, very few thermostable lipase crystals had been reported. Crystals from thermostable Bacillus stearothermophilus P1 lipase [12], Bacillus thermocatenulatus lipase 2 (BTL2) [13], and Bacillus stearothermophilus L1 (L1 lipase) [14] were obtained from hanging drop vapour diffusion method. Currently, only two thermostable lipases were reported grown by counter diffusion method, which are L42 lipase and L2 lipase. X-ray data of L42 lipase were collected at 2.0  and the crystal belonged to the monoclinic space group [15]. However thermostable L2 lipase crystal belonged to orthorhombic space group with resolution 2.7  [16]. Maruki et al. [17] have acquired a high resolution structure of a new crystal form of human triosephosphate isomerase (TIM) by employing counter diffusion method crystallization in microgravity environment. Figure 2 shows T1 lipase crystal obtained by counter diffusion method (ground and space grown). The dimensions of the crystals were approximately  mm and  mm for ground-grown crystal and space-grown crystal, respectively. Data sets for T1 lipase crystals were collected at the beamline of the BLX41 at SPring-8 Synchrotron radiation facilities, Japan, to a resolution of 1.3  for ground-grown crystal and 1.1  for space-grown crystal. Table 1 shows X-ray data statistics for T1 lipase crystals grown under microgravity and on the earth. The reflections of ground-grown crystal was indexed on a centred monoclinic lattice (C2 with unit cell parameters , , , , , and . Space-grown crystal belonged to the same system with unit cell parameters , , , , , and . Mosaicity measurements and the resolution of the diffraction are the most important indicators for crystal quality [18]. Space-grown crystals display a better signal to noise ratio (I/Sigma) rather than ground-grown crystals. It shows that microgravity can provide a better environment for growing good quality crystals which are difficult to obtain on the ground. The quantity and quality of the diffraction data will affect the electron density map. The difference in resolution will be reflected in the quality of the electron density map, where the model of three-dimensional protein structures is built. A higher resolution of the data results in higher resolution of the electron density map. It will give higher accuracy of the positions of the atoms in the structure. From the experiment, the final 2Fo-Fc electron density maps which are of excellent quality in both crystals are further improved in space-grown crystals. This is due to the higher resolution of space-grown crystal. The positions of the atoms are much better defined. Well-defined side chains of amino acids and the greater number of ordered water molecules are emphasized, thus resulting in more accurate structure [19]. Overall, the polypeptide backbone can be easily fit into the electron density observed with space data. Figure 3 shows the snapshots of electron density from space-grown and ground-grown crystal obtained after final refinement. All atoms in the protein model were best fit into the electron density map. The crystal has a predicted solvent content of 55.11% and a Matthews’s coefficient value () of 2.74 Da⁻¹ for ground grown crystal and solvent content of 55.01% and a Matthews coefficient value () of 2.73 Da⁻¹ for space grown crystal. The and solvent content values of T1 lipase are in the range, as Matthews observed that the solvent content in protein crystal ranged from 27% to 65% with an average of 43% [20]. According to Kantardjieff and Rupp, protein crystals with less solvent tend to diffract better [21]. It shows that tightly packed crystals tend to diffract better than the loosely packed ones. In early stages of crystal structure determination, it is important to know the solvent content which helps to determine the number of molecules in the asymmetric unit [20]. The atomic coordinates Geobacillus zalihae T1 lipase (Protein Data Bank under accession code 2DSN) without water molecules was used as a template for the refinement of the ground-grown and space-grown crystal structures. The structural refinement was performed using the molecular replacement method using CCP4 program suite (Collaborative Computational Project, number 4). Refinement was performed by alternating rounds of model building using the program COOT (Crystallographic Object-Oriented Toolkit) and crystallographic refinement using REFMAC5. The quality and accuracy of the refined model was evaluated using Ramachandran plot [22] and ERRAT. Analysis of B-factor will provide information of the protein dynamics, flexibility of amino acids, and protein stability [23]. Crystal packing interactions make a significant contribution to the value of B-factor [24]. The B-factor (atomic displacement parameter) in protein crystal structures reflects the fluctuation of an atom about its average position. The distribution of B-factors along a protein sequence is regarded as an important indicator of the protein’s structure, reflecting its flexibility and dynamicity, where the low value of B-factors correspond to well-defined parts of the structure, whereas high value of B-factors might indicate highly disordered parts of the structure or even misinterpreted parts of the model [25]. Figure 4 shows that Wilson plot generated from CCp4 software gives overall value B-factor of 10 A² and 5.6 A² for ground-grown and space-grown crystal structure, respectively. Overall, ground-grown crystal structure shows higher B-factor value, indicating that it has the higher flexibility compared to the space-grown structures. The distribution of average B-factors over the peptide main-chain for each molecule for both structures are shown in Figure 5. Table 2 shows refinement statistics for ground-grown and space-grown crystals. Both refined model contains two molecules (chain A and chain B) per asymmetric unit and each comprises 387 amino acids (residues Ser2 to Pro388) (Figure 3). Both solved structures show a common / hydrolase fold that comprises the conserved Ser-His-Asp catalytic triad (Ser113-His358-Asp317) (Figure 6). The arrangement of -helix 6 and -helix 7 was shown as “lid” that covers the active site. The presence of hydrophobic region of -helix 6 (Phe176, Phe180, Phe181, Leu183, Ala186, Val187, Leu188, and Ala190) provided a hydrophobic surroundings to the active site. This environment may play an important role in enzyme catalysis. The two molecules in the asymmetric unit are highly similar to each other. Main backbone of chain A and chain B for both ground- and space-grown crystal was superimposed giving root mean square deviation (RMSD) value of 0.2185  for chain A and 0.4214  for chain B, respectively (Figure 7).

Table 1: X-ray data collection statistics for T1 lipase.

Ground
Space

Space group

Unit cell parameters (Å)

Resolution range (Å)
30.0–1.30
40–1.10

0.103 (0.00)
0.063 (0.46)

Data completeness (%)
99.6 (99.1)
98.1 (94.1)

I/Sigma (I)
35.02 (3.93)
24.23 (2.0)

Mosaicity range
0.12–0.21
0.14–0.20

Molecule per asymmetric unit
2
2

Matthews coefficient (Å³/Da)
2.74
2.73

Solvent content (%)
55.11
55.01

Values in bracket refer to the highest resolution shell.

Table 2: Final statistics of refined structure for the ground-grown and space-grown T1 lipase crystals.

Ground-grown crystal
Space-grown crystal

Resolution range used in refinement (Å)
30–1.2
30–1.2

Number of reflections used in refinement
272755
271607

R factor
0.134
0.129

R free
0.162
0.150

Number of protein atoms
387
387

Number of water molecules
721
781

Number of metal ions
7
7

Number of glycerols
4
4

RMS deviation

Bond length (Å)
0.021
0.020

Angles (Deg)
1.882
1.822

Planes (Å)
0.014
0.014

Ramachandran plot details

Residues in favoured regions (%)
91.5
91.8

Residues in additional allowed regions (%)
7.9
7.6

Residues in generously allowed regions (%)
0.3
0.3

Residues in disallowed regions (%)
0.3
0.3

Average B-factor value (A 2)

Chain A
10.8
9.4

Chain B
12.1
10.5

Figure 1: Detection of T1 lipase on gel electrophoresis. (a) SDS-PAGE of T1 lipase after ion exchange chromatography (lane: M: marker, 1–4: purified T1 lipase). (b) The homogeneity was confirmed via Native-PAGE. (lane: 1–4: purified T1 lipase).

(a)

(b)

Figure 2: Ground-grown (a) and space-grown (b) crystals of T1 lipase in capillary using counter diffusion method.

Figure 3: Snapshots of ground-grown crystal (1.1 ) and space-grown crystal (1.3 ) electron density map countered at 1 sigma level displayed in COOT.

(a)

(b)

Figure 4: Wilson plot for (a) ground-grown structure and (b) space-grown structure.

(a)

(b)

Figure 5: Plot of the average B-factor per residue for main chain of each structure. (a) Ground-grown crystal structure and (b) space-grown crystal structure.

Figure 6: Ribbon diagram for overall structure. The structure of T1 lipase is shown as a ribbon diagram with -helices (blue), -strands (red), 3¹⁰ helix (yellow), turn (green), and coil (cyan). It comprises the catalytic triad Ser113, Asp317, and His358. (Comparison of molecules A and B for both structures shows that the overall structures of these two molecules are almost identical. Therefore, only one of them is shown.)

(a)

(b)

Figure 7: Superimpose of chains A and B for ground-grown and space-grown of T1 lipase. (a) Superimpose of chain A for both T1 lipase crystal (magenta for space crystal) with RMSD of 0.2185 . (b) Superimpose of chain B (gray for space crystal) with RMSD of 0.4214 .

4. Conclusion

Comparative crystallographic analysis revealed that enhanced quality of the diffraction data was derived from space-grown crystal. Compared to ground-grown crystal, the space-grown crystal diffracts better. As a result, an improved initial electron density map was obtained for modeling and refinement, contributing to clearer structural information of the enzyme. Strictly comparative crystallographic analysis reveals that space-grown crystals are better compared to ground-crystal. Analysis of data collection and refinement statistics showed that crystallization in counter diffusion method using microgravity environment have highly improved the internal order of crystals and thus gave a more precise three-dimensional structure. A more in depth analysis of both structures can help in the understanding of the enzyme.

Conflict of Interests

The authors hereby declare that this paper has no conflict of interests from any party.

Acknowledgments

The authors thank Professor Yoshiki Higuchi, Dr. Koji Inaka, and Dr. Naoki Shibata for the helpful discussion on data analysis and JAXA (Japan Aerospace Exploration Agency) for X-ray data collection. This work was performed under the International Collaborative Research Program of the Institute for Protein Research, Osaka University. This research was funded by the Ministry of Science, Technology, and Innovation, Malaysia (09-05-MGI-GMB001).

References

F. Hasan, A. A. Shah, and A. Hameed, “Industrial applications of microbial lipases,” Enzyme and Microbial Technology, vol. 39, no. 2, pp. 235–251, 2006.

F. Otalora, J. A. Gavira, J. D. Ng, and J. M. Garcia-Ruiz, “Counter diffusion methods applied to protein crystallization,” Progress Biophysic Molecular Biology, vol. 101, pp. 26–37, 2009.

J. M. García-Ruiz and A. Moreno, “Investigations on protein growth by the gel acupuncture method,” Acta Crystallographica D, vol. 50, pp. 484–490, 1994.

H. Tanaka, K. Inaka, S. Sugiyama et al., “A simplified counter diffusion method combined with a 1-D simulation program for optimizing crystallization conditions,” Journal of Synchrotron Radiation, vol. 11, no. 1, pp. 45–48, 2004.

T. Y. Miller, X.-M. He, and D. C. Carter, “A comparison between protein crystals grown with vapor diffusion methods in microgravity and protein crystals using a gel liquid-liquid diffusion ground-based method,” Journal of Crystal Growth, vol. 122, no. 1-4, pp. 306–309, 1992.

T. C. Leow, R. N. Z. R. A. Rahman, M. Basri, and A. B. Salleh, “A thermoalkaliphilic lipase of Geobacillus sp. T1,” Extremophiles, vol. 11, no. 3, pp. 527–535, 2007.

U. K. Laemmli, “Cleavage of structural proteins during the assembly of the head of bacteriophage T4,” Nature, vol. 227, no. 5259, pp. 680–685, 1970.

T. C. Leow, R. N. Z. R. A. Rahman, A. B. Salleh, and M. Basri, “High-temperature crystallization of thermostable T1 lipase,” Crystal Growth and Design, vol. 7, no. 2, pp. 406–410, 2007.

Z. Otwinowski and W. Minor, “Processing of X-ray diffraction data collected in oscillation mode,” Methods in Enzymology, vol. 276, pp. 307–326, 1997.

E. Krieger, G. Koraimann, and G. Vriend, “Increasing the precision of comparative models with YASARA NOVA—a self-parameterizing force field,” Proteins: Structure, Function and Genetics, vol. 47, no. 3, pp. 393–402, 2002.

M. Benvenuti and S. Mangani, “Crystallization of soluble proteins in vapor diffusion for x-ray crystallography,” Nature protocols, vol. 2, no. 7, pp. 1633–1651, 2007.

S. Sinchaikul, J. D. A. Tyndall, L. A. Fothergill-Gilmore et al., “Expression, purification, crystallization and preliminary crystallographic analysis of a thermostable lipase from Bacillus stearothermophilus P1,” Acta Crystallographica D, vol. 58, no. 1, pp. 182–185, 2002.

C. Carrasco-López, C. Godoy, B. De Las Rivas et al., “Crystallization and preliminary X-ray diffraction studies of the BTL2 lipase from the extremophilic microorganism Bacillus thermocatenulatus,” Acta Crystallographica F, vol. 64, no. 11, pp. 1043–1045, 2008.

S.-T. Jeong, H.-K. Kim, S.-J. Kim, J.-G. Pan, T.-K. Oh, and S.-E. Ryu, “Crystallization and preliminary X-ray analysis of a thermoalkalophilic lipase from Bacillus stearothermophilus L1,” Acta Crystallographica D, vol. 57, no. 9, pp. 1300–1302, 2001.

M. S. Khusaini, R. N. Z. R. A. Rahman, M. S. Mohamad Ali, T. C. Leow, M. Basri, and A. B. Salleh, “Crystallization and preliminary X-ray crystallographic analysis of a thermostable organic solvent-tolerant lipase from Bacillus sp. strain 42,” Acta Crystallographica F, vol. 67, no. 3, pp. 401–403, 2011.

F. M. Shariff, R. N. Z. R. A. Rahman, M. S. M. Ali, A. L. T. Chor, M. Basri, and A. B. Salleh, “Crystallization and preliminary X-ray crystallographic analysis of highly thermostable L2 lipase from the newly isolated Bacillus sp. L2,” Acta Crystallographica F, vol. 66, no. 6, pp. 715–717, 2010.

R. Maruki, M. Warizaya, T. Kinoshita, H. Nakajima, and S. Nishimura, “Structure of a high-resolution crystal form of human triosephosphate isomerase: Improvement of crystals using the gel-tube method,” Acta Crystallographica F, vol. 61, no. 4, pp. 346–349, 2005.

S. S. Terzyan, C. R. Bourne, P. A. Ramsland, P. C. Bourne, and A. B. Edmundson, “Comparison of the three-dimensional structures of a human Bence-Jones dimer crystallized on earth and aboard US Space Shuttle Mission STS-95,” Journal of Molecular Recognition, vol. 16, no. 2, pp. 83–90, 2003.

J. D. Ng, C. Sauter, B. Lorber, N. Kirkland, J. Arnez, and R. Giegé, “Comparative analysis of space-grown and earth-grown crystals of an aminoacyl-tRNA synthetase: space-grown crystals are more useful for structural determination,” Acta Crystallographica D, vol. 58, no. 4, pp. 645–652, 2002.

B. W. Matthews, “Solvent content of protein crystals,” Journal of Molecular Biology, vol. 33, no. 2, pp. 491–497, 1968.

K. A. Kantardjieff and B. Rupp, “Matthews coefficient probabilities: improved estimates for unit cell contents of proteins, DNA, and protein-nucleic acid complex crystals,” Protein Science, vol. 12, no. 9, pp. 1865–1871, 2003.

R. A. Laskowski, M. W. MacArthur, D. S. Moss, and J. M. Thornton, “PROCHECK—a program to check the stereochemical quality of protein structures,” Journal of Applied Crystallography, vol. 26, pp. 283–291, 1993.

Z. Yuan, J. Zhao, and Z.-X. Wang, “Flexibility analysis of enzyme active sites by crystallographic temperature factors,” Protein Engineering, vol. 16, no. 2, pp. 109–114, 2003.

D.-W. Li and R. Brüschweiler, “All-atom contact model for understanding protein dynamics from crystallographic B-factors,” Biophysical Journal, vol. 96, no. 8, pp. 3074–3081, 2009.

Z. Yuan, T. L. Bailey, and R. D. Teasdale, “Prediction of protein B-factor profiles,” Proteins, vol. 58, no. 4, pp. 905–912, 2005.