Abstract
We compared the topologies of protein and small molecule crystals, which have many common features â both are molecular crystals with intermolecular interactions much weaker than intramolecular interactions. They also have different features â a considerably large fraction of the volume of protein crystals is occupied by liquid water while no room is available to other molecules in small molecule crystals. We analyzed the overall and local topology and performed multilevel topological analyses (with the software package ToposPro) of carefully selected high quality sets of protein and small molecule crystal structures. Given the suboptimal packing of protein crystals, which is due the special shape and size of proteins, it would be reasonable to expect that the topology of protein crystals is different from the topology of small molecule crystals. Surprisingly, we discovered that these two types of crystalline compounds have strikingly similar topologies. This might suggest that molecular crystal formations share symmetry rules independent of molecular dimension.
Similar content being viewed by others
Introduction
According to McPherson1, the first observation of a protein crystal - earthworm hemoglobin - is about 150 years old2. However, during several decades, protein crystallization was just a technique for protein purification from complex extracts. Universal interest in protein crystallization arose stunningly only one century later, when Kendrew and Perutz proved that it was possible to determine crystal structures of proteins (sperm whale myoglobin and then hemoglobin)3,4,5. Thereinafter, thousands of protein crystals have been produced and about 120,000 protein crystal structures are now available in the Protein Data Bank6,7.
Protein crystallization is considered the bottleneck of protein crystallography. Proteins are extremely soluble, in physiological conditions, and an âevolutionary negative designâ has been proposed to have determined the difficulty to crystallize proteins8. However, protein crystallization in vivo has been observed and it allowed the determination of the first protein crystal structure in vivo 9.
Although they are crucial in macromolecular crystallography, protein crystals attracted only little attention in recent years. Few studies have been dedicated to the analysis of protein-protein interactions at crystal contacts and their comparison with physiological interactions10,11,12. Later, a survey on the frequency of crystal packing interactions appeared13. The importance of packing bridges â small molecules or metal cations that bridge two adjacent and symmetry related macromolecules â in protein crystallization was also discovered14,15.
Protein crystals are suboptimally packed molecular crystals, because of their special shape and size. They contain a substantial number of cavities and channels filled by liquid water. It is believed that, roughly, 20â80% of the crystal volume is made by liquid water16,17,18. This is the reason why enzymes can catalyze chemical reactions in the solid state19,20 and why protein crystals can be soaked to prepare heavy atom derivatives that can be used to solve the phase problem21. Therefore, beside water molecules that strongly interact with the protein molecules22,23,24, protein crystals contain a considerable amount of aqueous solution, containing several types of molecules that were present in the crystallization cocktail. Consequently, crystallization free energy of proteins is very modest, only 3â6âkcal/mole, relative to the dissolved state in solution25.
The topology of packing in protein crystals may be different from that of other molecular solids, where the packing is, in general, optimal and where there are no internal cavities filled with amorphous material. It is thus interesting to compare the topology of small organic and protein crystals. The most widespread model of molecular crystals is the Kitaigorodskiiâs model26, according to which each molecule in a homomolecular crystal tends to be surrounded by 12 other molecules; other coordination numbers (10 or 14) occur much rarer. However, a comprehensive treatment of all molecular structures stored in the Cambridge Structural Database using a computer-realized Voronoi approach27,28 showed that coordination number 14 is more abundant than 12. This was explained by the model of deformable molecules, in contrast to the model of rigid ones used by Kitaigorodskii.
Since the approach used in references27,28 can be applied to the molecules of any size and shape it was natural to use it for comparative analysis of packings of protein and small molecules.
Experimental
Data selection
Protein data
All protein crystal structures were taken from the Protein Data Bank6,7, according to the following criteria. We discarded structures containing nucleic acids and structures with an average B-factor smaller than 10 à 2 or larger than 40 à 2. We considered only structures with 50â500 amino acids in the asymmetric unit, with resolution better than 3âà , with working R-factor better than 0.25, and determined at 80â120âK. We discarded structures deposited without their experimental diffraction data. Redundancy was reduced to 30% sequence identity. Structures with missing protein atoms or with protein atoms deposited with zero occupancy were discarded as well as the structures where non-protein and non-water atoms are more than 5% of the total number of atoms.
We note that these are standard criteria for protein crystal structure selection from the Protein Data Bank29,30. In particular, resolution and R-factor thresholds ensure that unreliable structures are excluded from the analyzed data sets; moreover, we considered a homogeneous set containing only low temperature crystal structures, which are now the routine, since packing might depend on temperature; furthermore, limitations on the average B-factor and on protein dimension ensure that anomalous structures â much larger or smaller or much more or less flexible than customary proteins â are removed from the data sets; and eventually, redundancy reduction to 30% sequence identity, though very crude, is necessary to avoid distribution biases.
It is also important to observe that the exclusion of structures containing too many (more than 5%) hetero-atoms ensures that crystal contacts due to the presence of co-crystallized small molecules (the so called packing bridges) are minimized. Analogously, the exclusion of incomplete structures, where some of the protein atoms/residues were undetected, ensures that no protein-protein crystal packing contacts are neglected.
We prepared three ensembles of protein crystal structures. In the first, we collected only monomeric proteins that crystallized with only one molecule per asymmetric unit (monomer set; 394 structures). In the second, we grouped dimeric proteins that crystallized with only one dimer per asymmetric unit (dimer set; 207 structures). In the third and last ensemble, we pulled together structures of monomeric proteins with two molecules in the asymmetric unit (double set; 164 structures).
These data sets might seem small if compared to the number of entries deposited in the Protein Data Bank. However, it is necessary to consider that only 31% of the PDB files contain monomeric proteins with one molecule per asymmetric unit; only 7% contain dimeric proteins with one molecule per asymmetric unit; and only 9% contain monomeric proteins with two molecules per asymmetric unit. Moreover, redundancy reduction to 30% maximal sequence identity implies a 6-fold decrease of the number of PDB entries. Eventually, in the large majority of the protein crystal structures deposited in the Protein Data Bank there are residues that have not been seen experimentally31. Since they are usually at the protein surface, it is not possible to use these structures to perform reliable analyses of the crystal packing. We note also that the quality of the data extracted from the Protein Data Bank is becoming a crucial issue in structural bioinformatics29,32. It is also important to observe that although these are high quality and non-redundant datasets of protein crystal structures, it is impossible to estimate the extent to which they represent the entire protein âuniverseâ, a considerable fraction of which is constituted, for example, by membrane proteins that are under-represented in the Protein Data Bank33.
The identification codes of all these crystal structures are listed in the Supplementary Information (Table S1).
Small molecules data
Crystallographic data for organic molecular crystals consisting of chemically equivalent molecules were taken from the CSD (release 5.38). We have not analyzed molecular crystals of coordination or organometallic compounds, which have different chemical nature compared to proteins. In total, we have considered all 166,297 structures, which contain one independent molecule and 10,240 structures with two independent chemically equivalent molecules in the asymmetric unit. The entries containing incomplete, erroneous data, or disordered structures as well as those with Rfâ>â10% were excluded; no other restrictions were applied.
Crystal contacts identification
Protein data
The crystal contacts were identified as described previously13. The asymmetric unit was transformed according to all symmetry operations and translated up to three times along the axes a, b, and c in both negative and positive directions. Two molecules â one, the reference molecule, is in the asymmetric unit and the other, the satellite, is symmetry related to the first â were considered to be in contact if at least one atom of one of them is closer than 4.5âà from an atom of the other one. We refer to this method below as âDistanceâ method.
Small molecules data
All interatomic contacts were determined by means of the âDomainsâ method34 implemented into the ToposPro package35. This method is based on the Voronoi partition and accounts not only for interatomic distances, but also for atomic surrounding and screening effects. An interatomic contact A-B is assumed to exist if there is a face of an atomic Voronoi polyhedron with a solid angle, Ωâ>â1.5% of the total solid angle 4Ï steradian and the line A-B crosses the face (the contact is direct). To distinguish valence and non-valence contacts a set of additional criteria was used as described in reference34. The âDomainsâ method is designed to process large samples of crystal structures in an automated mode, which was crucial for our study. At the same time, this method was never used for proteins; thus, it was important to check if this method gives similar results for proteins as the âDistanceâ method.
The molecular coordination number was computed for the molecular Voronoi polyhedron, which is a union of Voronoi polyhedra of all atoms of the molecule (Fig. 1)27,28.
External faces of the molecular Voronoi polyhedron correspond to intermolecular contacts, and two molecules are considered interacting if there is at least one intermolecular contact between them (Fig. 2)27.
This method was also used to compute coordination numbers of some protein molecules for comparison with the simple geometrical approach described above.
Topological analysis
To analyze molecular packings of both proteins and small molecules we squeezed each molecule into its center of gravity to obtain the so-called underlying net, i.e. a net, which keeps the information on the connection of molecules throughout the crystal, and ignores the information about their geometrical properties and internal structure (Fig. 3)36. As a result, the crystal structure is represented as an infinite graph, whose vertices and edges correspond to molecules and intermolecular links, respectively. Such simplification of the structure allows us to put aside differences in size, shape and other geometrical features of molecular ensembles and to focus on the method of their local and overall connection, i.e. on the crystal architecture as a whole. This is a way to find correlations between packings of so different objects as small molecules and proteins.
A molecule of the bubble protein 1uoy50, its center of mass and centers of mass of four neighboring molecules (pink balls), which represent a fragment of the underlying net.
The packings represented by the same underlying net are topologically equal and belong to the same isoreticular series, which is designated by a topological type symbol. There are several nomenclatures for topological types; here we use two of them: the RCSR three-letter codes37 and the NDn symbols used in the collection of ToposPro topological databases (TTD)38, where N is the series of coordination numbers for all independent nodes, D is one of the letters C, L or T, indicating the dimensionality of the net (C, L, T for chain, layer, and three-periodic framework, respectively) and n is the number of the net with a given ND sequence in the TTD collection. For example, the 16T4 symbol means a three-periodic net with 16-coordinated nodes (i.e. any molecule has 16 neighbors in the packing) and it is the fourth 16-coordinated net in the TTD collection.
Since the strength of intermolecular contacts is quite different, the topology of a molecular packing strongly depends on which contacts are taken into account. Ignoring interactions of a particular level of strength, we generate different underlying nets for the same crystal structure. Each underlying net specifies the packing topology at a given level of interaction and corresponds to a structure representation.
For small molecule crystals, we estimated the strength of an intermolecular interaction by the value of Ωmol that is the part of the sum of solid angles of all contacts, which are formed by a given molecule; this part corresponds to all contacts between a given pair of molecules. ToposPro includes a special procedure, which generates all possible representations by subsequently breaking groups of weak intermolecular contacts with close Ωmol (the difference between neighboring Ωmol in the group does not exceed 1.5% of the total solid angle 4Ï steradian). As a result the structure is considered at different levels of intermolecular interactions. The underlying net was constructed for each level of intermolecular interaction to obtain the information about the way of molecular connection. Such a multilevel analysis, which we successfully used for elucidation of organic molecular crystals39, allows one to separate underlying motifs, which determine the structure architecture. These motifs can be the same in the molecular packings, which differ in details of weaker intermolecular contacts (Fig. 4). In our study we considered only the representations with molecular coordination numbers â¤14; there are 105,549 molecular packings with such connectivity among the 166,297 structures studied. The topologies of molecular packings with a higher connectivity are non-specific and can be usually reduced to a known topology with a smaller coordination number after discarding weak intermolecular contacts.
An illustration of the multilevel topological analysis for the 2-(1,2-diacetylhydrazino)acetamide (ABELAW)51: (a) 12 coordinated underlying net, hcp topology; (b) 10 coordinated underlying net, bct topology; (c) 8 coordinated underlying net, hex topology; (d) 6 coordinated underlying net, hxl topology; (e) 4 coordinated underlying net, sql topology. The origin molecule is in red, the molecules, connected at different levels of interactions are in yellow.
For protein crystals, strengths of crystal packing contacts were estimated based on their dissociation free energies as described below.
Crystal contacts thermodynamics
Since protein crystallization occurs in aqueous solutions and crystals are largely impregnated by aqueous solutions, the formation and stability of crystal packing contacts cannot be considered out of solvent context. For this reason, we used the program PISA40 to estimate the stability of individual crystal packing contacts formed by the reference molecule and each of its satellites, which were generated by applying all the symmetry operations to the asymmetric unit and by applying three translations along all the axes in the negative and in the positive directions. Intermolecular contacts were defined as pairs of atoms, one in the reference molecule and the other in a satellite molecule, closer than 4.5âà .
According to the PISA definitions, the dissociation free energy ÎG diss is defined as:
where ÎG int is the binding free energy of the two molecules (the reference molecule plus one of its satellites), T is the absolute temperature, and ÎS is the entropy change upon association of the two molecules. If ÎGdissâ>â0, the crystal packing contact is thermodynamically unstable in solution and is present only in the crystal as an experimental artifact. ÎG int is defined as:
where ÎG solv(AB) is the solvation energy of the dimeric entity AB (A is the reference molecule and B is one of its satellites), ÎG solv(A) is the solvation energy of the reference molecule, ÎG solv(B) is the solvation energy of the satellite molecule, ÎG cont(AB) is the contact energy between molecules A and B, and ÎG es(AB) is the electrostatic energy between molecules A and B. Krissinel and Henrick developed a strategy to compute all these components necessary to estimate the ÎG diss values. This strategy, in part, is based on heuristic guesses of some terms, for example for the atomic solvation parameter, the hydrogen bond, salt bridge, and disulfide bond contributions40.
Results and Discussion
Overall topologies of molecular packings
Since packing contacts were determined with two different approaches in protein crystals (âDistance methodâ, see Methods) and in small molecule crystals (âDomains methodâ), it is mandatory to compare the results obtained by following these two different approaches and to verify that they produce identical results. For this purpose we have computed intermolecular contacts for some proteins with the âDomains methodâ. In most cases, the results of both methods coincide (Fig. 5a), however, even if they differ, this concerns only weak long contacts, which the âDomainsâ method can consider in contrast to the âDistanceâ method (Fig. 5b). Since the weak contacts are not considered in the multilevel topological analysis, the topologies obtained by the two methods can be assumed to be comparable.
Underlying nets of crystal structures of (a) 4a3X52 and (b) 2xvs53 determined by âDistanceâ (left) and âDomainsâ (right) methods. The additional weak contacts, which are determined by the âDomainsâ method, are shown by red lines. Distances between the molecular centers of mass are given in à .
Table 1 shows the 20 most frequent underlying nets for monomeric protein crystals as well as for small molecule crystals. One can see that small molecules tend to dense packings with coordination numbers ranging from 12 to 14, in agreement with previous results27,28. On the contrary, coordination numbers are very often much lower in protein crystals, in agreement with the fact that a considerable fraction of their volume is occupied by liquid water. However, all underlying nets of protein crystals in Table 1 are well known and almost all of them (except nets chb, smt, yfh) are fcu and hcp subnets. Interestingly, the 14-coordinated body-centered cubic packing (bcu-x), which is the most frequent in small molecule crystals is rather infrequent in protein crystals (Table 1). Thus protein crystals follow the close packing topologies but with some gaps, which are obviously filled by the solvent molecules and decrease coordination numbers of protein molecules.
It is interesting to observe that while few topologies are extremely more frequent than others in small molecule crystals, this is not observed in protein crystals. 30% of the small molecule crystals are associated with the most frequent topology (bcu-x) and 65% of them have one of the four most common topologies (bcu-x, tcg-x, 14T3, gpu-x). On the contrary, only 6% of the protein crystals are associated with the most frequent topology (fcu) and only 19% of them have one of the four most common topologies (fcu, bct, hex, vcs). This suggests that the suboptimal packing of the protein molecules in the crystal state allows a wider number of topologies and none of them can be much more frequent than the others.
Similar results are observed for the crystals of dimeric proteins where there is one dimer per asymmetric unit (Table 1) and for the small molecule crystals where there are two independent molecules in the asymmetric unit (Table 2). On the contrary, the underlying nets of crystals of monomeric proteins where the asymmetric unit contains two independent molecules set are very diverse and it was impossible to separate any preferable motif (Table 2).
Local topologies: the coordination numbers
Here the attention is focused on the coordination numbers, independently of the overall topology. This simple analysis of the individual coordination numbers shows a considerable difference between protein crystals and small molecule crystals. The most frequent coordination numbers for protein molecules are 7 and 8, while small molecules prefer higher coordination numbers equal to or larger than 14 (Fig. 6; see also supplementary Table S2).
Distribution of the coordination numbers (%) for all the structures of proteins examined (see full data in Table S2) compared to the CN for small molecules.
Like for the overall topology analysis, we see that the most frequent coordination number of small molecules (14) is observed in a large fraction of crystals (52%) while the most frequent coordination number for proteins (7 or 8) is observed only in a smaller fraction of crystals (about 20%). Moreover, small molecules can adopt 25 different coordination number values (from 4 to 32) while proteins can have only 14 different coordination number values (ranging from 1 to 16 â coordination numbers 2 and 15 are never observed).
All these observations point out that proteins, which are suboptimally packed in their crystals, can hardly by surrounded by numerous other proteins and that the palette of their coordination numbers is considerably more limited.
Rare overall topologies
Although the overall topologies are quite similar in protein and in small molecule crystals, protein molecules can form packings of a topology that never occurs in packings of small molecules. Among 394 monomeric and 207 dimeric protein structures there are 31 and 10 underlying nets, respectively, which were not found in the packings of small molecules, and all of them have a unique topology (see Table S3 in the Supplementary Material). Among the crystal structures with two independent protein molecules per asymmetric unit, there are 44, which have unique overall topologies never revealed in the small molecule crystals (Table 2).
This means that from 5 to 26% of protein molecules essentially differ from small molecules by their ability to be packed.
Multilevel topological analysis
In order to reach a better understanding of the differences between the topologies of protein and small molecule crystals, we performed the multilevel analysis for the molecular packings in 394 monomeric protein structures and 105,549 structures of small molecules. As a result, 1,799 and 1,144,539 structure representations were generated, respectively.
The underlying net topologies are summarized in Table S4 (see Supplementary Material). In general, the underlying motifs that are frequent in protein crystals are also frequent in small molecules crystals. The frequency of the protein topologies and the frequency of the small molecule topologies shown in Table S4 correlate well (Pearson correlation coefficientâ=â0.625; Spearman correlation coefficientâ=â0.563) as shown in the histogram in Fig. 7.
It appears that the underlying motifs of packings in both protein and small molecule crystals have simple and quite common topologies, such as sql, hcb, dia, pcu, hex (Fig. 8). The role of close-packed motifs (fcu, hcp, Fig. 8) is expectedly higher for small molecules than for macromolecules. Since the multilevel analysis reveals the basic motifs (skeletons), which are formed by the strongest intermolecular contacts, the results obtained mean that that the skeletons in molecular crystals are mostly independent of the size and shape of the molecule, while the total packing can vary in proteins compared to small molecules.
Topology versus thermodynamics
To find correlations between thermodynamic parameters and molecular connectivity we examined all the 394 protein structures of the monomer set. A total of 3,156 crystal packing contacts were identified. We observed that the number of intermolecular contacts (NC) is strictly related to the ÎG diss . PISA does not compute the number of intermolecular contacts but computes the buried solvent area (BSA; the patch of solvent accessible surface of molecules A and B that become inaccessible to the solvent because of the interaction between A and B). There is a strict relationship between NC and BSA:
with Pearsonâs correlation coefficient of 0.945. This is obviously not surprising and it means that a single intermolecular contact buries nearly 7 à 2 of solvent accessible surface. This means that about 30% of the surface of an atom of radius equal to 1.4âà is buriable in a protein-protein interface, while the rest is buried intra-molecularly. Also the relationship between BSA and ÎG diss is strict:
with Pearsonâs correlation coefficient of 0.793. We note that this is a mere mathematical fit of the data, which has no physico-chemical meaning. Consequently, it is possible to fit the dependence of ÎG diss on NC as:
with Pearsonâs correlation coefficient of 0.770. The average absolute difference between the ÎG diss values computed by PISA and those computed on the basis of the NC values is only 2.00 (standard error 0.03) kcal/mol. We observe that the |ÎG diss - NC | values tend to increase slightly with NC. The average |ÎG diss - NC | values are only 1.46 (0.03) kcal/mol when 1ââ¤âNCââ¤â30 (808 observations) and they growth to 3.6(0.3) kcal/mol when 210ââ¤âNCââ¤â240 (74 observations). For larger NC values, there are too few observations to make reasonable statistics. Although this trend indicates that the fitting model is statistically inappropriate, for practical reasons it is sufficient to estimate ÎG diss values with an average error of about 2â3âkcal/mol in most cases.
Conclusions
The topologies of carefully selected high-quality protein crystal structures have been determined and compared to those of a large set of small molecule crystal structures. It is known that protein crystals are suboptimally packed, since they contain channels filled with liquid water. Therefore, the coordination numbers in protein structures are considerably lower than in small molecule crystals. Surprisingly, however, we have found their topologies very similar. Although it is impossible, based on the data presented in this communication, to provide a strict and definitive explanation for this surprising similarity, it is clear that some molecular mechanism must exist. It is possible to hypothesize, for example, that during the early stages of nucleation, only relatively few types of assemblies, with a specific stereochemistry, are sufficiently stable to be able to accumulate and to continue to grow into real crystals. At this regard, it might be interesting to remember that three main differences exist between protein and small molecule crystallization41,42,43,44. First, the kinetics of protein crystal nucleation and growth are, in general, two or three orders of magnitude slower than for small molecules, because of the protein larger size and lowered diffusivity45,46. Second, protein crystals nucleation occurs at very high levels of supersaturation, often two or three orders of magnitude greater than that required to sustain crystal growth. Third, proteins may assume several distinctive solid states that include amorphous precipitates, oils or gels as well as crystals, and most of these are kinetically favored. Further studies are necessary and additional data must be considered to find out the rationale why topology seems to be independent of packing efficiency and crystallization.
References
McPherson, A. A brief history of protein crystal growth. J. Cryst. Growth 110, 1â10 (1991).
Hunefeld, F. L. Die Chemismus in der thierischen Organisation. (Brockhaus, 1840).
Kendrew, J. C. & Perutz, M. F. X-ray studies of compounds of biological interest. Ann. Rev. Biochem. 26, 327â372 (1957).
Kendrew, J. C., Bodo, G., Dintzis, H. M. P. R. G., Wyckoff, H. & Phillips, D. C. A three-dimensional model of the myoglobin molecule obtained by x-ray analysis. Nature 181, 662â666 (1958).
Muirhead, H. & Perutz, M. F. Structure of haemoglobin. A three-dimensional Fourier synthesis of reduced human haemoglobin at 5.5 Ã resolution. Nature 199, 633â638 (1963).
Bernstein, F. C. et al. The Protein Data Bank: a computer-based archival file for macromolecular structures. J Mol Biol 112, 535â542 (1977).
Berman, H. M. et al. The Protein Data Bank. Nucleic Acids Res 28, 235â242 (2000).
Doye, J. P. K., Louis, A. A. & Vendruscolo, M. Inhibition of protein crystallization by evolutionary negative design. Phys. Biol. 1, 9â13 (2004).
Koopmann, R. et al. In vivo protein crystallization opens new routes in structural biology. Nat Meth 9, 259â262 (2012).
Janin, J. & Rodier, F. Protein-protein interaction at crystal contacts. Proteins 23, 580â587 (1995).
Carugo, O. & Argos, P. Protein-protein crystal-packing contacts. Protein Sci. 6, 2261â2263 (1997).
Janin, J. Specific versus non-specific contacts in protein crystals. Nat. Struct. Biol. 4, 973â974 (1997).
Carugo, O. & DjinoviÄ-Carugo, K. How many packing contacts are observed in protein crystals? J. Struct. Biol. 180, 96â100 (2012).
Carugo, O. & Djinovic-Carugo, K. Packing bridges in protein crystal structures. J. Appl. Cryst. 47, 458â461 (2014).
Djinovic-Carugo, K. & Carugo, O. Structural biology of the lanthanides - Mining rare earths in the Protein Data Bank. J. Inorg. Biochem. 143, 69â76 (2015).
Matthews, B. W. Solvent content of protein crystals. J. Mol. Biol. 33, 491â497 (1968).
Andersson, K. M. & Hovmoeller, S. The protein centent in crystals and packing coefficients in different space groups. Acta Cryst. D56, 789â790 (2000).
Kantardjieff, K. A. & Rupp, B. Matthews coefficient probabilities: Improved estimates for unit cell content of protein, DNA, and protein-nuclei acid complex crystals. Prot. Sci. 12, 1865â1871 (2003).
Quiocho, F. A., McMurray, C. H. & Lipscomb, W. N. Similarities between the conformation of arsanilazotyrosine 248 of carboxypeptidase A in the crystalline state and in solution. Proc. Natl. Acad. Sci. USA 69, 2850â2854 (1972).
Acheson, J. F., Bailey, L. J., Brunold, T. C. & Fox, B. G. In-crystal reaction cycle of a toluene-bound diiron hydrolase. Nature 544, 191â195 (2017).
Kooystra, P. J., Kalk, K. H. & Hol, W. G. Soaking in Cs2SO4 reveals a caesium-aromatic interaction in bovine-liver rhodanese. Eur. J. Biochem. 177, 345â349 (1988).
Carugo, O. Structure and function of water molecules buried in the protein core. Curr Protein Pept Sci 16, 259â265 (2015).
Carugo, O. Statistical survey of the buried waters in the Protein Data Bank. Amino Acids 48, 193â202 (2016).
Carugo, O. Protein hydration: Investigation of globular protein crystal structures. Int. J. Biol. Macromol. 99, 160â165 (2017).
Drenth, J. & Haas, C. Protein crystals and their stability. J. Cryst. Growth 122, 107â109 (1992).
Kitaigorodskii, A. I. Molecular Crystals and Molecules. (Academic Press, 1973).
Peresypkina, E. V. & Blatov, V. A. Molecular coordination numbers in crystal structures of organic compounds. Acta Cryst. B56, 501â511 (2000).
Peresypkina, E. V. & Blatov, V. A. Topology of molecular packings in organic crystals. Acta Cryst. B56, 1035â1045 (2000).
Djinovic-Carugo, K. & Carugo, O. Criteria to extract high quality Protein Data Bank subsets for structure users. Metods Mol. Biol. in the press (2016).
van Beusekom, B., Perrakis, A. & Joosten, R. P. Data mining of macromolecular structures. Methods. Mol. Biol. 1415, 107â138 (2016).
Djinovic Carugo, K. & Carugo, O. Missing strings of residues in protein crystal structures. Intr. Dis. Prot. 3, 1â7 (2015).
van Beusekom, B., Perrakis, A. & Joosten, R. P. Data mining of macromolecular structures. Methods. Mol. Biol. (in the press).
Doerr, A. Membrane protein structures. Nature Meth. 6, 35 (2009).
Blatov, V. A. A method for topological analysis of rod packings. Struct. Chem. 27, 1605â1611 (2016).
Blatov, V. A., Shevchenko, A. P. & Proserpio, D. M. Applied topological analysis of crystal structures with the program package ToposPro. Cryst. Growth Des. 14, 3576â3586 (2014).
Alexandrov, E. V., Blatov, V. A., Kochetkov, A. V. & Proserpio, D. M. Underlying nets in three-periodic coordination polymers: topology, taxonomy and prediction from a computer-aided analysis of the Cambridge Structural Database. CrystEngComm 13, 3947â3958 (2011).
OâKeeffe, M., Peskov, M. A., Ramsden, S. J. & Yaghi, O. M. The Reticular Chemistry Structure Resource (RCSR) Database of, and Symbols for, Crystal Nets. Acc. Chem. Res. 41, 1782â1789 (2008).
Blatov, V. A. & Proserpio, D. M. Topological relations between three-periodic nets. II. Binodal nets. Acta Cryst. A65, 202â212 (2009).
Aman, F. et al. Multilevel topological description of molecular packings in 1,2-benzothiazines. CrystEngComm 16, 1963â1970 (2014).
Krissinel, E. & Henrick, K. Inference of macromolecular assemblies from crystalline state. J. Mol. Biol. 372, 774â797 (2007).
Kuznetsov, Y. G., Malkin, A., Greenwood, A. & McPherson, A. Interferometric studies of growth kinetics and surface morphology in macromolecular crystal growth: canavalin, thaumatin, and turnip yellow mosaic virus. J. Struct. Biol. 114, 184â196 (1995).
Malkin, A. J., Kuznetsov, Y. G., Glantz, W. & McPherson, A. Atomic force microscopy studies of surface morphology and growth kinetics in thaumatin crystallization. J. Phys. Chem. 100, 11736â11743 (1995).
Malkin, A. J., Kuznetsov, Y. G. & McPherson, A. An in situ AFM investigation of catalase crystallization. Surf. Sci. 393, 95â107 (1997).
McPherson, A. & Gavira, J. A. Introduction to protein crystallization. Acta Cryst. F70, 2â20 (2014).
Chernov, A. A. Protein crystals and their growth. J. Struct. Biol. 142, 3â21 (2003).
Velikov, P. G. & Chernov, A. A. The physics of protein crystallization. Solid State Phys. 57, 2â147 (2002).
Pourayoubi, M., Keikhaa, M. & Parvezb, M. N, Nâ²-Bis(2-methylphenyl)-Nâ²â²-(2,2,2-trichloroacetyl) phosphoric triamide. Acta Cryst. E67, o2792 (2011).
Paesen, G. C. et al. A Tick Protein with a Modified Kunitz Fold Inhibits Human Tryptase. J. Mol. Biol. 368, 1172â1186 (2007).
Shen, Y. B., Wang, M., Liang, Q. K. & Luo, J. M. 11α,15 α -Dihydroxyandrost-4-ene-3,17. Acta Cryst. E67, o2752 (2011).
Olsen, J. G., Flensburg, C., Olsen, O., Bricogne, G. & Henriksen, A. Solving the Structure of the Bubble Protein Using the Anomalous Sulfur Signal from Single-Crystal in-House Cu Kalpha Diffraction Data Only. Acta Cryst. D60, 250â255 (2004).
Sarma, B. K., Yousufuddinb, M. & Kodadek, T. Acyl hydrazides as peptoid sub-monomers. Chem. Coomun. 47, 10590â10592 (2011).
Ielasi, F. S., Decanniere, K. & Willaert, R. G. The epithelial adhesin 1 (Epa1P) from the human-pathogenic yeast Candida Glabrata: Structural and functional study of the carbohydrate-binding domain. Acta Cryst. D68, 210â217 (2012).
Adams, C. J. et al. The P53 cofactor strap exhibits an unexpected Tpr motif and oligonucleotide-binding (Ob)-fold structure. Proc. Natl. Acad. Sci. USA 109, 3778â3783 (2012).
Acknowledgements
We thank E. Krissinel for his advices on the use of PISA and on the interpretation of the results. O.C. thanks Kristina DjinoviÄ for her kind hospitality in Vienna. O.A.B., E.O.M., V.A.B., D.M.P. are grateful to the Russian Government (Grant 14.B25.31.0005) for support.
Author information
Authors and Affiliations
Contributions
O.C., V.A.B. and D.M.P. conceived and planned the project; O.C., V.A.B., O.A.B. and E.O.M. carried out the computations; O.C., V.A.B. and D.M.P. wrote the manuscript. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Carugo, O., Blatova, O.A., Medrish, E.O. et al. Packing topology in crystals of proteins and small molecules: a comparison. Sci Rep 7, 13209 (2017). https://doi.org/10.1038/s41598-017-12699-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-12699-4
This article is cited by
-
Simplify to understand: how to elucidate crystal structures?
Structural Chemistry (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.