Analysis of conserved hydrophobic cores in proteins and supramolecular complexes

The conserved hydrophobic core is an important feature of a family of protein domains. We suggest a procedure for finding and the analysis of conserved hydrophobic cores. The procedure is based on using an original program called CluD (http://monkey.belozersky.msu.ru/CluD/cgi-bin/hftri.pl). Conserved hydrophobic cores of several families including homeodomains and interlock-containing domains are described. Hydrophobic clusters on some protein-DNA and protein-protein interfaces were also analyzed.     

Hydrophobic core in domains of immunoglobulin-like fold

This work analyzes proteins which contain an immunoglobulin fold, focusing on their hydrophobic core structure. The “fuzzy oil drop” model was used to measure the regularity of hydrophobicity distribution in globular domains belonging to proteins which exhibit the above-mentioned fold. Light-chain IgG domains are found to frequently contain regular hydrophobic cores, unlike the corresponding heavy-chain domains. Enzymes and DNA binding proteins present in the data-set are found to exhibit poor accordance with the hydrophobic core model.     

The immunoglobulin fold family: sequence analysis and 3D structure comparisons.

Fifty-two 3D structures of Ig-like domains covering the immunoglobulin fold family (IgFF) were compared and classified according to the conservation of their secondary structures. Members of the IgFF are distantly related proteins or evolutionarily unrelated proteins with a similar fold, the Ig fold. In this paper, a multiple structural alignment of the conserved common core is described and the correlation between corresponding sequences is discussed. While the members of the IgFF exhibit wide heterogeneity in terms of tissue and species distribution or functional implications, the 3D structures of these domains are far more conserved than their sequences. We define topologically equivalent residues in the Ig-like domains, describe the hydrophobic common cores and discuss the presence of additional strands. The disulfide bridges, not necessary for the stability of the Ig fold, may have an effect on the compactness of the domains. Based upon sequence and structure analysis, we propose the introduction of two new subtypes (C3 and C4) to the previous classifications, in addition to a new global structural classification. The very low mean sequence identity between subgroups of the IgFF suggests the occurrence of both divergent and convergent evolutionary processes, explaining the wide diversity of the superfamily. Finally, this review suggest that hydrophobic residues constituting the common hydrophobic cores are important clues to explain how highly divergent sequences can adopt a similar fold.     

Identification of compact, hydrophobically stabilized domains and modules containing multiple peptide chains.

Compactness has been used to locate discontinuous structural units containing one or more polypeptide chains in proteins of known structure. Rather than exhaustively calculating the compactness of all possible units, our procedure uses a screening algorithm to find discontinuous regions that are potentially compact. Precise calculations of compactness are restricted only to units in these regions. With our procedure, compactness can be used to discover discontinuous domains with virtually any number of disjoint peptides. Small, single-domain proteins may contain several compact regions: thus, compact regions do not always correspond to folding domains. Because a domain is an independent folding unit and should contain a hydrophobic core, compact units were further examined for the presence of hydrophobic clusters (Zehfus MH, 1995, Protein Sci 4:1188-1202). This added constraint limits the number of acceptable units and helps greatly in the location of the true structural domains. The larger hydrophobically stabilized compact units correspond to domains, while the smaller units may correspond to folding intermediates.     

A procedure for the automatic determination of hydrophobic cores in protein structures.

An algorithm is described for automatically detecting hydrophobic cores in proteins of known structure. Three pieces of information are considered in order to achieve this goal. These are: secondary structure, side-chain accessibility, and side-chain-side-chain contacts. Residues are considered to contribute to a core when they occur in regular secondary structure and have buried side chains that form predominantly nonpolar contacts with one another. This paper describes the algorithm's application to families of proteins with conserved topologies but low sequence similarities. The aim of this investigation is to determine the efficacy of the algorithm as well as to study the extent to which similar cores are identified within a common topology.     

A new protein folding algorithm based on hydrophobic compactness: Rigid Unconnected Secondary Structure Iterative Assembly (RUSSIA). II: Applications

The RUSSIA procedure (Rigid Unconnected Secondary Structure Iterative Assembly) produces structural models of cores of small- and medium-sized proteins. Loops are omitted from this treatment and regular secondary structures are reduced to points, the centers of their hydrophobic faces. This methodology relies on the maximum compactness of the hydrophobic residues, as described in detail in Part I. Starting data are the sequence and the predicted limits and natures of regular secondary structures (alpha or beta). Helices are treated as rigid cylinders, whereas beta-strands are collectively taken into account within beta-sheets modeled by helicoid surfaces. Strands are allowed to shift along their mean axis to allow some flexibility and the alpha-helices can be placed on either side of beta-sheets. Numerous initial conformations are produced by discrete rotations of the helices and sheets around the direction going from the center of their hydrophobic face to the global center of the protein. Selection of proposed models is based upon a criterion lying on the minimization of distances separating hydrophobic residues belonging to different regular secondary structures. The procedure is rapid and appears to be robust relative to the quality of starting data (nature and length of regular secondary structures). This dependence of the quality of the model on secondary structure prediction and in particular the beta-sheet topology, is one of the limits of the present algorithm. We present here some results for a set of 12 proteins (alpha, beta and alpha/beta classes) of lengths 40-166 amino acids. The r.m.s. deviations for core models with respect to the native proteins are in the range 1.4-3.7 A.     

A novel method for predicting transmembrane segments in proteins based on a statistical analysis of the SwissProt database: the PRED-TMR algorithm.

We present a novel method that predicts transmembrane domains in proteins using solely information contained in the sequence itself. The PRED-TMR algorithm described, refines a standard hydrophobicity analysis with a detection of potential termini ('edges', starts and ends) of transmembrane regions. This allows one both to discard highly hydrophobic regions not delimited by clear start and end configurations and to confirm putative transmembrane segments not distinguishable by their hydrophobic composition. The accuracy obtained on a test set of 101 non-homologous transmembrane proteins with reliable topologies compares well with that of other popular existing methods. Only a slight decrease in prediction accuracy was observed when the algorithm was applied to all transmembrane proteins of the SwissProt database (release 35). A WWW server running the PRED-TMR algorithm is available at http://o2.db.uoa. gr/PRED-TMR/     

Enrichment of integral membrane proteins for proteomic analysis using liquid chromatography-tandem mass spectrometry

An increasing number of proteomic strategies rely on liquid chromatography-tandem mass spectrometry (LC-MS/MS) to detect and identify constituent peptides of enzymatically digested proteins obtained from various organisms and cell types. However, sample preparation methods for isolating membrane proteins typically involve the use of detergents and chaotropes that often interfere with chromatographic separation and/or electrospray ionization. To address this problem, a sample preparation method combining carbonate extraction, surfactant-free organic solvent-assisted solubilization, and proteolysis was developed and demonstrated to target the membrane subproteome of Deinococcus radiodurans. Out of 503 proteins identified, 135 were recognized as hydrophobic on the basis of their calculated hydropathy values (GRAVY index), corresponding to coverage of 15% of the predicted hydrophobic proteome. Using the PSORT algorithm, 53 of the proteins identified were classified as integral outer membrane proteins and 215 were classified as integral cytoplasmic membrane proteins. All identified integral cytoplasmic membrane proteins had from 1 to 16 mapped transmembrane domains (TMDs), and 65% of those containing four or more mapped TMDs were identified by at least one hydrophobic membrane spanning peptide. The extensive coverage of the membrane subproteome (24%) by identification of highly hydrophobic proteins containing multiple TMDs validates the efficacy of the described sample preparation technique to isolate and solubilize hydrophobic integral membrane proteins from complex protein mixtures.     

An in-gel digestion procedure that facilitates the identification of highly hydrophobic proteins by electrospray ionization-mass spectrometry analysis

A procedure is described for in-gel tryptic digestion of proteins that allows the direct analysis of eluted peptides in electrospray ionization (ESI) mass spectrometers without the need of a postdigestion desalting step. It is based on the following principles: (a) a thorough desalting of the protein in-gel before digestion that takes advantage of the excellent properties of acrylamide polymers for size exclusion separations, (b) exploiting the activity of trypsin in water, in the absence of inorganic buffers, and (c) a procedure for peptide extraction using solvents of proven efficacy with highly hydrophobic peptides. Quality of spectra and sequence coverage are equivalent to those obtained after digestion in ammonium bicarbonate for hydrophilic proteins detected with Coomassie blue, mass spectrometry-compatible silver or imidazole-zinc but are significantly superior for highly hydrophobic proteins, such as membrane proteins with several transmembrane domains. ATPase subunit 9 (GRAVY 1.446) is a membrane protein channel, lipid-binding protein for which both the conventional in-gel digestion protocol and in solution digestion failed. It was identified with very high sequence coverage. Sample handling after digestion is notably simplified as peptides are directly loaded into the ESI source without postdigestion processing, increasing the chances for the identification of hydrophobic peptides.     

The crystal and molecular structure of human annexin V, an anticoagulant protein that binds to calcium and membranes

Human annexin V (PP4), a member of the family of calcium, membrane binding proteins, has been crystallized in the presence of calcium and analysed by crystallography by multiple isomorphic replacement at 3 A and preliminarily refined at 2.5 A resolution. The molecule has dimensions of 64 x 40 x 30 A3 and is folded into four domains of similar structure. Each domain consists of five alpha-helices wound into a right-handed superhelix yielding a globular structure of approximately 18 A diameter. The domains have hydrophobic cores whose amino acid sequences are conserved between the domains and within the annexin family of proteins. The four domains are folded into an almost planar array by tight (hydrophobic) pair-wise packing of domains II and III and I and IV to generate modules (II-III) and (I-IV), respectively. The assembly is symmetric with three parallel approximate diads relating II to III, I to IV and the module (II-III) to (I-IV), respectively. The latter diad marks a channel through the centre of the molecule coated with charged amino acid residues. The protein has structural features of channel forming membrane proteins and a polar surface characteristic of soluble proteins. It is a member of the third class of amphipathic proteins different from soluble and membrane proteins.     

Replacement of staphylococcal nuclease hydrophobic core residues with those from thermophilic homologues indicates packing is improved in some thermostable proteins

The importance of tight hydrophobic core packing in stabilizing proteins found in thermophilic organisms has been vigorously disputed. Here, portions of the cores found in three thermophilic homologues were transplanted into the core of staphylococcal nuclease, a protein of modest stability. Packing of the core was evaluated by comparing interaction energy of the three mutants to the comprehensive mutant library built up previously at these same sites in staphylococcal nuclease. It was found that the interaction energy of one thermophilic sequence is extraordinarily favorable and the interaction energies of other two transplanted thermophilic sequences are good, comparable to the interaction energies of mutant cores based on cores found in mesophilic homologues. As expected when transferring just a portion of the core sequence, the mutant proteins were destabilized overall relative to wild-type staphylococcal nuclease. The overall conclusion is that improvement of packing interactions is a mechanism to confer stability employed in some proteins from thermophiles, but not all.     

Protein structure alignment using environmental profiles

A new protein structure alignment procedure is described. An initial alignment is made by comparing a one-dimensional list of primary, secondary and tertiary structural features (profiles) of two proteins, without explicitly considering the three-dimensional geometry of the structures. The alignment is then iteratively refined in the second step, in which new alignments are found by three-dimensional superposition of the structures based on the current alignment. This new procedure is fast enough to do all-against-all structural comparisons routinely. The procedure sometimes finds an alignment that suggests an evolutionary relationship and which is not normally obtained if only geometry is considered. All pair-wise comparisons were made among 3539 protein structural domains that represent all known protein structures. The resulting 3539 z-scores were used to cluster the proteins. The number of main clusters increased continuously as the z-cutoff was raised, but the number of multiple-member clusters showed a maximum at z-cutoff values of 5.0 and 5.5. When a z-cutoff value of 5.0 was used, the total number of main clusters was 2043, of which only 336 clusters had more than one member.     

Specificity of the interaction between ubiquitin-associated domains and ubiquitin

Ubiquitin-associated (UBA) domains are found in a large number of proteins with diverse functions involved in ubiquitination, DNA repair, and signaling pathways. Recent studies have shown that several UBA domain proteins interact with ubiquitin (Ub), specifically p62, the phosphotyrosine-independent ligand of the SH2 domain of p56(lck); HHR23A, a human nucleotide excision repair protein; and DDI1, another damage-inducible protein. NMR chemical shift mapping reveals that Ub binds specifically but weakly to a conserved hydrophobic epitope on HHR23A UBA(1) and UBA(2) and that the UBA domains bind on the hydrophobic patch on the surface of the five-stranded beta-sheet of Ub. Models of the UBA(1)-Ub and UBA(2)-Ub complexes obtained from de novo docking reveal different orientations of the UBA domains on the Ub surface compared with those obtained by homology modeling with the related CUE domains, which also bind Ub. Our results suggest that UBA domains may interact with Ub as well as other proteins in more than one way while utilizing the same binding surface.     

Factors involved in the stability of isolated beta-sheets: Turn sequence, beta-sheet twisting, and hydrophobic surface burial

We have recently reported on the design of a 20-residue peptide able to form a significant population of a three-stranded up-and-down antiparallel beta-sheet in aqueous solution. To improve our beta-sheet model in terms of the folded population, we have modified the sequences of the two 2-residue turns by introducing the segment DPro-Gly, a sequence shown to lead to more rigid type II' beta-turns. The analysis of several NMR parameters, NOE data, as well as Deltadelta(CalphaH), DeltadeltaC(beta), and Deltadelta(Cbeta) values, demonstrates that the new peptide forms a beta-sheet structure in aqueous solution more stable than the original one, whereas the substitution of the DPro residues by LPro leads to a random coil peptide. This agrees with previous results on beta-hairpin-forming peptides showing the essential role of the turn sequence for beta-hairpin folding. The well-defined beta-sheet motif calculated for the new designed peptide (pair-wise RMSD for backbone atoms is 0.5 +/- 0.1 A) displays a high degree of twist. This twist likely contributes to stability, as a more hydrophobic surface is buried in the twisted beta-sheet than in a flatter one. The twist observed in the up-and-down antiparallel beta-sheet motifs of most proteins is less pronounced than in our designed peptide, except for the WW domains. The additional hydrophobic surface burial provided by beta-sheet twisting relative to a "flat" beta-sheet is probably more important for structure stability in peptides and small proteins like the WW domains than in larger proteins for which there exists a significant contribution to stability arising from their extensive hydrophobic cores.     

Hydrophobic Regions in Myelin Proteins Characterized Through Analysis of "Hydropathic"Profiles

Computer-generated "hydropathic" profiles were constructed for graphic comparison of the amino acid sequences for P2 protein, 18.5 kilodalton (kDa) myelin basic protein (BP), and myelin proteolipid protein (PLP). Profiles were also obtained for cytochrome b5, a membrane protein known to be capable of reversible association with lipid bilayers and of a size comparable to that of the myelin BPs. Analysis of the PLP sequence produced profiles generally compatible with the suggestions that PLP has three transbilayer and two bilayer intercalating segments. Profiles for P2 and 18.5 kDa BP were found to contain hydrophilic segments separated by relatively short hydrophobic regions. Whereas hydropathic indices in hydrophobic regions of P2, 18.5 kDa BP, and PLP fall in the value ranges recently reported for cores of globular proteins and intrabilayer domains of membrane proteins, hydrophobic sections of P2 and 18.5 kDa BP have hydropathic indices similar to those in the hydrophobic core (transprotein) regions of globular proteins. None of them are comparable to the region of cytochrome b5 known to anchor that protein in its membrane or to the segments of PLP sequence proposed as intrabilayer domains. This comparison suggests that neither BP has structural characteristics compatible with insertion into the hydrocarbon core of the myelin lipid bilayer, a conclusion that is consistent with a recently published study that identified the bilayer penetrating proteins of myelin with a hydrophobic probe. The above findings suggest an enhancement for some details of myelin architecture and a cautious approach to interpreting data for BP intercalation into bilayers.     

Hierarchic organization of domains in globular proteins

An automatic procedure is developed for the identification of domains in globular proteins from X-ray elucidated co-ordinates. Using this tool, domains are shown to be iteratively decomposable into subdomains, leading to a hierarchic molecular architecture.There is no convenient geometry that will fully characterize the atom by atom interdigitation at an interface between domains, and the strategy adopted here was devised to reduce this unwieldy three-dimensional problem to a closely approximating companion analysis in a plane. These analytically derived domain choices can be used subsequently to construct computer-generated, space-filling, color-coded views of the domains; and when this is done, the derived domains are seen to be completely resolved.The number of domains in a protein is a mathematically well-behaved function of the chain length, lending support to the supposition that the domains are an implicit structural consequence of the folding process. A spectrum of domains ranging in size from whole protein monomers to the individual units of secondary structure is apparent in each of the 22 proteins analyzed here.The hierarchic organization of structural domains is evidence in favor of an underlying protein folding process that proceeds by hierarchic condensation. In this highly constrained model, every pathway leading to the native state can be described by a tree of local folding interactions.     

Redesigning the hydrophobic core of a four-helix-bundle protein.

Rationally redesigned variants of the 4-helix-bundle protein Rop are described. The novel proteins have simplified, repacked, hydrophobic cores and yet reproduce the structure and native-like physical properties of the wild-type protein. The repacked proteins have been characterized thermodynamically and their equilibrium and kinetic thermal and chemical unfolding properties are compared with those of wild-type Rop. The equilibrium stability of the repacked proteins to thermal denaturation is enhanced relative to that of the wild-type protein. The rate of chemically induced folding and unfolding of wild-type Rop is extremely slow when compared with other small proteins. Interestingly, although the repacked proteins are more thermally stable than the wild type, their rates of chemically induced folding and unfolding are greatly increased in comparison to wild type. Perhaps as a consequence of this, their equilibrium stabilities to chemical denaturants are slightly reduced in comparison to the wild type.     

Molecular cloning of the CD9 antigen. A new family of cell surface proteins.

The CD9 antigen was described originally as a 24-kDa surface protein of non-T acute lymphoblastic leukemia cells and developing B-lymphocytes. It is also strongly expressed on platelets, among other cells, where it shows the property of mediating platelet activation and aggregation upon binding with mAbs. The primary structure has been elucidated by cloning the cDNA from a lambda gt11 expression vector library constructed with megakaryocytic mRNA. Monoclonal antibodies were used as probes with an APAAP amplification of the signal. The 5' region was further cloned in a lambda gt10 randomly primed cDNA library. The initiation codon was immediately followed by a sequence coding for the tetrapeptide corresponding to the NH2-terminal sequence identified in a microsequencing procedure. Only one species of mRNA was found with an estimated size of 1.4 kilobase. CD9 antigen appears to be a 227-amino acid molecule with four hydrophobic domains and one N-glycosylation site. Sequence and structural comparisons showed extensive similarity of the CD9 antigen with a 237-amino acid molecule described previously as the human melanoma-associated antigen ME491 and a Schistosoma mansoni membrane protein of 218 amino acids. These three proteins identify a new family of cell-surface proteins.     

Improved gel electrophoresis matrix for hydrophobic protein separation and identification

We propose an improved acrylamide gel for the separation of hydrophobic proteins. The separation strategy is based on the incorporation of N-alkylated and N,N′-dialkylated acrylamide monomers in the gel composition in order to increase hydrophobic interactions between the gel matrix and the membrane proteins. Focusing on the most efficient monomer, N,N′-dimethylacrylamide, the potentiality of the new matrix was evaluated on membrane proteins of the human colon HCT-116 cell line. Protein analysis was performed using an adapted analytical strategy based on FT-ICR tandem mass spectrometry. As a result of this comparative study, including advanced reproducibility experiments, more hydrophobic proteins were identified in the new gel (average GRAVY: −0.085) than in the classical gel (average GRAVY: −0.411). Highly hydrophobic peptides were identified reaching a GRAVY value up to 1.450, therefore indicating their probable locations in the membrane. Focusing on predicted transmembrane domains, it can be pointed out that 27 proteins were identified in the hydrophobic gel containing up to 11 transmembrane domains; in the classical gel, only 5 proteins containing 1 transmembrane domain were successfully identified. For example, multiple ionic channels and receptors were characterized in the hydrophobic gel such as the sodium/potassium channel and the glutamate or the transferrin receptors whereas they are traditionally detected using specific enrichment techniques such as immunoprecipitation. In total, membrane proteins identified in the classical gel are well documented in the literature, while most of the membrane proteins only identified on the hydrophobic gel have rarely or never been described using a proteomic-based approach.     

Single amino acids in the lumenal loop domain influence the stability of the major light-harvesting chlorophyll a/b complex

The major light-harvesting complex of photosystem II (LHCIIb) is one of the most abundant integral membrane proteins. It greatly enhances the efficiency of photosynthesis in green plants by binding a large number of accessory pigments that absorb light energy and conduct it toward the photosynthetic reaction centers. Most of these pigments are associated with the three transmembrane and one amphiphilic alpha helices of the protein. Less is known about the significance of the loop domains connecting the alpha helices for pigment binding. Therefore, we randomly exchanged single amino acids in the lumenal loop domain of the bacterially expressed apoprotein Lhcb1 and then reconstituted the mutant protein with pigments in vitro. The resulting collection of mutated recombinant LHCIIb versions was screened by using a 96-well-format plate-based procedure described previously [Heinemann, B., and Paulsen, H. (1999) Biochemistry 38, 14088-14093], enabling us to test several thousand mutants for their ability to form stable pigment-protein complexes in vitro. At least one-third of the positions in the loop domain turned out to be sensitive targets; i.e., their exchange abolished formation of LHCIIb in vitro. This confirms our earlier notion that the LHCIIb loop domains contribute more specifically to complex formation and/or stabilization than by merely connecting the alpha helices. Among the target sites, glycines and hydrophilic amino acids are more prominently represented than hydrophobic ones. Specifically, the exchange of any of the three acidic amino acids in the lumenal loop abolishes reconstitution of stable pigment-protein complexes, suggesting that ionic interactions with other protein domains are important for correct protein folding or complex stabilization. One hydrophobic amino acid, tryptophan in position 97, has been hit repeatedly in independent mutation experiments. From the LHCIIb structure and previous mutational analyses, we propose a stabilizing interaction between this amino acid and F195 near the C-proximal end of the third transmembrane helix.