X-ray structures of three interface mutants of gammaB-crystallin from bovine eye lens.

GammaB-crystallin consists of two domains each comprising two "Greek key" motifs. Both domains fold independently, and domain interactions contribute significantly to the stability of the C-terminal domain. In a previous study (Palme S et al., 1996, Protein Sci 6:1529-1636) it was shown that Phe56 from the N-terminal domain, a residue involved in forming a hydrophobic core at the domain interface, effects the interaction of the two domains, and therefore, the stability of the C-terminal domain. Ala or Asp at position 56 drastically decreased the stability of the C-terminal domain, whereas Trp had a more moderate effect. In this article we present the X-ray structures of these interface mutants and correlate them with the stability data. The mutations do not effect the overall structure of the molecule. No structural changes are observed in the vicinity of the replaced residue, suggesting that the local structure is too rigid to allow compensations for the amino acid replacements. In the mutants gammaB-F56A and -F56D, a solvent-filled groove accessible to the bulk solvent is created by the replacement of the bulky Phe side chain. In gammaB-F56W, the pyrrole moiety of the indole ring replaces the phenyl side chain of the wild type. With the exception of gammaB-F56W, there is a good correlation between the hydrophobicity of the amino acid at position 56 according to the octanol scale and the stability of the C-terminal domain. In gammaB-F56W, the C-terminal domain is less stable than estimated from the hydrophobicity, presumably because the ring nitrogen (Nepsilon1) has no partner to form hydrogen bonds. The data suggest that the packing of hydrophobic residues in the interface core is important for domain interactions and the stability of gammaB-crystallin. Apparently, for protein stability, the same principles apply for hydrophobic cores within domains and at domain interfaces.     

Amino acid sequence of bovine gamma E (IVa) lens crystallin.

When electrospray ionization mass spectrometry (ESMS) was used to analyze purified bovine gamma E (gamma IVa)-crystallin, it yielded a relative molecular mass (M(r)) of 20.955 +/- 5. This mass is significantly different from that calculated from the published sequence (M(r) 20.894) (White HE et al., 1989, J Mol Biol 207:217-235). Further, ES-MS analysis of the protein after it had been reduced and carboxymethylated indicated the presence of five cysteine residues, whereas the published sequence contains six (Kilby GW et al., 1995, Eur Mass Spectrom 1:203-208). The entire protein sequence of gamma E crystallin has therefore been studied via a combination of ES-MS, ES-MS/MS, and Edman amino acid sequencing. The corrected sequence gives an M(r) of 20.955.3, which matches that obtained by ES-MS analysis of the purified native protein. The corrected sequence is also in agreement with a recent cDNA sequence obtained for a bovine gamma-crystallin by R. Hay (pers. comm.).     

Three-dimensional model and quaternary structure of the human eye lens protein gamma S-crystallin based on beta- and gamma-crystallin X-ray coordinates and ultracentrifugation.

A 3-dimensional model of the human eye lens protein gamma S-crystallin has been constructed using comparative modeling approaches encoded in the program COMPOSER on the basis of the 3-dimensional structure of gamma-crystallin and beta-crystallin. The model is biased toward the monomeric gamma B-crystallin, which is more similar in sequence. Bovine gamma S-crystallin was shown to be monomeric by analytical ultracentrifugation without any tendency to form assemblies up to concentrations in the millimolar range. The connecting peptide between domains was therefore built assuming an intramolecular association as in the monomeric gamma-crystallins. Because the linker has 1 extra residue compared with gamma B and beta B2, the conformation of the connecting peptide was constructed by using a fragment from a protein database. gamma S-crystallin differs from gamma B-crystallin mainly in the interface region between domains. The charged residues are generally paired, although in a different way from both beta- and gamma-crystallins, and may contribute to the different roles of these proteins in the lens.     

Surface interactions of gamma-crystallins in the crystal medium in relation to their association in the eye lens

A comparative study of intermolecular interactions in crystals of two homologous low molecular weight proteins, gamma-II and gamma-IIIb crystallins, from calf eye lens was carried out. Crystal packings for these proteins are very different: intermolecular contact areas compose about 33% of the total accessible surface area of gamma-II as compared with 13% in gamma-III. Two key residues seem to be mainly responsible for the differences in protein association in the crystal medium. These are Ser 103 and Leu 155 in gamma-II, which are replaced by Met 103 and His 155 in gamma-IIb. A similar substitution of these residues is observed in different gene products of gamma-crystallins from a number of vertebrates. This is consistent with the existence of a genetically controlled mechanism for determining intermolecular association of gamma-crystallins in the native medium of the lens.     

Electrostatic interactions in the reconstitution of an SH2 domain from constituent peptide fragments

Fragment complementation has been used to delineate the essential recognition elements for stable folding in Src homology 2 (SH2) domains by using NMR spectroscopy, alanine scanning, and surface plasmon resonance. The unfolded 9-kD and 5-kD peptide fragments formed by limited proteolytic digestion of the N-terminal SH2 domain from the p85alpha subunit of phosphatidylinositol 3'-kinase fold into an active native-like structure on interaction with one another. The corresponding 5-kD fragment of the homologous Src protein, however, was not capable of structurally complementing the p85 9-kD fragment, indicating that fragment complementation among these SH2 domains is sensitive to the sequence differences between the Src and p85 domains. Partial complementation and folding activity could be recovered with hybrid sequences of these SH2 domains. Complete alanine scanning of the 5-kD p85 fragment was used to identify the sequence recognition elements required for complex formation. The alanine substitutions in the p85 5-kD fragment that abolished binding affinity with the cognate 9-kD fragment correlate well with highly conserved residues among SH2 domains that are either integrally involved in core packing or found at the interface between fragments. Surprisingly, however, mutation of a nonconserved surface-exposed aspartic acid to alanine was found to have a significant effect on complementation. A single additional mutation of arginine to aspartic acid allowed for recovery of native structure and increased the thermal stability of the designed Src-p85 chimera by 18 degrees C. This modification appears to relieve an unfavorable surface electrostatic interaction, demonstrating the importance of surface charge interactions in protein stability.     

Lens Crystallins and Their Microbial Homologs: Structure,Stability, and Function

Referee: Franz Schmid, Biochemicshes Laboratorium, Universitaet Bayeuth, D-95440 Bayeuth, Germany

abg-Crystallins are the major protein components in the vertebrate eye lens — a as a molecular chaperone and b and g as structural proteins. Surprisingly, the latter two share some structural characteristics with a number of microbial stress proteins. The common denominator is not only the Greek key topology of their polypeptide chains but also their high intrinsic stability, which, in certain microbial crystallin homologs, is further enhanced by high-affinity Ca²⁺-binding. Recent studies of natural and mutant vertebrate bg-crystallins as well as spherulin 3a from Physarum polycephalum and Protein S from Myxococcus xanthus allowed the correlation of structure and stability of crystallins to be elucidated in some detail. From the thermo-dynamic point of view, stability increments come from (1) local interactions involved in the close packing of the cooperative units, (2) the all-b secondary structure of the Greek-key motif, (3) intramolecular interactions between domains, (4) intermolecular domain interactions, including 3D domain swapping and (v) excluded volume effects due to “molecular crowding” at the high cellular protein concentrations. Apart from these contributions to the Gibbs free energy of stability, significant kinetic stabilization originates from the high activation energy barrier determining the rate of unfolding from the native to the unfolded state. From the functional point of view, the high stability is responsible for the long-term transparency of the eye lens, on the one hand, and the stress resistance of the microorganisms in their dormant state on the other. Local structural perturbations due to chemical modification, wrong protein interactions, or other irreversible processes may lead to protein aggregation. A leading cataract hypothesis is that only after a-crystallin, a member of the small heat-shock protein family, is titrated out does pathological opacity occur. Understanding the structural basis of protein stability in the healthy eye lens is the route to solve the enormous medical and economical problem of cataract.     

Stability and interactions of the amino-terminal domain of ClpB from Escherichia coli

ClpB is a member of a multichaperone system in Escherichia coli (with DnaK, DnaJ, and GrpE) that reactivates aggregated proteins. The sequence of ClpB contains two ATP-binding regions that are enclosed between the N- and C-terminal extensions. Whereas it has been found that the N-terminal region of ClpB is essential for the chaperone activity, the structure of this region is not known, and its biochemical properties have not been studied. We expressed and purified the N-terminal fragment of ClpB (residues 1-147). Circular dichroism of the isolated N-terminal region showed a high content of alpha-helical structure. Differential scanning calorimetry showed that the N-terminal region of ClpB is thermodynamically stable and contains a single folding domain. The N-terminal domain is monomeric, as determined by gel-filtration chromatography, and the elution profile of the N-terminal domain does not change in the presence of the N-terminally truncated ClpB (ClpBDeltaN). This indicates that the N-terminal domain does not form strong contacts with ClpBDeltaN. Consistently, addition of the separated N-terminal domain does not reverse an inhibition of ATPase activity of ClpBDeltaN in the presence of casein. As shown by ELISA measurements, full-length ClpB and ClpBDeltaN bind protein substrates (casein, inactivated luciferase) with similar affinity. We also found that the isolated N-terminal domain of ClpB interacts with heat-inactivated luciferase. Taken together, our results indicate that the N-terminal fragment of ClpB forms a distinct domain that is not strongly associated with the ClpB core and is not required for ClpB interactions with other proteins, but may be involved in recognition of protein substrates.     

The relationship between conservation, thermodynamic stability, and function in the SH3 domain hydrophobic core

To investigate the relationships between sequence conservation, protein stability, and protein function, we have measured the thermodynamic stability, folding kinetics, and in vitro peptide-binding activity of a large number of single-site substitutions in the hydrophobic core of the Fyn SH3 domain. Comparison of these data to that derived from an analysis of a large alignment of SH3 domain sequences revealed a very good correlation between the distinct pattern of conservation observed at each core position and the thermodynamic stability of mutants. Conservation was also found to correlate well with the unfolding rates of mutants, but not to the folding rates, suggesting that evolution selects more strongly for optimal native state packing interactions than for maximal folding rates. Structural analysis suggests that residue-residue core packing interactions are very similar in all SH3 domains, which provides an explanation for the correlation between conservation and mutant stability effects studied in a single SH3 domain. We also demonstrate a correlation between stability and the in vivo activity of mutants, and between conservation and activity. However, the relationship between conservation and activity was very strong only for the three most conserved hydrophobic core positions. The weaker correlation between activity and conservation seen at the other seven core positions indicates that maintenance of protein stability is the dominant selective pressure at these positions. In general, the pattern of conservation at hydrophobic core positions appears to arise from conserved packing constraints, and can be effectively utilized to predict the destabilizing effects of amino acid substitutions.     

The optimization of protein-solvent interactions: thermostability and the role of hydrophobic and electrostatic interactions.

Protein-solvent interactions were analyzed using an optimization parameter based on the ratio of the solvent-accessible area in the native and the unfolded protein structure. The calculations were performed for a set of 183 nonhomologous proteins with known three-dimensional structure available in the Protein Data Bank. The dependence of the total solvent-accessible surface area on the protein molecular mass was analyzed. It was shown that there is no difference between the monomeric and oligomeric proteins with respect to the solvent-accessible area. The results also suggested that for proteins with molecular mass above some critical mass, which is about 28 kDa, a formation of domain structure or subunit aggregation into oligomers is preferred rather than a further enlargement of a single domain structure. An analysis of the optimization of both protein-solvent and charge-charge interactions was performed for 14 proteins from thermophilic organisms. The comparison of the optimization parameters calculated for proteins from thermophiles and mesophiles showed that the former are generally characterized by a high degree of optimization of the hydrophobic interactions or, in cases where the optimization of the hydrophobic interactions is not sufficiently high, by highly optimized charge-charge interactions.     

Role of interchain alpha-helical hydrophobic interactions in Ca2+ affinity, formation, and stability of a two-site domain in troponin C.

We have previously shown that a 34-residue synthetic peptide representing the calcium-binding site III of troponin C formed a symmetric two-site dimer consisting of two helix-loop-helix motifs arranged in a head-to-tail fashion (Shaw, G.S., Hodges, R.S., & Sykes, B.D., 1990, Science 249, 280-283). In this study the hydrophobicities of the alpha-helices were altered by replacing L-98 and F-102 in the N-terminal region and/or I-121 and L-122 in the C-terminal region with alanine residues. Our results showed that substitution of hydrophobic residues either in the N- or C-terminal region have little effect on alpha-helix formation but resulted in a 100- and 300-fold decrease in Ca2+ affinity, respectively. Simultaneous substitution of both hydrophobes in the N- and C-terminal region resulted in a 1,000-fold decrease in Ca2+ affinity. Data from guanidine hydrochloride denaturation studies suggested that intermolecular interactions occur and that the less hydrophobic analogs had a lower overall conformational stability. These data support the contention that the hydrophobic residues are important in the formation of the two-site domain in troponin C, and this hydrophobic association stabilizes Ca2+ affinity.     

Contributions of hydrophobic domain interface interactions to the folding and stability of human gammaD-crystallin

Human gammaD-crystallin (HgammaD-Crys) is a monomeric eye lens protein composed of two highly homologous beta-sheet domains. The domains interact through interdomain side chain contacts forming two structurally distinct regions, a central hydrophobic cluster and peripheral residues. The hydrophobic cluster contains Met43, Phe56, and Ile81 from the N-terminal domain (N-td) and Val132, Leu145, and Val170 from the C-terminal domain (C-td). Equilibrium unfolding/refolding of wild-type HgammaD-Crys in guanidine hydrochloride (GuHCl) was best fit to a three-state model with transition midpoints of 2.2 and 2.8 M GuHCl. The two transitions likely corresponded to sequential unfolding/refolding of the N-td and the C-td. Previous kinetic experiments revealed that the C-td refolds more rapidly than the N-td. We constructed alanine substitutions of the hydrophobic interface residues to analyze their roles in folding and stability. After purification from E. coli, all mutant proteins adopted a native-like structure similar to wild type. The mutants F56A, I81A, V132A, and L145A had a destabilized N-td, causing greater population of the single folded domain intermediate. Compared to wild type, these mutants also had reduced rates for productive refolding of the N-td but not the C-td. These data suggest a refolding pathway where the domain interface residues of the refolded C-td act as a nucleating center for refolding of the N-td. Specificity of domain interface interactions is likely important for preventing incorrect associations in the high protein concentrations of the lens nucleus.     

Identification of the posttranslational modifications of bovine lens alpha B-crystallins by mass spectrometry.

A combination of mass spectrometric techniques has been used to investigate the amino acid sequence and post-translational modifications of alpha B-crystallin isolated from bovine lenses by gel filtration chromatography and reversed-phase high performance liquid chromatography. Chromatographic fractions were analyzed by electrospray ionization mass spectrometry to determine the homogeneity and molecular weights of proteins in the fractions. The alpha B-crystallin primary gene product, its mono- and diphosphorylated forms, its N- and C-terminal truncated forms, as well as other lens proteins unrelated to the alpha B-crystallins were identified by their molecular weights. Detailed information about the sites of phosphorylation, as well as evidence supporting reassignment of Asn to Asp at position 80, was obtained by analyzing proteolytic digests of these proteins by fast atom bombardment mass spectrometry. Results of this investigation indicate that alpha B-crystallin is phosphorylated in vivo at Ser 45, Ser 59, and either Ser 19 or 21. From the specificity of phosphorylation of alpha-crystallins, it appears that there may be two different kinases responsible for their phosphorylation.     

RIP death domain structural interactions implicated in TNF-mediated proliferation and survival

Death domain (DD)-containing proteins are involved in both apoptosis and survival/proliferation signaling induced by activated death receptors. Here, a phylogenetic and structural analysis was performed to highlight differences in DD domains and their key regulatory interaction sites. The phylogenetic analysis shows that receptor DDs are more conserved than DDs in adaptors. Adaptor DDs can be subdivided into those that activate or inhibit apoptosis. Modeling of six homotypic DD interactions involved in the TNF signaling pathway implicates that the DD of RIP (Receptor interacting protein kinase 1) is capable of interacting with the DD of TRADD (TNFR1-associated death domain protein) in two different, exclusive ways: one that subsequently recruits CRADD (apoptosis/inflammation) and another that recruits NFkappaB (survival/proliferation).     

Heat perturbation of bovine eye lens alpha-crystallin probed by covalently attached ratiometric fluorescent dye 4'-diethylamino-3-hydroxyflavone

Bovine eye lens alpha-crystallin was covalently labeled with 6-bromomethyl-4'-diethylamino-3-hydroxyflavone and studied under native-like conditions and at the elevated temperature (60 degrees C) that is known to facilitate alpha-crystallin chaperone-like activity. This novel SH-reactive two-band ratiometric fluorescent probe is characterized by two highly emissive N*- and T*-bands; the latter appears due to excited state intramolecular proton transfer reaction. The positions of these bands and the ratio of their intensities for the alpha-crystallin-dye conjugate are the sensitive indicators of polarity of the dye environment and its participation in intermolecular hydrogen bonding. Although we found that the dye labels both the SH and the NH2 groups in alpha-crystallin, a recently developed procedure allowed us to distinguish between the heat-induced spectral changes of the dye molecules attached to SH and NH2 groups. We observed that at elevated temperature the environment of the SH-attached dye becomes more polar and flexible. The number of H-bond acceptor groups in the vicinity of the dye decreases. Since alpha-crystallin contains a single Cys residue within the C-terminal domain of its (alpha)A subunit (the (alpha)B subunit contains none), we can attribute the observed effects to temperature-induced changes in the C-terminal domain of this protein.     

Comparative electrophoretic studies of muscle and eye lens proteins in freshwater fish of Iraq

Muscle myogens and eye lens proteins have been studied in ten species of freshwater fish from Iraq. The electrophoretic analysis revealed that the muscle myogens can be considered as a good taxonomic criterion to differentiate the family Mugilidae from the Cyprinodontidae and Cyprinidae, but not between the families Poeciliidae and Cyprinodontidae. Within the Cyprinidae the muscle myogens can be used to diferentiate Barbus grypus from the remaining species of this family. Eye lens proteins are not considered a good taxonomic criterion to differentiate the members of the four families studied, but can distinguish B. belawayei and B. grypus from the other Cyprinid species.     

Comparative electrophoretic studies of muscle, eye lens and heart protein in fishes from the Arabian Gulf

Muscle myogen, eye lens and heart protein have been studied in 31 species of fishes from Kuwait, Arabian Gulf. The electrophoretic analysis revealed that muscle myogen can be considered a good taxonomic criterion to differentiate the families Ariidae, Belonidae and Lutjanidae from the other fish families studied. Within the Lutjanidae and Pomadasyidae the muscle myogens can be used to differentiate Lutjanus kasmira and Pomadasys argenteus, respectively. The muscle myogens can also be used to differentiate Sardinella perforata, Platycephalus inducus, Drepane longemana and Psethodes erumei. Eye lens protein can be considered a good taxonomic criterion to differentiate between the fish species studied and especially among the families Belonidae, Lutjanidae, Pomadsyidae, Sparidae and Scaenidae. Within the Lutjanidae and Sparidae, the eye lens protein can be used to differentiate Lutjanus kasmira and Acanthopagrus berda, respectively. Heart protein is not considered a good taxonomic criterion to differentiate the species of fishes studied, but can distinguish Lutjanus coccineus from the remaining species studied.     

Predicting protein domain interactions from coevolution of conserved regions

The knowledge of protein and domain interactions provide crucial insights into their function within a cell. Several computational methods have been proposed to detect interactions between proteins and their constitutive domains. In this work, we focus on approaches based on correlated evolution (coevolution) of sequences of interacting proteins. In this type of approach, often referred to as the mirrortree method, a high correlation of evolutionary histories of two proteins is used as an indicator to predict protein interactions. Recently, it has been observed that subtracting the underlying speciation process by separating coevolution due to common speciation divergence from that due to common function of interacting pairs greatly improves the predictive power of the mirrortree approach. In this article, we investigate possible improvements and limitations of this method. In particular, we demonstrate that the performance of the mirrortree method that can be further improved by restricting the coevolution analysis to the relatively conserved regions in the protein domain sequences (disregarding highly divergent regions). We provide a theoretical validation of our results leading to new insights into the interplay between coevolution and speciation of interacting proteins.     

Finding biologically relevant protein domain interactions: conserved binding mode analysis

Proteins evolved through the shuffling of functional domains, and therefore, the same domain can be found in different proteins and species. Interactions between such conserved domains often involve specific, well-determined binding surfaces reflecting their important biological role in a cell. To find biologically relevant interactions we developed a method of systematically comparing and classifying protein domain interactions from the structural data. As a result, a set of conserved binding modes (CBMs) was created using the atomic detail of structure alignment data and the protein domain classification of the Conserved Domain Database. A conserved binding mode is inferred when different members of interacting domain families dock in the same way, such that their structural complexes superimpose well. Such domain interactions with recurring structural themes have greater significance to be biologically relevant, unlike spurious crystal packing interactions. Consequently, this study gives lower and upper bounds on the number of different types of interacting domain pairs in the structure database on the order of 1000-2000. We use CBMs to create domain interaction networks, which highlight functionally significant connections by avoiding many infrequent links between highly connected nodes. The CBMs also constitute a library of docking templates that may be used in molecular modeling to infer the characteristics of an unknown binding surface, just as conserved domains may be used to infer the structure of an unknown protein. The method's ability to sort through and classify large numbers of putative interacting domain pairs is demonstrated on the oligomeric interactions of globins.     

Sequence-structure signals of 3D domain swapping in proteins

Three-dimensional domain swapping occurs when two or more identical proteins exchange identical parts of their structure to generate an oligomeric unit. It affects proteins with diverse sequences and structures, and is expected to play important roles in evolution, functional regulation and even conformational diseases. Here, we search for traces of domain swapping in the protein sequence, by means of algorithms that predict the structure and stability of proteins using database-derived potentials. Regions whose sequences are not optimal with regard to the stability of the native structure, or showing marked intrinsic preferences for non-native conformations in absence of tertiary interactions are detected in most domain-swapping proteins. These regions are often located in areas crucial in the swapping process and are likely to influence it on a kinetic or thermodynamic level. In addition, cation-pi interactions are frequently observed to zip up the edges of the interface between intertwined chains or to involve hinge loop residues, thereby modulating stability. We end by proposing a set of mutations altering the swapping propensities, whose experimental characterization would contribute to refine our in silico derived hypotheses.     

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.