A conserved hydrophobic core at Bcl-xL mediates its structural stability and binding affinity with BH3-domain peptide of pro-apoptotic protein

Bcl-2 family proteins regulate apoptosis through their homo- and heterodimerization. By protein sequence analysis and structural comparison, we have identified a conserved hydrophobic core at the BH1 and BH2 domains of Bcl-2 family proteins. The hydrophobic core is stabilized by hydrophobic interactions among the residues of Trp137, Ile140, Trp181, Ile182, Trp188 and Phe191 in Bcl-x_L. Destabilization of the hydrophobic core can promote the protein unfolding and pore formation in synthetic lipid vesicles. Interestingly, though the hydrophobic core does not participate in binding with BH3 domain of pro-apoptotic proteins, disruption of the hydrophobic core can reduce the affinity of Bcl-x_L with BH3-domain peptide by changing the conformation of Bcl-x_L C-terminal residues that are involved in the peptide interaction. The BH3-domain peptide binding affinity and pore forming propensity of Bcl-x_L were correlated to its death-repressor activity, which provides new information to help study the regulatory mechanism of anti-apoptotic proteins. Meanwhile, as the tryptophans are conserved in the hydrophobic core, in vitro binding assay based on FRET of “Trp → AEDANS” can be devised to screen for new modulators targeting anti-apoptotic proteins as well as “multi-BH domains” pro-apoptotic proteins.     

Conserved residues and their role in the structure, function, and stability of acyl-coenzyme A binding protein

In the family of acyl-coenzyme A binding proteins, a subset of 26 sequence sites are identical in all eukaryotes and conserved throughout evolution of the eukaryotic kingdoms. In the context of the bovine protein, the importance of these 26 sequence positions for structure, function, stability, and folding has been analyzed using single-site mutations. A total of 28 mutant proteins were analyzed which covered 17 conserved sequence positions and three nonconserved positions. As a first step, the influence of the mutations on the protein folding reaction has been probed, revealing a folding nucleus of eight hydrophobic residues formed between the N- and C-terminal helices [Kragelund, B. B., et al. (1999) Nat. Struct. Biol. (In press)]. To fully analyze the role of the conserved residues, the function and the stability have been measured for the same set of mutant proteins. Effects on function were measured by the extent of binding of the ligand dodecanoyl-CoA using isothermal titration calorimetry, and effects on protein stability were measured with chemical denaturation followed by intrinsic tryptophan and tyrosine fluorescence. The sequence sites that have been conserved for direct functional purposes have been identified. These are Phe5, Tyr28, Tyr31, Lys32, Lys54, and Tyr73. Binding site residues are mainly polar or charged residues, and together, four of these contribute approximately 8 kcal mol-1 of the total free energy of binding of 11 kcal mol-1. The sequence sites conserved for stability of the structure have likewise been identified and are Phe5, Ala9, Val12, Leu15, Leu25, Tyr28, Lys32, Gln33, Tyr73, Val77, and Leu80. Essentially, all of the conserved residues that maintain the stability are hydrophobic residues at the interface of the helices. Only one conserved polar residue, Gln33, is involved in stability. The results indicate that conservation of residues in homologous proteins may result from a summed optimization of an effective folding reaction, a stable native protein, and a fully active binding site. This is important in protein design strategies, where optimization of one of these parameters, typically function or stability, may influence any of the others markedly.     

Sequence of gene malG in E. coli K12: homologies between integral membrane components from binding protein-dependent transport systems.

The MalG protein is needed for the transport of maltose in Escherichia coli K12. We present the sequence of gene malG. The deduced amino acid sequence corresponds to a protein of 296 amino acid residues (mol. wt. = 32 188 daltons). This protein is largely hydrophobic (hydrophobic index = 0.83) and is thus presumably an integral inner membrane protein which could span the membrane through six hydrophobic segments. We provide direct evidence from fusion proteins for the translation frame and we also identified the in vitro made MalG protein. We have found a sequence which is highly conserved between MalG and MalF, the other integral inner membrane protein of the maltose transport system. This conserved sequence is also present in all known integral membrane proteins of binding protein-dependent transport systems, always at the same distance (approximately 90 residues) from their COOH terminus. We discuss briefly this finding.     

The subunit interfaces of oligomeric enzymes are conserved to a similar extent to the overall protein sequences.

It is well established that, within families of homologous enzymes, amino acid residues that are involved in the chemistry of the reaction are highly conserved. To determine if residues at the subunit interface of oligomeric enzymes with shared active sites are also conserved, comparative analysis of five enzyme families was undertaken. For the chosen enzyme families, sequence data were available for a large number of proteins and a three-dimensional structure was known for at least two members of each family. The analysis indicates that the subunit interface and the hydrophobic core of proteins from all five families have diverged to a similar extent to the overall protein sequences.     

Conserved residue clustering and protein structure prediction

Protein residues that are critical for structure and function are expected to be conserved throughout evolution. Here, we investigate the extent to which these conserved residues are clustered in three-dimensional protein structures. In 92% of the proteins in a data set of 79 proteins, the most conserved positions in multiple sequence alignments are significantly more clustered than randomly selected sets of positions. The comparison to random subsets is not necessarily appropriate, however, because the signal could be the result of differences in the amino acid composition of sets of conserved residues compared to random subsets (hydrophobic residues tend to be close together in the protein core), or differences in sequence separation of the residues in the different sets. In order to overcome these limits, we compare the degree of clustering of the conserved positions on the native structure and on alternative conformations generated by the de novo structure prediction method Rosetta. For 65% of the 79 proteins, the conserved residues are significantly more clustered in the native structure than in the alternative conformations, indicating that the clustering of conserved residues in protein structures goes beyond that expected purely from sequence locality and composition effects. The differences in the spatial distribution of conserved residues can be utilized in de novo protein structure prediction: We find that for 79% of the proteins, selection of the Rosetta generated conformations with the greatest clustering of the conserved residues significantly enriches the fraction of close-to-native structures.     

Amino acid sequence determination and three-dimensional modelling of thioredoxin from the photosynthetic bacterium Rhodobacter sphaeroides Y

The complete primary structure of thioredoxin from Rhodobacter sphaeroides Y has been determined by analysis of peptides after cleavage with cyanogen bromide, chymotrypsin and trypsin. Peptides were separated by HPLC and analyzed by liquid-phase and gas-phase sequencer degradations. The protein consists of 105 residues (Mr = 11,180); its amino acid sequence shows a clear homology to the five known thioredoxins from plant or bacterial sources, with 40-56% residue identity when the proteins are aligned at the active-site disulfide. Not only the active-site regions are conserved, but also residues which belong to the hydrophobic surface suggested to be important for binding of procaryote thioredoxins in redox interactions with other proteins (residues 75-76; 91-93 in Escherichia coli). A three-dimensional model of Rb. sphaeroides thioredoxin has been derived from the E. coli crystallographic structure with computer graphics. This model indicates that the overall structures as well as the active sites are closely similar; however, the residue substitutions allow both proteins to adopt different local folding as shown in the hydrophobic core.     

NMR structure of the conserved hypothetical protein TM0487 from Thermotoga maritima: implications for 216 homologous DUF59 proteins

The NMR structure of the conserved hypothetical protein TM0487 from Thermotoga maritima represents an alpha/beta-topology formed by the regular secondary structures alpha1-beta1-beta2-alpha2-beta3-beta4-alpha3- beta5-3(10)-alpha4, with a small anti-parallel beta-sheet of beta-strands 1 and 2, and a mixed parallel/anti-parallel beta-sheet of beta-strands 3-5. Similar folds have previously been observed in other proteins, with amino acid sequence identity as low as 3% and a variety of different functions. There are also 216 sequence homologs of TM0487, which all have the signature sequence of domains of unknown function 59 (DUF59), for which no three-dimensional structures have as yet been reported. The TM0487 structure thus presents a platform for homology modeling of this large group of DUF59 proteins. Conserved among most of the DUF59s are 13 hydrophobic residues, which are clustered in the core of TM0487. A putative active site of TM0487 consisting of residues D20, E22, L23, T51, T52, and C55 is conserved in 98 of the 216 DUF59 sequences. Asp20 is buried within the proposed active site without any compensating positive charge, which suggests that its pK(a) value may be perturbed. Furthermore, the DUF59 family includes ORFs that are part of a conserved chromosomal group of proteins predicted to be involved in Fe-S cluster metabolism.     

Repacking of hydrophobic residues in a stable mutant of apocytochrome b562 selected by phage-display and proteolysis

To test a hydrophobic core-directed protein design approach, we previously have used phage-display and proteolysis to select stably folded proteins from a library of mutants of apocytochrome b562. The consensus sequence of the selected mutants has hydrophilic residues at two of the three positions that are designed to form a hydrophobic core. To understand this unexpected result, we determined the high-resolution structure of one of the selected mutants using multi-dimensional nuclear magnetic resonance (NMR). The structure shows that the two hydrophilic residues in the consensus sequence were on the surface of the structure. Instead, two of their neighboring hydrophobic residues reorganized their side-chain conformations and formed the hydrophobic core. This result suggests that the hydrophobic core-directed protein design by phage-display and proteolysis is a valid method in general but alternative hydrophobic packing needs to be considered in the initial design. The unexpected repacking of the hydrophobic residues also highlights the plastic nature of protein structures.     

Detection of secondary structure elements in proteins by hydrophobic cluster analysis.

Hydrophobic cluster analysis (HCA) is a protein sequence comparison method based on alpha-helical representations of the sequences where the size, shape and orientation of the clusters of hydrophobic residues are primarily compared. The effectiveness of HCA has been suggested to originate from its potential ability to focus on the residues forming the hydrophobic core of globular proteins. We have addressed the robustness of the bidimensional representation used for HCA in its ability to detect the regular secondary structure elements of proteins. Various parameters have been studied such as those governing cluster size and limits, the hydrophobic residues constituting the clusters as well as the potential shift of the cluster positions with respect to the position of the regular secondary structure elements. The following results have been found to support the alpha-helical bidimensional representation used in HCA: (i) there is a positive correlation (clearly above background noise) between the hydrophobic clusters and the regular secondary structure elements in proteins; (ii) the hydrophobic clusters are centred on the regular secondary structure elements; (iii) the pitch of the helical representation which gives the best correspondence is that of an alpha-helix. The correspondence between hydrophobic clusters and regular secondary structure elements suggests a way to implement variable gap penalties during the automatic alignment of protein sequences.     

De novo design of the hydrophobic cores of proteins.

We have developed and experimentally tested a novel computational approach for the de novo design of hydrophobic cores. A pair of computer programs has been written, the first of which creates a "custom" rotamer library for potential hydrophobic residues, based on the backbone structure of the protein of interest. The second program uses a genetic algorithm to globally optimize for a low energy core sequence and structure, using the custom rotamer library as input. Success of the programs in predicting the sequences of native proteins indicates that they should be effective tools for protein design. Using these programs, we have designed and engineered several variants of the phage 434 cro protein, containing five, seven, or eight sequence changes in the hydrophobic core. As controls, we have produced a variant consisting of a randomly generated core with six sequence changes but equal volume relative to the native core and a variant with a "minimalist" core containing predominantly leucine residues. Two of the designs, including one with eight core sequence changes, have thermal stabilities comparable to the native protein, whereas the third design and the minimalist protein are significantly destabilized. The randomly designed control is completely unfolded under equivalent conditions. These results suggest that rational de novo design of hydrophobic cores is feasible, and stress the importance of specific packing interactions for the stability of proteins. A surprising aspect of the results is that all of the variants display highly cooperative thermal denaturation curves and reasonably dispersed NMR spectra. This suggests that the non-core residues of a protein play a significant role in determining the uniqueness of the folded structure.     

Conservation Throughout Vertebrate Evolution of the Predicted β-Strands in Myelin Basic Protein

To identify functionally important parts of the 18.5-kDa myelin basic protein (MBP), the amino acid sequences from 10 species ranging from shark to human were aligned using the SEQHP computer program. The residues that are invariant or very conservatively substituted (Arg/Lys, Ser/Thr, Ile/Leu, Asp/Glu) among all 10 proteins were scored. Of the 72 conserved residues in the 170-residue human protein (42% conserved), 32 are found within the five beta-strands previously predicted (45 residues, 71% conserved), 23 within the small-loops region (42 residues, 55% conserved), but only 17 within the large-loops region (83 residues, 20% conserved). Of the 22 hydrophobic residues within the predicted beta-sheet of human MBP, 20 hydrophobic residues remain in the shark protein, 19 of them in the same positions. In contrast, there are 10 hydrophobic residues elsewhere in the human protein, but only 7 remain in the shark protein and only 1 of them is in the same position. The triprolyl sequence found in all mammalian MBPs and in the chicken MBP is not conserved in the shark protein. The four alternately spliced forms of mouse MBP can be accommodated by the beta-structural model, but not the 17-kDa human MBP, which lacks exon 5. These findings confirm the crucial role of the hydrophobic residues in the predicted beta-sheet for the structure and function of the protein. It seems likely that the conserved portions of the protein make an important contribution to the highly ordered lamellar structure of myelin.     

Amino acid and cDNA sequences of a methionine-rich 2S protein from sunflower seed (Helianthus annuus L.)

The amino acid sequence of a methionine-rich 2S seed protein from sunflower (Helianthus annuus L.) and the sequence of a cDNA clone which codes for the entire primary translation product have been determined. The mature protein consists of a single polypeptide chain of 103 amino acids (molecular mass 12133 Da) which contains 16 residues of methionine and 8 residues of cysteine. The cDNA sequence established that the protein is synthesized as a precursor of 141 residues with a typical hydrophobic signal sequence of 25 residues followed by a further 13-residue hydrophobic pro-sequence which is presumably removed by post-translational cleavage. The sequence of the mature protein and that deduced from the cDNA were identical with no evidence of processing at the C-terminus. Comparison of the sunflower methionine-rich protein sequence with sequences of other seed 2S proteins from dicotyledons and monocotyledons showed limited but distinct sequence similarities; in particular the arrangement of the cysteine residues was conserved. The sunflower protein shows 34% identity with the methionine-rich Brazil nut 2S protein and the prepro regions of the precursors of these two proteins show about 50% identity. This similarity indicates that these methionine-rich 2S proteins have diverged as a subclass of the 2S superfamily of proteins which contain only 2-3% methionine. While the related 2S proteins from other dicotyledons are processed to a small and large subunit, the sunflower protein is not cleaved in this way.     

A new class of hexahelical insect proteins revealed as putative carriers of small hydrophobic ligands

BACKGROUND: THP12 is an abundant and extraordinarily hydrophilic hemolymph protein from the mealworm Tenebrio molitor and belongs to a group of small insect proteins with four highly conserved cysteine residues. Despite their sequence homology to odorant-binding proteins and pheromone-binding proteins, the function of these proteins is unclear. RESULTS: The first three-dimensional structure of THP12 has been determined by multidimensional NMR spectroscopy. The protein has a nonbundle helical structure consisting of six alpha helices. The arrangement of the alpha helices has a 'baseball glove' shape. In addition to the hydrophobic core, electrostatic interactions make contributions to the overall stability of the protein. NMR binding studies demonstrated the binding of small hydrophobic ligands to the single hydrophobic groove in THP12. Comparing the structure of THP12 with the predicted secondary structure of homologs reveals a common fold for this new class of insect proteins. A search with the program DALI revealed extensive similarity between the three-dimensional structure of THP12 and the N-terminal domain (residues 1-95) of recoverin, a member of the family of calcium-binding EF-hand proteins. CONCLUSIONS: Although the biological function of this new class of proteins is as yet undetermined, a general role as alpha-helical carrier proteins for small hydrophobic ligands, such as fatty acids or pheromones, is proposed on the basis of NMR-shift perturbation spectroscopy.     

Self peptides bound by HLA class I molecules are derived from highly conserved regions of a set of evolutionarily conserved proteins

An evolutionary analysis of self peptides reported to be bound by HLA class I molecules showed that these peptides are largely derived from proteins that have been highly conserved in the history of mammals. These proteins also often have universal tissue expression and have a higher than average frequency of highly hydrophilic residues. The peptides themselves are generally still more highly conserved than the source proteins and have a higher frequency of highly hydrophobic residues, evidently often derived from conserved hydrophobic cores of the source proteins. These results suggest that the mechanism by which peptides are derived for MHC presentation may preferentially select peptides from conserved protein regions. In the case of parasite-derived peptides, such a mechanism would be adaptive in that it would reduce the likelihood of escape mutants.     

Predicted structure for the calcium-dependent membrane-binding proteins p35, p36, and p32

A new family of proteins (annexins) that bind to membranes at micromolar free Ca2+ has been recognized. Its members include an EGF-receptor kinase substrate (p35), a retroviral tyrosine kinase substrate (p36), the liver protein endonexin (p32) and an electric ray protein, calelectrin. Each protein contains four sequence repeats with a further 2-fold internal homology. Using the predicted secondary structure and pattern of conserved hydrophobic residues in each repeat, we have built a three-dimensional model that is largely isostructural with the known molecular conformation of bovine intestinal calcium-binding protein. The final (energy-refined) model had a core formed from the conserved hydrophobic residues. It differed from ICaBP principally in the length of the two Ca2+-binding loops with only one loop being able to bind. The model suggests a mechanism for interaction of these new Ca2+-binding proteins with phospholipid bilayers.     

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

Mutational analysis of the Pyrococcus furiosus holliday junction resolvase hjc revealed functionally important residues for dimer formation, junction DNA binding, and cleavage activities

The Holliday junction cleavage protein, Hjc resolvase of Pyrococcus furiosus, is the first Holliday junction resolvase to be discovered in Archaea. Although the archaeal resolvase shares certain biochemical properties with other non-archaeal junction resolvases, no amino acid sequence similarity has been identified. To investigate the structure-function relationship of this new Holliday junction resolvase, we constructed a series of mutant hjc genes using site-directed mutagenesis targeted at the residues conserved among the archaeal orthologs. The products of these mutant genes were purified to homogeneity. With analysis of the activity of the mutant proteins to bind and cleave synthetic Holliday junctions, one acidic residue, Glu-9, and two basic residues, Arg-10 and Arg-25, were found to play critical roles in enzyme action. This is in addition to the three conserved residues, Asp-33, Glu-46, and Lys-48, which are also conserved in the motif found in the type II restriction endonuclease family proteins. Two aromatic residues, Phe-68 and Phe-72, are important for the formation of the homodimer probably through hydrophobic interactions. The results of these studies have provided insights into the structure-function relationships of the archaeal Holliday junction resolvase as well as the universality and diversity of the Holliday junction cleavage reaction.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.     

Structure of sterol carrier protein 2 at 1.8 A resolution reveals a hydrophobic tunnel suitable for lipid binding

Sterol carrier protein 2, also known as nonspecific lipid transfer protein is a ubiquitous, small, basic protein of 13 kDa found in animals. Its primary structure is highly conserved between different species, and it has been implicated in the intracellular transport of lipids and in a wide range of other in vitro functions related to sterol and fatty acid metabolism. Sterol carrier protein 2 deficiency in mice leads to elevated concentrations of phytanic acid in the serum and causes hepatocarcinogenesis. However, its actual physiological role is still unknown. Although sterol carrier protein 2 has been studied extensively in the past 20 years, very little is known concerning its three-dimensional structure. The crystal structure of rabbit sterol carrier protein 2, determined at 1.8 A resolution with the MIRAS method, shows a unique alpha/beta-fold. The core of this protein forms a five-stranded antiparallel beta-sheet flanked by five helices. A C-terminal segment (residues 114-123), together with part of the beta-sheet and four alpha-helices, form a hydrophobic tunnel providing the environment for apolar ligands such as fatty acids and fatty acyl-coenzyme As. Structurally well-characterized nonspecific lipid transfer proteins from plants have hydrophobic tunnel-like cavities, which were identified as the binding site for fatty acids and related apolar ligands. Despite the fact that plant nonspecific lipid transfer proteins are smaller proteins than sterol carrier protein 2, show no sequence homology to sterol carrier protein 2, and are structurally unrelated, the cavities of these two classes of proteins are very similar with respect to size, shape, and hydrophobicity, suggesting a common functional role.