Crystallographic characterization of two novel crystal forms of human insulin induced by chaotropic agents and a shift in pH
© Norrman and Schluckebier; licensee BioMed Central Ltd. 2007
Received: 21 August 2007
Accepted: 19 December 2007
Published: 19 December 2007
Insulin is a therapeutic protein that is widely used for the treatment of diabetes. Its biological function was discovered more than 80 years ago and it has since then been characterized extensively. Crystallization of the insulin molecule has always been a key activity since the protein is often administered by subcutaneous injections of crystalline insulin formulations. Over the years, insulin has been crystallized and characterized in a number of crystal systems.
Interestingly, we have now discovered two new crystal forms of human insulin. The crystals were obtained when the two chaotropic agents, urea and thiocyanate were present in the crystallization experiments, and their structures were determined by X-ray crystallography. The crystals belong to the orthorhombic and monoclinic crystal systems, with space groups C2221 and C2 respectively. The orthorhombic crystals were obtained at pH 6.5 and contained three insulin hexamers in R6 conformation in the asymmetric unit whilst the monoclinic C2 crystals were obtained at pH 7.0 and contained one R6 hexamer in the asymmetric unit. Common for the two new crystals is a hexamer-hexamer interaction that has not been found in any of the previous crystal forms of insulin. The contacts involve a tight glutamate-glutamate interaction with a distance of 2.3 Å between groups. The short distance suggests a low barrier hydrogen bond. In addition, two tyrosine-tyrosine interactions occupying a known phenol binding pocket contribute to the stabilization of the contacts. Within the crystals, distinct binding sites for urea were found, adding further to the discussion on the role of urea in protein denaturation.
The change in space group from C2221 to C2 was primarily caused by an increase in pH. The fewer number of hexamer-hexamer interactions comprising the short hydrogen bond in the C2 space group suggest that pH is the driving force. In addition, the distance between the two glutamates increases from 2.32 Å in the C2221 crystals to 2.4 Å in the C2 crystals. However, in both cases the low barrier hydrogen bond and the tyrosine-tyrosine interaction should contribute to the stability of the crystals which is crucial when used in pharmaceutical formulations.
The therapeutic hormone insulin is a small protein used daily in the medical treatment of diabetes by millions of people. The primary administration route is by subcutaneous injections of microcrystals or mixtures of microcrystals and amorphous protein. After subcutaneous injection, the insulin crystals dissolve slowly, thus leading to a slow intermediate release of insulin into the blood stream. Ever since the biological function of insulin was discovered in the 1920s , the molecule has been widely characterized both biophysically and structurally. The crystallographic structure was one of the first protein structures determined . It has since then been crystallized in a number of space groups of which the most common belong to the monoclinic, rhombohedral, cubic and tetragonal crystal forms. The type, size and morphology of the crystals affect how fast insulin is released, which is why crystallization of insulin has been studied extensively. Alternative administration routes are presently a rapidly expanding research field and insulin microcrystals may be well suited for other delivery methods, including pulmonary delivery or sustained release formulations [3–5]. The insulin molecule consists of two chains, A and B, with 21 and 30 residues respectively. Chain A is built up by two helical fragments separated by a short loop linked to one of the helices by an intra-chain disulfide bond. Two additional disulfide bonds link chain A to the larger chain B. In the biologically active form, insulin exists as a monomer in which chain B contains a central helical region flanked by two elongated parts. In the presence of divalent ions like zinc, the monomers assemble into hexamers , where each of the two central zinc ions is coordinated by three histidine residues. In the hexameric form, chain B has been shown to exist in two allosteric states denoted T and R . The R state has two allosteric binding sites usually referred to as the phenolic binding site and HisB10 anion site. The T → R state transition and the two different B chain conformations, referred to as T6 and R6, have been described by spectroscopic and crystallographic studies. The T6 conformation, which is characterized by an extended conformation of residues 1–8 of chain B, is obtained at low chloride concentrations and in absence of phenol derivatives [8, 9]. Phenolic derivatives are used as preservative in insulin pharmaceutical formulations. The most commonly used are phenol, meta-cresol, resorcinol and methylparaben. The R conformation is obtained in presence of these derivatives and at high chloride concentrations. In this form the first eight residues of the B chain adopt a helical conformation, which together with the central helical segment gives a continuous helix which includes residues B1 to B19 [10–13]. This transition from an extended to an alpha-helical conformation causes the first residue of chain B, PheB1 to undergo a ~30 Å shift in position. Although chloride is the most commonly used anion, other anions such as SCN-, OCN-, CN-, N3- and NO2- have also been shown to be useful [7, 14, 15]. Like chloride, in the absence of phenolic derivatives and at high concentrations, these anions are able to induce the R state in three of six monomers in a hexamer. The remaining three monomers have an extended conformation (T state) in the region including residues B1 to B8. The R state of the first three monomers is incomplete with residues B4 to B8 in a helical conformation, while residues B1 to B3 have an extended conformation. This hexamer configuration is denoted T3R3f, where the 'f' indicates a frayed R conformation [16, 17].
We here present a study which shows that certain chaotropic additives can induce two novel types of insulin crystals, and that the type of crystals formed depends on the charge state of insulin, i.e. is pH dependent. The structures and crystal packing interaction of the two new crystal forms have been analyzed and compared to crystal packing interactions in other previously known insulin crystals.
Results and discussion
Crystallization of insulin is of high importance in pharmaceutical formulations and in insulin manufacturing and has been systematically investigated since the 1920s [6, 18–20]. By introducing chaotropic agents in the crystallization experiments, we succeeded in identifying two new crystal forms of native human insulin. The crystals were found using two different crystallization experiments. The first crystals were obtained in a crystallization screen with varying concentrations of urea and sodium chloride in presence of zinc and resorcinol. The crystals were initially characterized by X-ray powder diffraction and were shown to have a powder pattern differing from previously known insulin crystal forms . Further optimizations of crystallization conditions resulted in crystals suitable for single crystal analysis and determination of crystal system which was found to be orthorhombic in space group C2221. The crystals appeared in the pH range 6.0 – 6.5 while a second type of crystals, characterized as monoclinic with space group C2 grew in the pH range 6.5 – 7.0. In the overlapping pH interval around pH 6.5, the C2221 crystals were present in wells with lower salt and urea concentrations. A few drops contained a mixture of the two crystal types. In a parallel experiment, the urea and sodium chloride were substituted for thiocyanate. Interestingly, the same two crystal types appeared here, at the same pH intervals, with a clearer pH distinction at pH 7.0. Crystallization with thiocyanate or chloride ions without a phenol derivative has previously been shown to stabilize the T3R3f form of hexameric insulin in a rhombohedral crystal system [16, 17]. In our case when resorcinol and thiocyanate were present, the orthorhombic C2221 and monoclinic C2 crystals appeared. The well characterized monoclinic crystals in space group P21 [10, 11] were present at pH values above 7.0 in wells with low salt and high urea concentration and increased in frequency as the pH was raised to become the dominating crystal form at pH ≥ 7.5. The crystals obtained in presence of urea will be referred to as C2221urea and C2urea while the two crystal forms obtained with thiocyanate are referred to as C2221scn and C2scn.
Structure of insulin in the orthorhombic lattice
Structure of insulin in the monoclinic lattice
The crystals obtained at slightly higher pH (pH 7.0) belong to the monoclinic space group C2 with cell dimensions a = 100 Å, b = 60 Å, c = 62 Å, β = 116°. They contain one hexamer with R6 conformation in the asymmetric unit and have a solvent content of 50%. The crystal packing is shown in Figure 1c. Both the C2urea and C2scn structures have two zinc atoms/hexamer located 14.7 Å and 15.2 Å apart respectively. The zinc coordination is identical to the C2221 crystals. In the C2urea structure, two additional resorcinol molecules could be fitted into the electron density. The location of the first is very close to binding site II, described in . At this site, the first hydroxy group of the resorcinol molecule is hydrogen bonded to the OG atom of SerB9 (2.9 Å) in an alternating conformation. The other hydroxyl group interacts with the carbonyl oxygen of GluB13 (2.5 Å) and a water molecule (3.1 Å). The water molecule, in turn, makes a hydrogen bond to the carbonyl oxygen of SerB9 (3.1 Å). In contrast to the phenolic binding interactions observed earlier in the PDB entry 1ZEG, where the phenolic oxygen hydrogen binds to HisB5, the angular orientation of the HisB5 in the C2urea structure does not seem to permit any interaction with the resorcinol molecule. The second additional resorcinol is located at the surface of the insulin in a solvent channel between two monomers, surrounded by water molecules. At the end of refinement, one glycerol molecule from the cryo-solution was added to both C2 structures. The glycerol molecule in the C2scn structure was found at a corresponding position as in the C2221scn structure, while the glycerol molecule in the C2urea structure was found in the solvent channel leading towards one of the zinc atoms, where it interacted with surrounding water molecules.
Data processing and refinement statistics for the orthorhombic and monoclinic insulin crystallized in presence of urea and thiocyanate respectively.
Urea containing crystals
Na-SCN containing crystals
Values in parenthesis are for the highest resolution shell
Cell axis a, b, c (Å)
58.9, 219.3, 223.7
100.2, 60.2, 62.9
59.0, 219.5, 224.5
100.6, 60.8, 62.1
Cell angles α β γ (deg)
90.0 90.0 90.0
90.0 116.2 90.0
90.0 90.0 90.0
90.0 116.1 90.0
Diffraction limit (Å)
Highest resolution shell (Å)
2.1 – 2.05
1.75 – 1.70
2.0 – 1.97
1.75 – 1.704
No. of observations
450 276 (23593)
72 333 (3616)
885 954 (9682)
143 773 (4037)
91 251 (6331)
34 195 (1952)
102 732 (3997)
34 924 (1890)
Rmerge a (%)
Values in parenthesis are for the highest resolution shell
Resolution range (Å)
28.31 – 2.05
19.57 – 1.70
19.99 – 1.97
19.78 – 1.70
Highest resolution shell (Å)
2.10 – 2.05
1.74 – 1.70
2.02 – 1.97
1.75 – 1.704
No. of reflections
86 749 (6376)
32 471 (1621)
97 508 (6629)
33 178 (1700)
R value (%) b
Rfree value (%) b
Bond length (Å)
Bond angles (deg)
B factors (Å 2) d
Average all atoms
All four insulin molecules are structurally very similar. Pair-wise superposition and comparison of the C2 structures results in a root-mean-square (r.m.s.) distance between corresponding Cα atoms of 0.35 Å and of 0.88 Å when all common atoms are included. For the C2221 structures the same r.m.s. distances are 0.26 Å and 0.58 Å.
In contrast to the T3R3f conformation, where B1–B3 have an extended conformation, only resides B1–B2 have a non-helical state. A similar, non-helical conformation of the PheB1 residue has previously been observed for one of the B-chains in an R6 insulin in complex with resorcinol (PDB ID: 1EVR) . In that case, the carbonyl oxygen of PheB1 is coordinating a sodium ion which was further coordinated by the C terminal AsnA21 of a symmetry-related molecule. In our case, PheB1 is stabilized in a non-helical conformation by a hydrogen bond between amide nitrogen of PheB1 to the carbonyl oxygen of ThrA8 in a neighboring molecule, or a hydrogen bond between the carbonyl oxygen of PheB1 and the amide nitrogen of AsnB3 in the same chain. There are further interactions with symmetry-related molecules, such as the PheB1 amide nitrogen interactions with the OH group of a symmetry-related TyrA14, or the carbonyl oxygen of CysA20 and AsnA18.
In close proximity to the PheB1 residue of the three B chains with non-helical conformation in the C2urea structure there was an electron density peak with a height of 5 σ in the 2Fo-Fc map and ~5 σ in a Fo-Fc difference map. The location of the peak was close to the position where the carbonyl oxygen of PheB1 would be located if the conformation was a-helical. Given the observed electron density, a chloride ion was fitted into this position. It is coordinated to the amide nitrogen of HisB5 (3.2 Å) and two or three water molecules at an average distance of 3.3 Å. The corresponding sites in the C2221urea structure were too disordered to be interpreted in a similar manner.
Location of main-chain and side-chain atoms was ambiguous for the residues LysB29 and ThrB30 in most of the chains in the four structures. Furthermore, the following residues were modeled with alternating side chain conformations; C2urea structure: GlnB4.2, SerB9.5, ValB18.5, LeuB17.6; C2scn structure: GluB13.3, ValB18.4; C2221urea structure: LeuB17.I.1, ArgB22.I.4, GlnB4.I.6, ValB18.II.4 AsnB3.III.1; C2221scn structure: ValB18.I.2, ValB18.II.4, ValB18.II.5 (the roman numerals refer to the hexamer number while the single integer following a punctuation indicates monomer).
A special crystal interaction at the dimer-dimer interface
As the pH is increased to 7.0, the second crystal form C2 appears. In this crystal form, there is only one crystal packing interaction comprising the tight glutamate-glutamate and tyrosine-tyrosine interaction, Figure 3. The increased pH could be the reason for the smaller number of such contacts. At higher pH, the shared hydrogen between the two glutamates becomes more delocalized and the repulsive forces will dominate. Consequently, the distance between the carboxylates is longer, 2.40 Å, versus 2.32 Å for the C2221 structures, indicating a weaker interaction at this pH. At pH values above 7.5, only the monoclinic P21 crystal form  is observed, in which no such interface exists.
Analysis of crystal contact surfaces
Properties of crystal packing contact surfaces of insulin hexamers in seven different space groups.
P4 3 2 1 2
Solvent content. (%)
protein atoms (Å2)
exposed atoms (Å2)
contact atoms (Å2)
Tot SASA (Å2) a
fract of tot SASA (%)
Atom type specific area as a fraction of A cont
C area (%)
O area (%)
N area (%)
S area (%)
C aromatic (%)
C aliphatic (%)
Pos charged Area (%)
Neg charged Area (%)
The contact surfaces were characterized as either polar (oxygen and nitrogen atoms, including ionisable groups) or non-polar (carbons). The four structures presented in this study constitute a group with a high fraction of non-polar contact surface, ranging from 53% to 56% of the total contact area, compared to 41% to 50% for the other crystal forms. The monoclinic P21 crystal form is the most hydrophilic, with a 40/60 distribution between hydrophobic and hydrophilic contact area. This analysis is limited in that bound water molecules were not considered in the crystal contact interactions since the criteria for modeling water molecules may vary among crystallographers and are also dependent on data quality. Several interactions could however involve hydrogen bonds to water molecules. Side chains with missing atoms were rebuilt in order to use the surface with an atom composition representing the true surface for the property calculations. They were however rebuilt automatically and could potentially be in a wrong orientation.
Comparing the residue identity of the crystal contacts for the crystals presented in Table 2 shows that seven of the interface residues are common for all crystal forms (GlnA5, ThrA8, TyrA14, GlnA15, AsnA18, TyrA19 and PheB1). Altogether, the contact sites for the six crystal types compared in this study cover almost the entire surface of an insulin hexamer. A comparison of the exposed residues with the residues involved in crystal contacts shows that all residues with an exposure of more than 20% participate in some contact interface. The degree of exposure was calculated according to . A number of studies, where crystal packing contacts have been systematically investigated [26, 27] conclude that atomic composition within crystal contacts is indistinguishable from that of the protein surface and is rather non-specific. Studies of pancreatic ribonuclease  and cutinase , crystallized in a number of space groups, showed in accordance with the present study that virtually the entirely protein surface can be involved in crystal contacts.
Insulin has been shown to be tolerant of high concentrations of urea and other denaturants [30, 31] and urea has previously been used to increase its solubility. One example is the inclusion of urea to promote growth of larger crystals of an insulin-protamine complex . No significant conformational changes were detected in our structures as a result of direct urea interactions. In spite of its common use as a chemical denaturant of proteins, the molecular mechanism of urea-mediated unfolding is not known. Identification of denaturant interactions with proteins may give insight into the early stages of protein unfolding . For a denaturant to be effective, protein-solvent interactions must be disturbed and this is thought to happen either through a direct or an indirect mechanism . A direct mechanism would involve binding of urea molecules to the protein surface and thus compete with water-protein interactions and enhance the solubility of hydrophobic residues. Indirect urea denaturation would involve disruption of solvent-mediated hydrophobic interactions which would destabilize the protein structure. In addition, studies have shown that urea and guanidine hydrochloride at sub-denaturant concentrations stabilizes proteins at a sub global level in a mechanism called protein stiffening [35, 36]. The present study shows that urea at concentrations ~3 M has one specific binding site on the surface of the insulin molecule, interacting with backbone carbonyl groups of primarily GlnA5 but also of SerA9 and IleA10 residues. Given the high concentrations of urea present in the crystallization experiments, we would expect to detect even weak binding sites with a Kd of several hundred mM. Thus it seems unlikely that insulin denaturation occurs via a direct mechanism which requires binding of several urea molecules. On the other hand we see no signs of partial unfolding in our structures which would be indicative of an indirect mechanism. A recent study suggests that the denaturant effect of urea is neither due to a direct or indirect mechanism but rather an effect of a reduction of ion pairing between ionic and polar groups at aggregate surfaces , something which also could explain the relative higher fraction of non-polar surface at crystal contacts in our crystal forms.
In spite of the extensive research on insulin crystallization during the last 80 years we could surprisingly identify two new crystal forms of human insulin. The major factor behind the change of space group from C2221 to C2 was an increase in pH from ~6.5 to 7.0. A comparison of the crystal contacts in the two space groups showed that the interfaces are very similar. The most obvious difference and the most remarkable contact interaction was the double TyrA14-TyrA14 interaction combined with a tight GluA17-GluA17 interaction. When taking all symmetry-related contacts into account it was found that this contact type was more frequent in the C2221 crystals when compared to the C2 crystals. This crystal packing interaction has not been observed in any of the previously known insulin crystals. The addition of chaotropes such as urea and thiocyanate to the crystallization buffer can have an effect on the protein charge properties by charge screening. This could explain altered pKa values of the glutamic acids involved in the crystal contacts and the higher fraction of hydrophobic crystal contacts in the present crystal forms compared to previously known insulin crystals. The short carboxyl-carboxyl interaction indicates the presence of a shared proton between the two groups and would be a strong low barrier hydrogen bond which should contribute to the stability of the crystals.
In the case of insulin, much effort has been put into the modification of the dimerization interactions [38, 39] and hexamer formations [40, 41], but less focus has been on the inter-hexamer interactions in solution or within crystals. Modifications of surface residues can induce changes in crystal packing due to breaking of existing interactions and/or formation of new ones [42, 43]. Engineering of the protein surface to specifically induce a change in the crystal form or improve stability in lattice contacts may produce a better diffracting crystal . Since insulin is a therapeutic protein administered also in crystalline form, the discovery and analysis of new polymorphic forms has implications beyond providing improved crystals for structural studies. The inter-hexamer interaction found in the present structures provides an interesting and novel interface that is specific for these two crystal forms. A single additional hydrophobic or several polar interactions may increase the half-life of a protein by several orders of magnitude . Thus, additional inter-hexamer interactions can increase the stability and thus the shelf life of crystalline insulin formulations. The structures presented here provide a framework for further site-directed mutagenesis studies of the residues involved in inter-hexamer interactions, aimed at providing improved formulations useable within the rapidly advancing field of alternate delivery routes of crystalline biopharmaceuticals .
This study also demonstrates the usefulness of X-ray powder diffraction (XRPD) on protein samples. The small size of the initially obtained micro crystals made visual analysis and single crystal X-ray diffraction difficult. However, the combination of XRPD and principal component analysis (PCA) facilitated the identification of a new crystal form . Since the use of proteins as therapeutic agents is a growing field, applications of protein XRPD, similar to the present study, will have an important role during discovery and development of therapeutic protein formulations.
Human insulin was obtained from Novo Nordisk A/S (Denmark). Crystals were grown by hanging-drop vapor diffusion technique at 291 K. The crystals were obtained in crystallization experiments with varying reservoir concentrations of NaCl and urea. A protein solution containing 6.9 mg/ml of human insulin, with zinc content corresponding to two zinc ions per hexamer and 50 mM resorcinol were mixed with equal volumes of reservoir solution. The protein solution was filtered through a 0.22-micrometer centrifugal filter (Ultrafree-MC, Millipore, USA) prior to crystallization. The best diffracting crystals were obtained from the following conditions. For the C2221 crystals: 2 M NaCl, 3 M urea, 100 mM phosphate buffer pH 6.5; for the C2 crystals: 2.5 M NaCl, 4 M urea, 100 mM phosphate buffer pH 7.0. The same two crystal forms were obtained when including 15 mM NaSCN in the protein solution in absence of urea and using a reservoir solution containing 5% (v/v) ethanol and 200 mM phosphate buffer at pH 6.5 and 7.0, respectively. Crystals with dimensions of about 0.25 mm on one edge were detected after two days. The two orthorhombic crystals in space group C2221 diffracted to a resolution of 2.0 Å while data from the crystals in space group C2 could be collected to 1.7 Å.
Data collection and refinement
Data sets were collected from a single crystal of each type at 100 K using synchrotron radiation (Maxlab synchrotron, Lund, Sweden, beamline 911-2 and 911-3 ), with a MarMosaic 225 CCD detector (MarResearch, Evanston, USA). The urea containing crystals were soaked in a cryo-solution containing 23% glycerol and 77% reservoir prior to freezing in liquid nitrogen. For the NaSCN crystals a cryo-solution containing 30% glycerol was needed. All data sets were processed and scaled using the XDS package .
For the C2221 crystals with urea, an additional low-resolution data set was collected from the same crystal and merged with the high-resolution data. An in house structure of a hexamer with R6 conformation, excluding all non-protein atoms except zinc, was used as search model for molecular replacement in Molrep . Three hexamers were found in the asymmetric unit, corresponding to a solvent content of 64%. During refinement in Refmac , 5% of the data was excluded and used for calculation of the R-free value. The initial 2Fo-Fc map clearly indicated positions of chloride ions and resorcinol molecules (six in each hexamer). After several rounds of refinement using the maximum likelihood option in Refmac and manual adjustments of main-chain and side-chain atoms in Coot , TLS refinement  was employed, with each monomer defined as a separate TLS group. Subsequently, water was added by the find-water function in Coot. At the end of the refinement, urea molecules were modeled into the electron density using positive Fo-Fc peaks, where the shape of the 2Fo-Fc density was flat and reminiscent of the triangular shape of a urea molecule. 14 such positions were found with an average B-factor of 53.1 Å 2. The final number of water molecules was 634 with an average B-factor of 41.7 Å 2.
The C2 structure obtained with urea was solved with the same search model as the C2221 data set and the procedure for structure solving and refinement followed the same route. In the beginning of the refinement, there was a clear 2Fo-Fc density for an extra resorcinol binding site, and at later stages of refinement, a second additional position showed density resembling a resorcinol molecule. Seven urea molecules (average B-factor = 39.3 Å 2) and 257 water molecules were built into the density at the end of refinement.
For the two structures co-crystallized with NaSCN, one high and one low resolution data set were collected and subsequently merged. The structures were solved using the urea containing structures as search models in molecular replacement rounds (excluding non-protein atoms except zinc). Based on the experience of flexible residues in the B-chain terminals, the search models were truncated at both ends to reduce bias (PheB1, ValB2, LysB29 and ThrB30). Refinement followed the same scheme as for the urea crystals. In total, six resorcinol molecules were modeled into the C2 structure and 18 in the C2221 structure (six in each hexamer). 313 water molecules were fitted into the C2 structure and 755 into the C2221. Based on the shape of the electron density, each zinc atom was found to interact with one thiocyanate molecule. Data collection details and refinement statistics for all four structures are summarized in Table 1. In this paper, the crystals obtained in presence of urea will be referred to as C2221urea and C2urea whilst the two forms obtained with thiocyanate are referred to as C2221scn and C2scn.
Analysis of crystal packing
Symmetry-related molecules were generated in Pymol  for analysis of crystal contacts. A rather strict criterion for identification of symmetry contacts was used by searching for symmetry-related atoms within 4 Å from the reference protein. In order to limit the analysis to protein-protein interactions, non-protein atoms such as water and urea molecules were removed prior to searches. Residues with missing side-chain atoms were reconstructed in Swiss-PdbViewer  using the 'auto reconstruct residues with missing atom' function. When dual conformations of residue side chains were present, the one with highest occupancy, or when equal, the conformation denoted with an 'A' in the pdb file was chosen. The contact area (Acont) between the reference molecule and neighboring symmetry-related molecules was defined as the solvent accessible area buried by symmetry-related molecules. The solvent accessible surface area (SASA) was calculated in areaImol (within the CCP4 program package ) that utilizes the algorithm of Lee & Richards . A water probe of 1.4 Å was used. Acont for each atom was obtained by taking the SASA for the reference molecule alone, minus the SASA when contact atoms were present. All atoms with reduced SASA were assumed to be involved in contacts with a symmetry molecule. Basic physicochemical properties related to hydrophobic surface area were calculated for the atoms of the reference molecule involved in contacts. Contact areas were split into hydrophobic (carbons) and hydrophilic (nitrogen and oxygens), aliphatic-, aromatic carbons and positive-, negative charged area. For the sake of simplicity and ease of comparison, crystals without a full hexamer in the asymmetric unit were complemented with the required symmetry molecules by applying appropriate symmetry operations to generate a hexamer prior to calculations. A hexamer was considered as one molecule.
Coordinates with structure factors have been deposited to the Protein Data Bank (PDB)  with the accession codes 2OLY, 2OLZ, 2OM0 and 2OM1.
We thank Lene Drube and Berit Bergerud Hansen for technical assistance and Salam Al-Karadaghi and Eva Johansson for valuable input during the preparation of this manuscript. The work was supported by the VTU (Ministry of Science, Technology and Innovation), Denmark and the Novo Nordisk CORA Training and research program.
- Banting FG, Best CH: Insulin its chemistry and physiology. J Lab Clin Med 1922, 7: 251–266.Google Scholar
- Adams MJ, Blundell TL, Dodson EJ, Dodson GG, Vijayan M, Baker EN, Harding MM, Hodgkin DC, Rimmer B, Sheat S: Structure of Rhombohedral 2 Zinc Insulin Crystals. Nature 1969, 224: 491–495. 10.1038/224491a0View ArticleGoogle Scholar
- Basu SK, Govardhan CP, Jung CW, Margolin AL: Protein crystals for the delivery of biopharmaceuticals. Expert Opin Biol Ther 2004, 4: 301–317. 10.1517/14712522.214.171.1241View ArticlePubMedGoogle Scholar
- DeFelippis MR, Dobbins MA: Insulin crystals for pulmonary administration. In Patent. Eli Lilly; 2001:PCT/US00/15026.Google Scholar
- Havelund S: Zinc free insulin crystals for use in pulmonary compositions. In Patent. Novo Nordisk A/S; 1998:PCT/DK98/00109.Google Scholar
- Schlichtkrull J: Chemical and biological studies on insulin crystals and insulin zinc suspensions. In Thesis. University of Copenhagen; 1958.Google Scholar
- Kaarsholm NC, Ko HC, Dunn MF: Comparison of Solution Structural Flexibility and Zinc Binding Domains for Insulin, Proinsulin, and Miniproinsulin. Biochemistry 1989, 28: 4427–4435. 10.1021/bi00436a046View ArticlePubMedGoogle Scholar
- Baker EN, Blundell TL, Cutfield JF, Cutfield SM, Dodson EJ, Dodson GG, Hodgkin DM, Hubbard RE, Isaacs NW, Reynolds CD: The structure of 2Zn pig insulin crystals at 1.5 A resolution. Philos Trans R Soc London Ser B 1988, 319: 369–456. 10.1098/rstb.1988.0058View ArticleGoogle Scholar
- Smith GD, Pangborn WA, Blessing RH: The structure of T6 human insulin at 1.0 A resolution. Acta Crystallogr D Biol Crystallogr 2003, 59(Pt 3):474–482. 10.1107/S0907444902023685View ArticlePubMedGoogle Scholar
- Derewenda U, Derewenda Z, Dodson EJ, Dodson GG, Reynolds CD, Smith GD, Sparks C, Swenson D: Phenol stabilizes more helix in a new symmetrical zinc insulin hexamer. Nature 1989, 338: 594–596. 10.1038/338594a0View ArticlePubMedGoogle Scholar
- Smith GD, Ciszak E, Magrum LA, Pangborn WA, Blessing RH: R6 hexameric insulin complexed with m-cresol or resorcinol. Acta Crystallogr D Biol Crystallogr 2000, 56(Pt 12):1541–1548. 10.1107/S0907444900012749View ArticlePubMedGoogle Scholar
- Smith GD, Dodson GG: Structure of a rhombohedral R6 insulin/phenol complex. Proteins 1992, 14: 401–408. 10.1002/prot.340140309View ArticlePubMedGoogle Scholar
- Whittingham JL, Edwards DJ, Antson AA, Clarkson JM, Dodson GG: Interactions of Phenol and m-Cresol in the Insulin Hexamer, and Their Effect on the Association Properties of B28 Pro Asp Insulin Analogues. Biochemistry 1998, 37: 11516–11523. 10.1021/bi980807sView ArticlePubMedGoogle Scholar
- Roy M, Brader ML, Lee RW, Kaarsholm NC, Hansen JF, Dunn MF: Spectroscopic signatures of the T to R conformational transition in the insulin hexamer. J Biol Chem 1989, 264: 19081–19085.PubMedGoogle Scholar
- Brader ML, Dunn MF: Insulin hexamers: new conformations and applications. Trends Biochem Sci 1991, 16: 341–345. 10.1016/0968-0004(91)90140-QView ArticlePubMedGoogle Scholar
- Whittingham JL, Chaudhuri S, Dodson EJ, Moody PCE, Dodson G: X-ray Crystallographic Studies on Hexameric Insulins in the Presence of Helix-S tabilizing Agents, Thiocyanate, Methylparaben, and Phenol. Biochemistry 1995, 34: 15553–15563. 10.1021/bi00047a022View ArticlePubMedGoogle Scholar
- Ciszak E, Smith GD: Crystallographic evidence for dual coordination around zinc in the T3R3 human insulin hexamer. Biochemistry 1994, 33: 1512–1517. 10.1021/bi00172a030View ArticlePubMedGoogle Scholar
- Abel JJ: Crystalline Insulin. Proc Natl Acad Sci U S A 1926, 12: 132–136. 10.1073/pnas.12.2.132PubMed CentralView ArticlePubMedGoogle Scholar
- Scott DA: Crystalline insulin. Biochem J 1934, 28: 1592–151602.PubMed CentralView ArticlePubMedGoogle Scholar
- Krayenbuhl C, Rosenberg T: Crystalline protamine insulin. Rep Steno Mem Hosp Nord Insulinlab 1946, 1: 60–73.Google Scholar
- Norrman M, Stahl K, Schluckebier G, Al Karadaghi S: Characterization of insulin microcrystals using powder diffraction and multivariate data analysis. J Appl Cryst 2006, 39: 391–400. 10.1107/S0021889806011058View ArticleGoogle Scholar
- Cleland WW, Frey PA, Gerlt JA: The Low Barrier Hydrogen Bond in Enzymatic Catalysis. J Biol Chem 1998, 273: 25529–25532. 10.1074/jbc.273.40.25529View ArticlePubMedGoogle Scholar
- Wohlfahrt G: Analysis of pH-dependent elements in proteins: Geometry and properties of pairs of hydrogen-bonded carboxylic acid side-chains. Proteins: Structure, Function, and Bioinformatics 2005, 58: 396–406. 10.1002/prot.20328View ArticleGoogle Scholar
- Flocco MM, Mowbray SL: Strange Bedfellows: Interactions between Acidic Side-chains in Proteins. J Mol Biol 1995, 254: 96–105. 10.1006/jmbi.1995.0602View ArticlePubMedGoogle Scholar
- Miller S, Janin J, Lesk AM, Chothia C: Interior and surface of monomeric proteins. J Mol Biol 1987, 196: 641–656. 10.1016/0022-2836(87)90038-6View ArticlePubMedGoogle Scholar
- Carugo O, Argos P: Protein-protein crystal-packing contacts. Protein Sci 1997, 6: 2261–2263.PubMed CentralView ArticlePubMedGoogle Scholar
- Janin J, Rodier F: Protein-protein interaction at crystal contacts. Proteins 1995, 23: 580–587. 10.1002/prot.340230413View ArticlePubMedGoogle Scholar
- Crosio MP, Janin J, Jullien M: Crystal packing in six crystal forms of pancreatic ribonuclease. J Mol Biol 1992, 228: 243–251. 10.1016/0022-2836(92)90503-CView ArticlePubMedGoogle Scholar
- Jelsch C, Longhi S, Cambillau C: Packing forces in nine crystal forms of cutinase. Proteins: Structure, Function and Genetics 1998, 31: 320–333. Publisher Full Text 10.1002/(SICI)1097-0134(19980515)31:3%3C320::AID-PROT8%3E3.0.CO;2-MView ArticleGoogle Scholar
- Bryant C, Strohl M, Green LK, Long HB, Alter LA, Pekar AH, Chance RE, Brems DN: Detection of an equilibrium intermediate in the folding of a monomeric insulin analog. Biochemistry 1992, 31: 5692–5698. 10.1021/bi00140a002View ArticlePubMedGoogle Scholar
- Brems DN, Brown PL, Heckenlaible LA, Frank BH: Equilibrium denaturation of insulin and proinsulin. Biochemistry 1990, 29: 9289–9293. 10.1021/bi00491a026View ArticlePubMedGoogle Scholar
- Balschmidt P, Hansen FB, Dodson EJ, Dodson GG, Korber F: Structure of porcine insulin cocrystallized with clupeine Z. Acta Crystallogr B 1991, 47(Pt 6):975–986. 10.1107/S010876819100842XView ArticlePubMedGoogle Scholar
- Dobson CM, Evans PA, Radford SE: Understanding how proteins fold: the lysozyme story so far. Trends Biochem Sci 1994, 19: 31–37. 10.1016/0968-0004(94)90171-6View ArticlePubMedGoogle Scholar
- Robinson DR, Jencks WP: The effect of compounds of the urea-guanidinium class on the activity coefficient of acetyltetraglycine ethyl ester and related compounds. J Am Chem Soc 1965, 87: 2462–2470. 10.1021/ja01089a028View ArticlePubMedGoogle Scholar
- Bhuyan AK: Protein Stabilization by Urea and Guanidine Hydrochloride. Biochemistry 2002, 41: 13386–13394. 10.1021/bi020371nView ArticlePubMedGoogle Scholar
- Kumar R, Prabhu NP, Yadaiah M, Bhuyan AK: Protein Stiffening and Entropic Stabilization in the Subdenaturing Limit of Guanidine Hydrochloride. Biophys J 2004, 87: 2656–2662. 10.1529/biophysj.104.044701PubMed CentralView ArticlePubMedGoogle Scholar
- Dias LG, Florenzano FH, Reed WF, Baptista MS, Souza SMB, Alvarez EB, Chaimovich H, Cuccovia IM, Amaral CLC, Brasil CR, Romsted LS, Politi MJ: Effect of Urea on Biomimetic Systems: Neither Water 3-D Structure Rupture nor Direct Mechanism, Simply a More "Polar Water". Langmuir 2002, 18: 319–324. 10.1021/la010176uView ArticleGoogle Scholar
- Brange J: The new era of biotech insulin analogues. Diabetologia 1997, 40 Suppl 2: S48-S53. 10.1007/s001250051400View ArticlePubMedGoogle Scholar
- Brange J, Ribel U, Hansen JF, Dodson G, Hansen MT, Havelund S, Melberg SG, Norris F, Norris K, Snel L, Sørensen AR, Voigt HO: Monomeric insulins obtained by protein engineering and their medical implications. Nature 1988, 333: 679–682. 10.1038/333679a0View ArticlePubMedGoogle Scholar
- Bentley GA, Brange J, Derewenda Z, Dodson EJ, Dodson GG, Markussen J, Wilkinson AJ, Wollmer A, Xiao B: Role of B13 glu in insulin assembly - the hexamer structure of recombinant mutant (B13 glu->gln) insulin. J Mol Biol 1992, 228: 1163–1176. 10.1016/0022-2836(92)90323-CView ArticlePubMedGoogle Scholar
- Vajo Z, Duckworth WC: Genetically Engineered Insulin Analogs: Diabetes in the New Millenium. Pharmacol Rev 2000, 52: 1–10.PubMedGoogle Scholar
- D'Arcy A, Stihle M, Kostrewa D, Dale G: Crystal engineering: a case study using the 24 kDa fragment of the DNA gyrase B subunit from Escherichia coli . Acta Crystallogr D Biol Crystallogr 1999, 55(Pt 9):1623–1625. 10.1107/S0907444999008136View ArticlePubMedGoogle Scholar
- Xiao B, Dodson GG, Dodson EJ, Clarkson J, Turkenburg J, Derewenda U, Brady L, Brange J: Mutant insulin crystal packing. J Cryst Growth 1992, 122: 144–151. 10.1016/0022-0248(92)90237-DView ArticleGoogle Scholar
- Derewenda ZS, Vekilov PG: Entropy and surface engineering in protein crystallization. Acta Crystallogr D Biol Crystallogr 2006, 62(Pt 1):116–124.View ArticlePubMedGoogle Scholar
- Daniel RM: The upper limits of enzyme thermal stability. Enzyme Microb Technol 1996, 19: 74–79. 10.1016/0141-0229(95)00174-3View ArticleGoogle Scholar
- Mammen CB, Ursby T, Cerenius Y, Thunnissen M, Als-Nielsen J, Larsen S, Liljas A: Design of a 5-Station Macromolecular Crystallography Beamline at MAX-Lab . Acta Physica Polonica A 2002, 101: 595–602.Google Scholar
- Kabsch W: Automatic processing of rotation diffraction data from crystals of initially unknown symmetry and cell constants. J Appl Cryst 1993, 26: 795–800. 10.1107/S0021889893005588View ArticleGoogle Scholar
- Vagin A, Teplyakov A: MOLREP: an Automated Program for Molecular Replacement. J Appl Cryst 1997, 30: 1022–1025. 10.1107/S0021889897006766View ArticleGoogle Scholar
- Murshudov GN, Vagin AA, Dodson EJ: Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D Biol Crystallogr 1997, 53(Pt 3):240–255. 10.1107/S0907444996012255View ArticlePubMedGoogle Scholar
- Emsley P, Cowtan K: Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr 2004, 60(Pt 12 Pt 1):2126–2132. 10.1107/S0907444904019158View ArticlePubMedGoogle Scholar
- Winn MD, Isupov MN, Murshudov GN: Use of TLS parameters to model anisotropic displacements in macromolecular refinement. Acta Crystallogr D Biol Crystallogr 2001, 57(Pt 1):122–133. 10.1107/S0907444900014736View ArticlePubMedGoogle Scholar
- DeLano WL: The PyMOL Molecular Graphics System. San Carlos, CA, USA, DeLano Scientific; 2002.Google Scholar
- Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis 1997, 18: 2714–2723. 10.1002/elps.1150181505View ArticlePubMedGoogle Scholar
- Collaborative Computational Project, Number 4: The CCP4 suite: programs for protein crystallography. Acta Crystallogr D Biol Crystallogr 1994, 50(Pt 5):760–763. 10.1107/S0907444994003112View ArticleGoogle Scholar
- Lee B, Richards FM: The interpretation of protein structures: estimation of static accessibility. J Mol Biol 1971, 55: 379–400. 10.1016/0022-2836(71)90324-XView ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucl Acids Res 2000, 28: 235–242. 10.1093/nar/28.1.235PubMed CentralView ArticlePubMedGoogle Scholar