Crystal packing in six crystal forms of pancreatic ribonuclease

A coarse-grained xDLVO model for colloidal protein-protein interactions

Colloidal protein-protein interactions (PPIs) of attractive and repulsive nature modulate the solubility of proteins, their aggregation, precipitation and crystallization. Such interactions are very important for many biotechnological processes, but are complex and hard to control, therefore, difficult to be understood in terms of measurements alone. In diluted protein solutions, PPIs can be estimated from the osmotic second virial coefficient, B22, which has been calculated using different methods and levels of theory. The most popular approach is based on the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory and its extended versions, i.e. xDLVO. Despite much efforts, these models are not fully quantitative and must be fitted to experiments, which limits their predictive value. Here, we report an extended xDLVO-CG model, which extends existing models by a coarse-grained representation of proteins and the inclusion of an additional ion-protein dispersion interaction term. We demonstrate for four proteins, i.e. lysozyme (LYZ), subtilisin (Subs), bovine serum albumin (BSA) and immunoglobulin (IgG1), that semi-quantitative agreement with experimental values without the need to fit to experimental B22 values. While most likely not the final step in the nearly hundred years of research in PPIs, xDLVO-CG is a step towards predictive PPIs calculations that are transferable to different proteins.

Determination of Protein-Protein Interactions in a Mixture of Two Monoclonal Antibodies

The coformulation of monoclonal antibody (mAb) mixtures provides an attractive route to achieving therapeutic efficacy where the targeting of multiple epitopes is necessary. Controlling and predicting the behavior of such mixtures requires elucidating the molecular basis for the self- and cross-protein-protein interactions and how they depend on solution variables. While self-interactions are now beginning to be well understood, systematic studies of cross-interactions between mAbs in solution do not exist. Here, we have used static light scattering to measure the set of self- and cross-osmotic second virial coefficients in a solution containing a mixture of two mAbs, mAbA and mAbB, as a function of ionic strength and pH. mAbB exhibits strong association at a low ionic strength, which is attributed to an electrostatic attraction that is enhanced by the presence of a strong short-ranged attraction of nonelectrostatic origin. Under all solution conditions, the measured cross-interactions are intermediate self-interactions and follow similar patterns of behavior. There is a strong electrostatic attraction at higher pH values, reflecting the behavior of mAbB. Protein-protein interactions become more attractive with an increasing pH due to reducing the overall protein net charges, an effect that is attenuated with an increasing ionic strength due to the screening of electrostatic interactions. Under moderate ionic strength conditions, the reduced cross-virial coefficient, which reflects only the energetic contribution to protein-protein interactions, is given by a geometric average of the corresponding self-coefficients. We show the relationship can be rationalized using a patchy sphere model, where the interaction energy between sites i and j is given by the arithmetic mean of the i-i and j-j interactions. The geometric mean does not necessarily apply to all mAb mixtures and is expected to break down at a lower ionic strength due to the nonadditivity of electrostatic interactions.

Understanding the fabric of protein crystals: computational classification of biological interfaces and crystal contacts

Modern structural biology still draws the vast majority of information from crystallography, a technique where the objects being investigated are embedded in a crystal lattice. Given the complexity and variety of those objects, it becomes fundamental to computationally assess which of the interfaces in the lattice are biologically relevant and which are simply crystal contacts. Since the mid-1990s, several approaches have been applied to obtain high-accuracy classification of crystal contacts and biological protein–protein interfaces. This review provides an overview of the concepts and main approaches to protein interface classification: thermodynamic estimation of interface stability, evolutionary approaches based on conservation of interface residues, and co-occurrence of the interface across different crystal forms. Among the three categories, evolutionary approaches offer the strongest promise for improvement, thanks to the incessant growth in sequence knowledge. Importantly, protein interface classification algorithms can also be used on multimeric structures obtained using other high-resolution techniques or for protein assembly design or validation purposes. A key issue linked to protein interface classification is the identification of the biological assembly of a crystal structure and the analysis of its symmetry. Here, we highlight the most important concepts and problems to be overcome in assembly prediction. Over the next few years, tools and concepts of interface classification will probably become more frequently used and integrated in several areas of structural biology and structural bioinformatics. Among the main challenges for the future are better addressing of weak interfaces and the application of interface classification concepts to prediction problems like protein–protein docking.

Supplementary information: Supplementary data are available at Bioinformatics online.

Contact:guido.capitani@psi.ch

SdsA polymorph isolation and improvement of their crystal quality using nonconventional crystallization techniques

SdsA, a sodium dodecyl sulfate hydrolase, fromPseudomonas aeruginosawas crystallized in three different crystal polymorphs and their three-dimensional structure was determined. The different polymorphs present different crystal packing habits. One of the polymorphs suggests the existence of a tetramer, an oligomeric state not observed previously, while the crystal packing of the remaining two polymorphs obstructs the active site entrance but stabilizes flexible regions of the protein. Nonconventional crystallization methods that minimize convection, such as counterdiffusion in polyvinyl alcohol gel coupled with the influence of a 500 MHz (10.2 T) magnetic field, were necessary to isolate the poorest diffracting polymorph and increase its internal order to determine its structure by X-ray diffraction. The results obtained show the effectiveness of nonconventional crystallographic methods to isolate different crystal polymorphs.

Simulation of protein association: Kinetic pathways towards crystal contacts

We conducted molecular dynamics simulations combined with distance-based umbrella sampling and forward flux sampling to investigate the early stages of protein crystallization. Formation of contacts with long-range interactions and/or an exposed position on the protein surface was kinetically preferred over more stable hydrophobic contacts with a shorter attractive range, while the thermodynamic stability of the protein crystal was provided by hydrophobic interactions. Contacts with a large interaction area showed complex dissociation pathways that were not detected by distance-based umbrella sampling. Instead, forward flux sampling simulations of contact dissociation identified long-range attractive interactions.

Effects of protein engineering and rational mutagenesis on crystal lattice of single chain antibody fragments

Protein crystallization is dependent upon, and sensitive to, the intermolecular contacts that assist in ordering proteins into a three-dimensional lattice. Here we used protein engineering and mutagenesis to affect the crystallization of single chain antibody fragments (scFvs) that recognize the EE epitope (EYMPME) with high affinity. These hypercrystallizable scFvs are under development to assist difficult proteins, such as membrane proteins, in forming crystals, by acting as crystallization chaperones. Guided by analyses of intermolecular crystal lattice contacts, two second-generation anti-EE scFvs were produced, which bind to proteins with installed EE tags. Surprisingly, although noncomplementarity determining region (CDR) lattice residues from the parent scFv framework remained unchanged through the processes of protein engineering and rational design, crystal lattices of the derivative scFvs differ. Comparison of energy calculations and the experimentally-determined lattice interactions for this basis set provides insight into the complexity of the forces driving crystal lattice choice and demonstrates the availability of multiple well-ordered surface features in our scFvs capable of forming versatile crystal contacts.

Depletion of IK causes mitotic arrest through aberrant regulation of mitotic kinases and phosphatases

IK is known to inhibit the expression of major histocompatibility complex (MHC) class II antigen, but other cellular functions of IK remain to be uncovered. In this study, IK depletion caused misalignment of chromosomes through an increase in Aurora A and PLK1 phosphorylation, which was mediated by a decrease in PP1 and PP2A activities. On the other hand, the treatment of a dual inhibitor against CDK and Aurora kinases overrode IK depletion-induced mitotic arrest through the activation of phosphatase activity. These findings imply that IK is an essential protein for achieving correct mitotic progress through the regulation of mitotic kinases and phosphatases.

Polymorphic Protein Crystal Growth: Influence of Hydration and Ions in Glucose Isomerase

Crystal polymorphs of glucose isomerase were examined to characterize the properties and to quantify the energetics of protein crystal growth. Transitions of polymorph stability were measured in poly(ethylene glycol)/NaCl solutions, and one transition point was singled out for more detailed quantitative analysis. Single crystal x-ray diffraction was used to confirm space groups and identify complementary crystal structures. Crystal polymorph stability was found to depend on the NaCl concentration, with stability transitions requiring > 1 M NaCl combined with a low concentration of PEG. Both salting-in and salting-out behavior was observed and was found to differ for the two polymorphs. For NaCl concentrations above the observed polymorph transition, the increase in solubility of the less stable polymorph together with an increase in the osmotic second virial coefficient suggests that changes in protein hydration upon addition of salt may explain the experimental trends. A combination of atomistic and continuum models was employed to dissect this behavior. Molecular dynamics simulations of the solvent environment were interpreted using quasi-chemical theory to understand changes in protein hydration as a function of NaCl concentration. The results suggest that protein surface hydration and Na+ binding may introduce steric barriers to contact formation, resulting in polymorph selection.

Self-interaction chromatography: a tool for the study of protein-protein interactions in bioprocessing environments

We describe a new protein characterization technique called self-interaction chromatography (SIC), which exploits the specificity of protein-protein interactions that is common to protein aggregates and enables the rapid screening of protein formulation additives as physical stabilizers against aggregation. This technique also enables the identification of specific interaction sites and the determination of their relative importance for self-association. Mannitol, glycine, and dextran 40 were tested for their stabilizing effect toward the model protein lysozyme. Dextran 40 exhibited a poor stabilizing effect. While mannitol stabilized both the native and acid-denatured forms of lysozyme, glycine stabilized the native form with respect to the denatured species. These results are in good agreement with findings in the formulation literature. The SIC shows tremendous potential as a rapid formulation development tool. We also screened two putative interaction sites for involvement in the self-association of lysozyme and estimated the associated binding energies using a binding isotherm model that we developed. The sites screened consisted of residues 41-48 and 125-128 and were selected based on their apparent importance in forming crystal contacts in several different crystal forms of lysozyme. Of the two sites, only residues 125-128 were found to influence self-association under the conditions we employed. Because the success of this technique depends on the exploitation of self-interactions between native species, several important applications are also suggested such as separating native from misfolded or variant species and probing site utilization in aggregation versus crystallization phenomena.

Regulation of protein function: crystal packing interfaces and conformational dimerization

The accepted view of interprotein electron transport involves molecules diffusing between donor and acceptor redox sites. An emerging alternative hypothesis is that efficient long-range electron transport can be achieved through proteins arranged in supramolecular assemblies. In this study, we have investigated the crystal packing interfaces in three crystal forms of plastocyanin, an integral component of the photosynthetic electron transport chain, and discuss their potential relevance to in vivo supramolecular assemblies. Symmetry-related protein chains within these crystals have Cu-Cu separations of <25 A, a distance that readily supports electron transfer. In one structure, the plastocyanin molecule exists in two forms in which a backbone displacement coupled with side chain rearrangements enables the modulation of protein-protein interfaces.

Modification of protein crystal packing by systematic mutations of surface residues: implications on biotemplating and crystal porosity

Bioinspired nano-scale biotemplating for the development of novel composite materials has recently culminated in several demonstrations of nano-structured hybrid materials. Protein crystals, routinely prepared for the elucidation of protein 3D structures by X-ray crystallography, present an ordered and highly accurate 3D array of protein molecules. Inherent to the 3D arrangement of the protein "building blocks" in the crystal, a complementary 3D array of interconnected cavities--voids array, exhibiting highly ordered porosity is formed. The porous arrays of protein crystal may serve as a nano-structured, accurate biotemplate by a "filling" process. These cavities arrays are shaped by the mode of protein packing throughout the crystallization process. Here we propose and demonstrate feasibility of targeting site specific mutations to modify protein's surface to affect protein crystal packing, enabling the generation of a series of protein crystal "biotemplates" all originating from same parent protein. The selection of these modification sites was based on in silico analysis of protein-protein interface contact areas in the parent crystal. The model protein selected for this study was the N-terminal type II cohesin from the cellulosomal scaffold in ScaB subunit of Acetivibrio cellulolyticus and mutations were focused on lysine residues involved in protein packing as prime target. The impact of systematically mutating these lysine residues on protein packing and its resulting interconnected cavities array were found to be most significant when surface lysine residues were substituted to tryptophan residues. Our results demonstrate the feasibility of using pre-designed site directed mutations for the generation of a series of protein crystal biotemplates from a "parent" protein.

Molecular dynamics characterization of protein crystal contacts in aqueous solutions

We employ nonequilibrium molecular dynamics simulation to characterize the effective interactions between lysozyme molecules involved in the formation of two hydrophobic crystal contacts. We show that the effective interactions between crystal contacts do not exceed a few kT, the range of the attractive part of the potential is less than 4 angstroms, and, within this range, there is a significant depletion of water density between two protein contacts. Our findings highlight the different natures of protein crystallization and protein recognition processes.

Protein–protein interaction and quaternary structure

Protein–protein recognition plays an essential role in structure and function. Specific non-covalent interactions stabilize the structure of macromolecular assemblies, exemplified in this review by oligomeric proteins and the capsids of icosahedral viruses. They also allow proteins to form complexes that have a very wide range of stability and lifetimes and are involved in all cellular processes. We present some of the structure-based computational methods that have been developed to characterize the quaternary structure of oligomeric proteins and other molecular assemblies and analyze the properties of the interfaces between the subunits. We compare the size, the chemical and amino acid compositions and the atomic packing of the subunit interfaces of protein–protein complexes, oligomeric proteins, viral capsids and protein–nucleic acid complexes. These biologically significant interfaces are generally close-packed, whereas the non-specific interfaces between molecules in protein crystals are loosely packed, an observation that gives a structural basis to specific recognition. A distinction is made within each interface between a core that contains buried atoms and a solvent accessible rim. The core and the rim differ in their amino acid composition and their conservation in evolution, and the distinction helps correlating the structural data with the results of site-directed mutagenesis and in vitro studies of self-assembly.

Crystallographic characterization of two novel crystal forms of human insulin induced by chaotropic agents and a shift in pH

Background

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.

Results

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 C222₁ and C2 respectively. The orthorhombic crystals were obtained at pH 6.5 and contained three insulin hexamers in R₆ conformation in the asymmetric unit whilst the monoclinic C2 crystals were obtained at pH 7.0 and contained one R₆ 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.

Conclusion

The change in space group from C222₁ 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 C222₁ 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.

Peptide segments in protein-protein interfaces

An important component of functional genomics involves the understanding of protein association. The interfaces resulting from protein-protein interactions - (i) specific, as represented by the homodimeric quaternary structures and the complexes formed by two independently occurring protein components, and (ii) non-specific, as observed in the crystal lattice of monomeric proteins - have been analysed on the basis of the length and the number of peptide segments. In 1000 A2 of the interface area, contributed by a polypeptide chain, there would be 3.4 segments in homodimers, 5.6 in complexes and 6.3 in crystal contacts. Concomitantly, the segments are the longest (with 8.7 interface residues) in homodimers. Core segments (likely to contribute more towards binding) are more in number in homodimers (1.7) than in crystal contacts (0.5), and this number can be used as one of the parameters to distinguish between the two types of interfaces. Dominant segments involved in specific interactions, along with their secondary structural features, are enumerated.

Conservation and relative importance of residues across protein-protein interfaces

A core region surrounded by a rim characterizes biological interfaces. We ascertain the importance of the core by showing the sequence entropies of the residues comprising the core to be smaller than those in the rim. Such a distinction is not seen in the 2-fold-related, nonphysiological interfaces formed in crystal lattices of monomeric proteins, thereby providing a procedure for characterizing the oligomeric state from crystal structures of protein molecules. This method is better than those that rely on the comparison of the sequence entropies in the interface and the rest of the protein surface, especially in cases where the surface harbors additional binding sites. To a good approximation there is a correlation between the accessible surface area lost because of complexation and DeltaDeltaG values obtained through alanine-scanning mutagenesis (26-38 cal per A(2) of the surface buried) for residues located in the core, a relationship that is not discernable for rim residues. If, however, a residue participates in hydrogen bonding across the interface, the extent of stabilization is 52 cal/mol per 1 A(2) of the nonpolar surface area buried by the residue. As opposed to an amino acid classification used earlier, an environment-based grouping of residues yields a better discrimination in the sequence entropy between the core and the rim.

The Limit of Accuracy of Protein Modeling: Influence of Crystal Packing on Protein Structure

The size of the protein database (PDB) makes it now feasible to arrive at statistical conclusions regarding structural effects of crystal packing. These effects are relevant for setting upper practical limits of accuracy on protein modeling. Proteins whose crystals have more than one molecule in the asymmetric unit or whose structures were determined at least twice by X-ray crystallography were paired and their differences analyzed. We demonstrate a clear influence of crystal environment on protein structure, including backbone conformations, hinge-like motions and side-chain conformations. The positions of surface water molecules tend to be variable in different crystal environments while those of ligands are not. Structures determined by independent groups vary more than structures determined by the same authors. The use of different refinement methods is a major source for this effect. Our pair-wise analysis derives a practical limit to the accuracy of protein modeling. For different crystal forms, the limit of accuracy (C(alpha), root-mean-square deviation (RMSD)) is approximately 0.8A for the entire protein, which includes approximately 0.3A due to crystal packing. For organized secondary elements, the upper limit of C(alpha) RMSD is 0.5-0.6A while for loops or protein surface it reaches 1.0A. Twenty percent of exposed side- chains exhibit different chi(1+2) conformations with approximately half of the effect also resulting from crystal packing. A web based tool for analysis and graphic presentation of surface areas of crystal contacts is available (http://ligin.weizmann.ac.il/cryco).

A Dissection of Specific and Non-specific Protein–Protein Interfaces

We compare the geometric and physical-chemical properties of interfaces involved in specific and non-specific protein-protein interactions in crystal structures reported in the Protein Data Bank. Specific interactions are illustrated by 70 protein-protein complexes and by subunit contacts in 122 homodimeric proteins; non-specific interactions are illustrated by 188 pairs of monomeric proteins making crystal-packing contacts selected to bury more than 800 A2 of protein surface. A majority of these pairs have 2-fold symmetry and form "crystal dimers" that cannot be distinguished from real dimers on the basis of the interface size or symmetry. The chemical and amino acid compositions of the large crystal-packing interfaces resemble the protein solvent-accessible surface. These interfaces are less hydrophobic than in homodimers and contain much fewer fully buried atoms. We develop a residue propensity score and a hydrophobic interaction score to assess preferences seen in the chemical and amino acid compositions of the different types of interfaces, and we derive indexes to evaluate the atomic packing, which we find to be less compact at non-specific than at specific interfaces. We test the capacity of these parameters to identify homodimeric proteins in crystal structures, and show that a simple combination of the non-polar interface area and the fraction of buried interface atoms assigns the quaternary structure of 88% of the homodimers and 77% of the monomers in our data set correctly. These success rates increase to 93-95% when the residue propensity score of the interfaces is taken into consideration.

Designing proteins to crystallize through beta-strand pairing

Inherent difficulties in growing protein crystals are major concerns within structural biology and particularly in structural proteomics. Here, we describe a novel approach of engineering target proteins by surface mutagenesis to increase the odds of crystallizing the molecules. To this end, we have exploited our recent triad-hypothesis using proteins with crystallographically defined beta-structures as the principal models. Crystal packing analyses of 182 protein structures belonging to 21 different superfamilies implied that the propensities to crystallize could be engineered into target proteins by replacing short segments, 5-6 residues, of their beta-strands with 'cassettes' of suitable packing motifs. These packing motifs will generate specific crystal packing interactions that promote crystallization. Key features of the primary and tertiary structures of such packing motifs have been identified for immunoglobulins. Further, packing motifs have been engineered successfully into six model antibodies without disturbing their capabilities to be produced, their immunoreactivity and their overall structure. Preliminary crystallization analyses have also been performed. Taken together, the procedures outline a rational protocol for crystallizing proteins by surface mutagenesis. The importance of these findings is discussed in relation to the crystallization of proteins in general.

Structure of the C-terminal domains of merozoite surface protein-1 from Plasmodium knowlesi reveals a novel histidine binding site

The protozoan parasite Plasmodium causes malaria, with hundreds of millions of cases recorded annually. Protection against malaria infection can be conferred by antibodies against merozoite surface protein (MSP)-1, making it an attractive vaccine candidate. Here we present the structure of the C-terminal domains of MSP-1 (known as MSP-1(19)) from Plasmodium knowlesi. The structure reveals two tightly packed epidermal growth factor-like domains oriented head to tail. In domain 1, the molecule displays a histidine binding site formed primarily by a highly conserved tryptophan. The protein carries a pronounced overall negative charge primarily due to the large number of acidic groups in domain 2. To map protein binding surfaces on MSP-1(19), we have analyzed the crystal contacts in five different crystal environments, revealing that domain 1 is highly preferred in protein-protein interactions. A comparison of MSP-1(19) structures from P. knowlesi, P. cynomolgi, and P. falciparum shows that, although the overall protein folds are similar, the molecules show significant differences in charge distribution. We propose the histidine binding site in domain 1 as a target for inhibitors of protein binding to MSP-1, which might prevent invasion of the merozoite into red blood cells.

Conservation helps to identify biologically relevant crystal contacts

Some crystal contacts are biologically relevant, most are not. We assess the utility of combining measures of size and conservation to discriminate between biological and non-biological contacts. Conservation and size information is calculated for crystal contacts in 53 families of homodimers and 65 families of monomers. Biological contacts are shown to be usually conserved and typically the largest contact in the crystal. A range of neural networks accepting different combinations and encodings of this information is used to answer the following questions: (1) is a given crystal contact biological, and (2) given all crystal contacts in a homodimer, which is the biological one? Predictions for (1) are performed on both homodimer and monomer datasets. The best performing neural network combined size and conservation inputs. For the homodimers, it correctly classified 48 out of 53 biological contacts and 364 out of 366 non-biological contacts, giving a combined accuracy of 98.3 %. A more robust performance statistic, the phi-coefficient, which accounts for imbalances in the dataset, gave a value of 0.92. Taking all 535 non-biological contacts from the 65 monomers, this predictor made erroneous classifications only 4.3 % of the time. Predictions for (2) were performed on homodimers only. The best performing network achieved a prediction accuracy of 98.1 % using size information alone. We conclude that in answering question (1) size and conservation combined discriminate biological from non-biological contacts better than either measure alone. For answering question (2), we conclude that in our dataset size is so powerful a discriminant that conservation adds little predictive benefit.

The analysis of the human high affinity IgE receptor Fc epsilon Ri alpha from multiple crystal forms

We have solved the structure of the human high affinity IgE receptor, Fc epsilon RI alpha, in six different crystal forms, showing the structure in 15 different chemical environments. This database of structures shows no change in the overall shape of the molecule, as the angle between domains 1 and 2 (D1 and D2) varies little across the ensemble. However, the receptor has local conformational variability in the C' strand of D2 and in the BC loop of D1. In every crystal form, a residue inserts between tryptophan residues 87 and 110, mimicking the position of a proline from the IgE ligand. The different crystal forms reveal a distribution of carbohydrates lining the front and back surfaces of the structure. An analysis of crystal contacts in the different forms indicates regions where the molecule interacts with other proteins, and reveals a potential new binding site distal to the IgE binding site. The results of this study point to new directions for the design of molecules to inhibit the interaction of Fc epsilon RI alpha with its natural ligand and thus to prevent a primary step in the allergic response.

The BPTI decamer observed in acidic pH crystal forms pre-exists as a stable species in solution

Bovine pancreatic trypsin inhibitor (BPTI) crystallizes under acidic pH conditions in the presence of thiocyanate, chloride and sulfate ions, yielding three different polymorphs in P2(1), P6(4)22 and P6(3)22 space groups, respectively. In all three crystal forms, the same decamer is found in the packing (ten BPTI molecules organized through two perpendicular 2-fold and 5-fold axes as a well-defined and compact object) in contrast to the monomeric crystal forms observed at basic pH conditions. The crystallization of BPTI under acidic conditions (pH 4.5) was investigated by small angle X-ray scattering with both under- and supersaturated BPTI solutions. Data showed the oligomerization of BPTI molecules under all investigated conditions. Accordingly, various mixtures of discrete oligomers (n=1 to 10) were considered. Calculated scattering curves were obtained using models based on the crystallographic structures, and the experimental patterns were analyzed as a linear combination of the model curves using a non-linear curve fitting procedure. The results, confirmed by gel filtration experiments, unambiguously demonstrate the co-existence of two different BPTI particles in solution: a monomer and a decamer, with no evidence of any other intermediates. Moreover, using both approaches, the fraction of decamers was found to increase with increasing salt concentration, even beyond the solubility curve. We therefore propose that at acidic pH, BPTI crystallizes following a two step process: decamers are first built in under- and supersaturated solutions, upon which crystal growth proceeds by decamer stacking. Indeed, those BPTI crystals should best be described as "BPTI decamer" crystals.

A soft, mean-field potential derived from crystal contacts for predicting protein-protein interactions

We derive a series of novel mean-field potentials from statistical analyses of protein-protein contact regions in crystal structures. These potentials are parameterized in terms of the number of contacts made by an atom in an interface region. Such an explicit number dependence avoids the pairwise assumption and is intrinsically softer than distance-based approaches. It appears well suited to protein-protein docking applications, for which detailed interface geometry is generally lacking. In tests including protein complex reconstitution and docking of independently determined protein structures, we show that a hydrophobic potential of this type performs remarkably well, identifying native-like complexes by their favourable potential energies and in several cases demonstrating a recognition energy gap of 4-8 kcal/mol according to the system.

Molecular origins of osmotic second virial coefficients of proteins

The thermodynamic properties of protein solutions are determined by the molecular interactions involving both solvent and solute molecules. A quantitative understanding of the relationship would facilitate more systematic procedures for manipulating the properties in a process environment. In this work the molecular basis for the osmotic second virial coefficient, B22, is studied; osmotic effects are critical in membrane transport, and the value of B22 has also been shown to correlate with protein crystallization behavior. The calculations here account for steric, electrostatic, and short-range interactions, with the structural and functional anisotropy of the protein molecules explicitly accounted for. The orientational dependence of the protein interactions is seen to have a pronounced effect on the calculations; in particular, the relatively few protein-protein configurations in which the apposing surfaces display geometric complementarity contribute disproportionately strongly to B22. The importance of electrostatic interactions is also amplified in these high-complementarity configurations. The significance of molecular recognition in determining B22 can explain the correlation with crystallization behavior, and it suggests that alteration of local molecular geometry can help in manipulating protein solution behavior. The results also have implications for the role of protein interactions in biological self-organization.

Molecular packing in virus crystals: geometry, chemistry, and biology

An automated procedure was developed to determine the geometrical and chemical interactions of crystalline virus particles using the crystal parameters, particle position, orientation, and atomic coordinates for an icosahedral asymmetric unit. Two applications of the program are reported: (1) An analysis of a novel pseudo-rhombohedral (R32) symmetry present in the monoclinic crystal lattices of both Nodamura Virus (NOV) and Coxsackie virus B3 (CVB3). The study shows that in both cases the interactions between particles is substantially increased by minor deviations from exact R32 symmetry and that only particles with the proper ratio of dimensions along twofold and fivefold symmetry axes (such as southern bean mosaic virus) can achieve comparable buried surface area in the true R32 space group. (2) An attempt was made to correlate biological function with remarkably conserved interparticle contact regions found in different crystal forms of three members of the nodavirus family, NOV, Flock House Virus (FHV), and Black Beetle Virus (BBV). Mutational evidence implicates the quasi-threefold region on the viral surface in receptor binding in nodaviruses and this region is dominant in particle contacts in all three virus crystals. Examination of particle contacts in numerous crystal structures of viruses in the picornavirus super-family showed that portions of the capsid surface known to interact with a receptor or serve as an epitope for monoclonal antibodies frequently stabilize crystal contacts.

A structural comparison of the colicin immunity proteins Im7 and Im9 gives new insights into the molecular determinants of immunity-protein specificity

We report the first detailed comparison of two immunity proteins which, in conjunction with recent protein engineering data, begins to explain how these structurally similar proteins are able to bind and inhibit the endonuclease domain of colicin E9 (E9 DNase) with affinities that differ by 12 orders of magnitude. In the present work, we have determined the X-ray structure of the Escherichia coli colicin E7 immunity protein Im7 to 2.0 A resolution by molecular replacement, using as a trial model the recently determined NMR solution structure of Im9. Whereas the two proteins adopt similar four-helix structures, subtle structural differences, in particular involving a conserved tyrosine residue critical for E9 DNase binding, and the identity of key residues in the specificity helix, lie at the heart of their markedly different ability to bind the E9 DNase. Two other crystal structures were reported recently for Im7; in one, Im7 was a monomer and was very similar to the structure reported here, whereas in the other it was a dimer to which functional significance was assigned. Since this previous work suggested that Im7 could exist either as a monomer or a dimer, we used analytical ultracentrifugation to investigate this question further. Under a variety of solution conditions, we found that Im7 only ever exists in solution as a monomer, even up to protein concentrations of 15 mg/ml, casting doubt on the functional significance of the crystallographically observed dimer. This work provides a structural framework with which we can understand immunity-protein specificity, and in addition we believe it to be the first successfully refined crystal structure solved by molecular replacement using an NMR trial model with less than 100% sequence identity.

No salting-in of lysozyme chloride observed at low ionic strength over a large range of pH

Solubility of lysozyme chloride was determined in the absence of added salt and in the presence of 0.05-1.2 M NaCl, starting from isoionic lysozyme, which was then brought to pH values from 9 to 3 by addition of HCl. The main observation is the absence of a salting-in region whatever the pH studied. This is explained by a predominant electrostatic screening of the positively charged protein and/or by adsorption of chloride ions by the protein. The solubility increases with the protein net charge at low ionic strength, but the reverse is observed at high ionic strength. The solubility of lysozyme chloride seems to become independent of ionic strength at pH approximately 9.5, which is interpreted as a shift of the isoionic pH (10.8) to an isoelectric pH due to chloride binding. The crystallization at very low ionic strength, where lysozyme crystallizes at supersaturation values as low as 1.1, amplifies the effect of pH on protein solubility. Understanding the effect of the net charge and of ionic strength on protein-protein interactions is valuable not only for protein crystal growth but more generally also for the formation of protein-protein or protein-ligand complexes.

Protein-protein crystal-packing contacts

Protein-protein contacts in monomeric protein crystal structures have been analyzed and compared to the physiological protein-protein contacts in oligomerization. A number of features differentiate the crystal-packing contacts from the natural contacts occurring in multimeric proteins. The area of the protein surface patches involved in packing contacts is generally smaller and its amino acid composition is indistinguishable from that of the protein surface accessible to the solvent. The fraction of protein surface in crystal contacts is very variable and independent of the number of packing contacts. The thermal motion at the crystal packing interface and that of the protein core, even for large packing interfaces, though the tendency is to be closer to that of the core. These results suggest that protein crystallization depends on random protein-protein interactions, which have little in common with physiological protein-protein recognition processes, and that the possibility of engineering macromolecular crystallization to improve crystal quality could be widened.

Extent and nature of contacts between protein molecules in crystal lattices and between subunits of protein oligomers

A survey was compiled of several characteristics of the intersubunit contacts in 58 oligomeric proteins, and of the intermolecular contracts in the lattice for 223 protein crystal structures. The total number of atoms in contact and the secondary structure elements involved are similar in the two types of interfaces. Crystal contact patches are frequently smaller than patches involved in oligomer interfaces. Crystal contacts result from more numerous interactions by polar residues, compared with a tendency toward nonpolar amino acids at oligomer interfaces. Arginine is the only amino acid prominent in both types of interfaces. Potentials of mean force for residue-residue contacts at both crystal and oligomer interfaces were derived from comparison of the number of observed residue-residue interactions with the number expected by mass action. They show that hydrophobic interactions at oligomer interfaces favor aromatic amino acids and methionine over aliphatic amino acids; and that crystal contacts form in such a way as to avoid inclusion of hydrophobic interactions. They also suggest that complex salt bridges with certain amino acid compositions might be important in oligomer formation. For a protein that is recalcitrant to crystallization, substitution of lysine residues with arginine or glutamine is a recommended strategy.

Protein crystallization

Crystallization is necessary to obtain the three-dimensional structure of proteins and nucleic acids; it often represents the bottleneck in structure determination. Our understanding of crystallization mechanisms is still incomplete. In this review, we emphasize fundamental aspects of the crystallization process. Protein-protein contacts in crystals are complex, involving a delicate balance of specific and nonspecific interactions. Depending on solution conditions, these interactions can lead to nucleation of crystals or to amorphous aggregation; this stage of crystallization has been successfully studied by light scattering. Post-nucleation crystal growth may proceed by mechanisms involving crystal defects or two-dimensional nucleation, as observed by atomic force and interference microscopy. Cessation of growth has been observed but remains incompletely understood. Impurities may play important roles during all stages of crystallization. Phase diagrams can guide optimization of conditions for nucleation and subsequent crystal growth; a theoretical understanding relating these to the intermolecular interactions is beginning to develop.

Three quaternary structures for a single protein

The structure of a multisubunit protein (immunoglobulin light chain) was solved in three crystal forms, differing only in the solvent of crystallization. The three structures were obtained at high ionic strength and low pH, high ionic strength and high pH, and low ionic strength and neutral pH. The three resulting "snapshots" of possible structures show that their variable-domain interactions differ, reflecting their stabilities under specific solvent conditions. In the three crystal forms, the variable domains had different rotational and translational relationships, whereas no alteration of the constant domains was found. The critical residues involved in the observed effect of the solvent are tryptophans and histidines located between the two variable domains in the dimeric structure. Tryptophan residues are commonly found in interfaces between proteins and their subunits, and histidines have been implicated in pH-dependent conformation changes. The quaternary structure observed for a multisubunit protein or protein complex in a crystal may be influenced by the interactions of the constituents within the molecule or complex and/or by crystal packing interactions. The comparison of buried surface areas and hydrogen bonds between the domains forming the molecule and between the molecules forming the crystals suggest that, for this system, the interactions within the molecule are most likely the determining factors.

Simulations of reversible protein aggregate and crystal structure

We simulated the structure of reversible protein aggregates as a function of protein surface characteristics, protein-protein interaction energies, and the entropic penalty accompanying the immobilization of protein in a solid phase. These simulations represent an extension of our previous work on kinetically irreversible protein aggregate structure and are based on an explicit accounting of the specific protein-protein interactions that occur within reversible aggregates and crystals. We considered protein monomers with a mixture of hydrophobic and hydrophilic surface regions suspended in a polar solvent; the energetic driving force for aggregation is provided by the burial of solvent-exposed hydrophobic surface area. We analyzed the physical properties of the generated aggregates, including density, protein-protein contact distributions, solvent accessible surface area, porosity, and order, and compared our results with the protein crystallization literature as well as with the kinetically irreversible case. The physical properties of reversible aggregates were consonant with those observed for the irreversible aggregates, although in general, reversible aggregates were more stable energetically and were more crystal-like in their order content than their irreversible counterparts. The reversible aggregates were less dense than the irreversible aggregates, indicating that the increased energetic stability is derived primarily from the optimality rather than the density of the packing in the solid phase. The extent of hydrophobic protein-protein contacts and solvent-exposed surface area within the aggregate phase depended on the aggregation pathway: reversible aggregates tended to have a greater proportion of hydrophobic-hydrophobic contacts and a smaller fraction of hydrophobic solvent-exposed surface area. Furthermore, the arrangement of hydrophobic patches on the protein surface played a major role in the distribution of protein contacts and solvent content. This was readily reflected in the order of the aggregates: the greater the contiguity of the hydrophobic patches on the monomer surface, the less ordered the aggregates became, despite the opportunities for rearrangement offered by a reversible pathway. These simulations have enhanced our understanding of the impact of protein structural motifs on aggregate properties and on the demarcation between aggregation and crystallization.

Protein-protein interaction at crystal contacts

Packing contacts are crystal artifacts, yet they make use of the same forces that govern specific recognition in protein-protein complexes and oligomeric proteins. They provide examples of a nonspecific protein-protein interaction which can be compared to biologically relevant ones. We evaluate the number and size of pairwise interfaces in 152 crystal forms where the asymmetric unit contains a monomeric protein. In those crystal forms that have no element of 2-fold symmetry, we find that molecules form 8 to 10 pairwise interfaces. The total area of the surface buried on each molecule is large, up to 4400 A2. Pairwise interfaces bury 200-1200 A2, like interfaces generated at random in a computer simulation, and less than interfaces in protease-inhibitor or antigen-antibody complexes, which bury 1500 A2 or more. Thus, specific contacts occurring in such complexes extend over a larger surface than nonspecific ones. In crystal forms with 2-fold symmetry, pairwise interfaces are fewer and larger on average than in the absence of 2-fold symmetry. Some bury 1500-2500 A2, like interfaces in oligomeric proteins, and create "crystal oligomers" which may have formed in the solution before crystallizing.

Principles of protein-protein recognition from structure to thermodynamics

Specific recognition is illustrated by X-ray structures of protease-inhibitor, antigen-antibody and other high affinity complexes including five electron transfer complexes. We attempt to give a physical definition to affinity and specificity on the basis of these data. In a protein-protein complex, specific recognition results from the assembly of complementary surfaces into well-packed interfaces that cover about 1500 A2 and contain about ten hydrogen bonds. These interfaces are larger than between molecules in protein crystals, and smaller than between subunits in oligomeric proteins. We relate the size and chemical nature of interfaces in complexes to the thermodynamical parameters that characterize affinity: the heat capacity and free enthalpy (Gibbs energy) of dissociation at equilibrium, the activation free enthalpy for the dissociation reaction. The same structural and thermodynamical parameters are inadequate for representing the specificity of recognition. We propose instead to describe specificity with the help of statistical physics, and we illustrate the application of the random energy model to antigen-antibody recognition by analyzing results of computer simulations by docking.

The structures of RNase A complexed with 3'-CMP and d(CpA): active site conformation and conserved water molecules

The interactions of RNase A with cytidine 3'-monophosphate (3'-CMP) and deoxycytidyl-3',5'-deoxyadenosine (d(CpA)) were analyzed by X-ray crystallography. The 3'-CMP complex and the native structure were determined from trigonal crystals, and the d(CpA) complex from monoclinic crystals. The differences between the overall structures are concentrated in loop regions and are relatively small. The protein-inhibitor contacts are interpreted in terms of the catalytic mechanism. The general base His 12 interacts with the 2' oxygen, as does the electrostatic catalyst Lys 41. The general acid His 119 has 2 conformations (A and B) in the native structure and is found in, respectively, the A and the B conformation in the d(CpA) and the 3'-CMP complex. From the present structures and from a comparison with RNase T1, we propose that His 119 is active in the A conformation. The structure of the d(CpA) complex permits a detailed analysis of the downstream binding site, which includes His 119 and Asn 71. The comparison of the present RNase A structures with an inhibitor complex of RNase T1 shows that there are important similarities in the active sites of these 2 enzymes, despite the absence of any sequence homology. The water molecules were analyzed in order to identify conserved water sites. Seventeen water sites were found to be conserved in RNase A structures from 5 different space groups. It is proposed that 7 of those water molecules play a role in the binding of the N-terminal helix to the rest of the protein and in the stabilization of the active site.

NMR study of the positions of His-12 and His-119 in the ribonuclease A-uridine vanadate complex

The binding of uridine vanadate to ribonuclease A has been investigated by one- and two-dimensional 1H NMR. The homonuclear Nuclear Overhauser and exchange spectroscopy spectrum of the uridine vanadate/RNase A complex exhibits cross peaks between both the C5H and C6H protons of uridine vanadate and the H epsilon 1 proton of His-12 of ribonuclease A. These cross peaks suggest that the H epsilon 1 proton of His-12 is in the vicinity of the uracil base of uridine vanadate, as observed in the crystallographic structure of the uridine vanadate/RNase A complex. However, no cross peaks are observed between the C5H and C6H protons of uridine vanadate and the H epsilon 1 proton of His-119 of ribonuclease A, although they were predicted based upon the distances calculated from coordinates of the crystallographic structure of the complex. These results suggest that there is a significant difference between the positioning of the His-119 side chain in the solution and in the crystallographic structures.

Analysis of electrostatic interactions in ribonuclease A monoclinic crystal

The electrostatic field in ribonuclease A monoclinic crystal is computed. A calculation of the electrostatic free energy of interaction between one molecule and its neighbours in the crystal lattice is performed in three crystal planes. It is shown that at the pH value of crystallization for the polymorphic form under investigation the total electrostatic free energy in the monoclinic crystal of ribonuclease A already crystallized is equal to zero. However, the electrostatic free energy along the crystal axis has different values. Electrostatic forces play a negative role for the crystal stability along the alpha-axis while along the other axis the electrostatic interaction supports the crystal state. The electrostatic potential distribution of a ribonuclease A monoclinic crystal at pH 6 is plotted and it is demonstrated that these interactions have a different character in each crystallographic plane.

Comparison of the structures and the crystal contacts of trypanosomal triosephosphate isomerase in four different crystal forms

Triosephosphate isomerase (TIM) is a dimeric enzyme consisting of 2 identical subunits. Trypanosomal TIM can be crystallized in 4 different spacegroups: P2(1)2(1)2(1), C2(big cell), C2(small cell), and P1. The P1 crystal form only grows in the presence of 1.4 M DMSO; there are 2 DMSO binding sites per subunit. The structures have been refined at a resolution of 1.83 A, 2.10 A, 2.13 A, and 1.80 A, respectively. In the 4 different spacegroups the TIM subunit can be observed in the context of 7 different crystallographic environments. In the C2 cells, the dimer 2-fold axis coincides with a crystallographic 2-fold axis. The similarities and differences of the 7 subunits are discussed. In 6 subunits the flexible loop (loop 6) is open, whereas in the P2(1)2(1)2(1) cell, the flexible loop of subunit 2 is in an almost closed conformation. The crystal contacts in the 4 different crystal forms are predominantly generated by polar residues in loops. A statistical analysis of the residues involved in crystal contacts shows that, in particular, serines are frequently involved in these interactions; 19% of the exposed serines are involved in crystal contacts.