Variation in structural location and amino acid conservation of functional sites in protein domain families

Background  

The functional sites of a protein present important information for determining its cellular function and are fundamental in drug design. Accordingly, accurate methods for the prediction of functional sites are of immense value. Most available methods are based on a set of homologous sequences and structural or evolutionary information, and assume that functional sites are more conserved than the average. In the analysis presented here, we have investigated the conservation of location and type of amino acids at functional sites, and compared the behaviour of functional sites between different protein domains.     

Degradation of protein disulphide bonds in dilute alkali.

The degradation of S--S bonds in 0.2 M-NaOH at 25 degrees C was studied for a series of proteins and simple aliphatic disulphide compounds, by using cathodic stripping voltammetry, ion-selective-electrode potentiometry, spectrophotometry and ultrafiltration. The disulphide bonds that dissociated in 0.2 M-NaOH were usually those that are solvent accessible and that can be reduced by mild chemical reductants. Some unexpected differences were found between similar proteins, both in the number of S--S bonds dissociated and in their rates of decomposition. Chymotrypsin has one S--S bond attacked, whereas chymotrypsinogen and trypsinogen have two. Ribonuclease A has two S--S bonds dissociated, but ribonuclease S and S-protein have three. Denaturation in 6 M-guanidine hydrochloride before alkaline digestion caused the loss of an additional S--S bond in ribonuclease A and insulin, and increased the rate of dissociation of the S--S bonds of some other proteins. The initial product of S--S bond dissociation in dilute alkali is believed to be a persulphide intermediate formed by a beta-elimination reaction. This intermediate is in mobile equilibrium with bisulphide ion, HS-, and decomposes at a mercury electrode or in acid solution to yield a stoichiometric amount of sulphide. Rate constants and equilibrium constants were measured for the equilibria between HS- and the intermediates involved in the alkaline dissociation of several proteins. Elemental sulphur was not detected in any of the protein digests. It is suggested that formation of HS- from a persulphide intermediate involves a hydrolysis reaction to yield a sulphenic acid derivative. The small polypeptides glutathione and oxytocin gave only a low yield of persulphide, and their alkaline decomposition must proceed by a mechanism different from that of the proteins.     

Using intramolecular disulfide bonds in tau protein to deduce structural features of aggregation-resistant conformations

Because tau aggregation likely plays a role in a number of neurodegenerative diseases, understanding the processes that affect tau aggregation is of considerable importance. One factor that has been shown to influence the aggregation propensity is the oxidation state of the protein itself. Tau protein, which contains two naturally occurring cysteine residues, can form both intermolecular disulfide bonds and intramolecular disulfide bonds. Several studies suggest that intermolecular disulfide bonds can promote tau aggregation in vitro. By contrast, although there are data to suggest that intramolecular disulfide bond formation retards tau aggregation in vitro, the precise mechanism underlying this observation remains unclear. While it has been hypothesized that a single intramolecular disulfide bond in tau leads to compact conformations that cannot form extended structure consistent with tau fibrils, there are few data to support this conjecture. In the present study we generate oxidized forms of the truncation mutant, K18, which contains all four microtubule binding repeats, and isolate the monomeric fraction, which corresponds to K18 monomers that have a single intramolecular disulfide bond. We study the aggregation propensity of the oxidized monomeric fraction and relate these data to an atomistic model of the K18 unfolded ensemble. Our results argue that the main effect of intramolecular disulfide bond formation is to preferentially stabilize conformers within the unfolded ensemble that place the aggregation-prone tau subsequences, PHF6* and PHF6, in conformations that are inconsistent with the formation of cross-β-structure. These data further our understanding of the precise structural features that retard tau aggregation.     

Using a library of structural templates to recognise catalytic sites and explore their evolution in homologous families

Catalytic site structure is normally highly conserved between distantly related enzymes. As a consequence, templates representing catalytic sites have the potential to succeed at function prediction in cases where methods based on sequence or overall structure fail. There are many methods for searching protein structures for matches to structural templates, but few validated template libraries to use with these methods. We present a library of structural templates representing catalytic sites, based on information from the scientific literature. Furthermore, we analyse homologous template families to discover the diversity within families and the utility of templates for active site recognition. Templates representing the catalytic sites of homologous proteins mostly differ by less than 1A root mean square deviation, even when the sequence similarity between the two proteins is low. Within these sets of homologues there is usually no discernible relationship between catalytic site structure similarity and sequence similarity. Because of this structural conservation of catalytic sites, the templates can discriminate between matches to related proteins and random matches with over 85% sensitivity and predictive accuracy. Templates based on protein backbone positions are more discriminating than those based on side-chain atoms. These analyses show encouraging prospects for prediction of functional sites in structural genomics structures of unknown function, and will be of use in analyses of convergent evolution and exploring relationships between active site geometry and chemistry. The template library can be queried via a web server at and is available for download.     

The disulphide bonds in the catalytic domain of BACE are critical but not essential for amyloid precursor protein processing activity

beta-Site APP-cleaving enzyme (BACE) initiates the processing of the amyloid precursor protein (APP) leading to the generation of beta-amyloid, the main component of Alzheimer's disease senile plaques. BACE (Asp2, memapsin 2) is a type I transmembrane aspartic protease responsible for the beta-secretase cleavage of APP producing a soluble form of the ectodomain (sAPPbeta) and the membrane-bound, carboxy-terminal intermediates C99 and C89. BACE maturation involves cysteine bridge formation, N -glycosylation and propeptide removal. We investigated variants of BACE in which the disulphide bonds of the catalytic domain spanning between Cys216/Cys420, Cys278/Cys443 and Cys330/Cys380 were removed by mutagenesis. When transfected in cultured cells, these mutants showed impaired maturation. Nevertheless, a fraction of mutated protein retained both the competence to mature as well as the activity to process APP. For the generation of a functional enzyme the conserved Cys330/Cys380 bond was the most critical, whereas the two bonds between Cys216/Cys420 and Cys278/Cys443, which are typical for the membrane-bound BACE, appeared to be less important.     

Coordinated amino acid changes in homologous protein families

In the tobamovirus coat protein family, amino acid residues at some spatially close positions are found to be substituted in a coordinated manner [Altschuh et al. (1987) J. Mol. Biol., 193, 693]. Therefore, these positions show an identical pattern of amino acid substitutions when amino acid sequences of these homologous proteins are aligned. Based on this principle, coordinated substitutions have been searched for in three additional protein families: serine proteases, cysteine proteases and the haemoglobins. Coordinated changes have been found in all three protein families mostly within structurally constrained regions. This method works with a varying degree of success depending on the function of the proteins, the range of sequence similarities and the number of sequences considered. By relaxing the criteria for residue selection, the method was adapted to cover a broader range of protein families and to study regions of the proteins having weaker structural constraints. The information derived by these methods provides a general guide for engineering of a large variety of proteins to analyse structure-function relationships.     

Formation and transfer of disulphide bonds in living cells

Protein disulphide bonds are formed in the endoplasmic reticulum of eukaryotic cells and the periplasmic space of prokaryotic cells. The main pathways that catalyse the formation of protein disulphide bonds in prokaryotes and eukaryotes are remarkably similar, and they share several mechanistic features. The recent identification of new redox-active proteins in humans and yeast that mechanistically parallel the more established redox-active enzymes indicates that there might be further uncharacterized redox pathways throughout the cell.     

Exhaustive enumeration of protein domain families

Domains are considered as the basic units of protein folding, evolution, and function. Decomposing each protein into modular domains is thus a basic prerequisite for accurate functional classification of biological molecules. Here, we present ADDA, an automatic algorithm for domain decomposition and clustering of all protein domain families. We use alignments derived from an all-on-all sequence comparison to define domains within protein sequences based on a global maximum likelihood model. In all, 90% of domain boundaries are predicted within 10% of domain size when compared with the manual domain definitions given in the SCOP database. A representative database of 249,264 protein sequences were decomposed into 450,462 domains. These domains were clustered on the basis of sequence similarities into 33,879 domain families containing at least two members with less than 40% sequence identity. Validation against family definitions in the manually curated databases SCOP and PFAM indicates almost perfect unification of various large domain families while contamination by unrelated sequences remains at a low level. The global survey of protein-domain space by ADDA confirms that most large and universal domain families are already described in PFAM and/or SMART. However, a survey of the complete set of mobile modules leads to the identification of 1479 new interesting domain families which shuffle around in multi-domain proteins. The data are publicly available at ftp://ftp.ebi.ac.uk/pub/contrib/heger/adda.     

A comparison of sequence and structure protein domain families as a basis for structural genomics.

MOTIVATION: Protein families can be defined based on structure or sequence similarity. We wanted to compare two protein family databases, one based on structural and one on sequence similarity, to investigate to what extent they overlap, the similarity in definition of corresponding families, and to create a list of large protein families with unknown structure as a resource for structural genomics. We also wanted to increase the sensitivity of fold assignment by exploiting protein family HMMs. RESULTS: We compared Pfam, a protein family database based on sequence similarity, to Scop, which is based on structural similarity. We found that 70% of the Scop families exist in Pfam while 57% of the Pfam families exist in Scop. Most families that occur in both databases correspond well to each other, but in some cases they are different. Such cases highlight situations in which structure and sequence approaches differ significantly. The comparison enabled us to compile a list of the largest families that do not occur in Scop; these are suitable targets for structure prediction and determination, and may be useful to guide projects in structural genomics. It can be noted that 13 out of the 20 largest protein families without a known structure are likely transmembrane proteins. We also exploited Pfam to increase the sensitivity of detecting homologs of proteins with known structure, by comparing query sequences to Pfam HMMs that correspond to Scop families. For SWISSPROT+TREMBL, this yielded an increase in fold assignment from 31% to 42% compared to using FASTA only. This method assigned a structure to 22% of the proteins in Saccharomyces cerevisiae, 24% in Escherichia coli, and 16% in Methanococcus jannaschii.     

Formation, isomerisation and reduction of disulphide bonds during protein quality control in the endoplasmic reticulum

Efficient protein folding and quality control in the endoplasmic reticulum (ER) require that disulphide bonds are formed in nascent proteins, isomerised during assisted folding and reduced in terminally misfolded molecules. Recent findings in yeast and mammalian cells indicate that specific protein-protein interactions underlie redox control in the ER, allowing these competing reactions to occur simultaneously during protein quality control.     

C-terminal, endoplasmic reticulum-lumenal domain of prosurfactant protein C - structural features and membrane interactions

Surfactant protein C (SP-C) constitutes the transmembrane part of prosurfactant protein C (proSP-C) and is alpha-helical in its native state. The C-terminal part of proSP-C (CTC) is localized in the endoplasmic reticulum lumen and binds to misfolded (beta-strand) SP-C, thereby preventing its aggregation and amyloid fibril formation. In this study, we investigated the structure of recombinant human CTC and the effects of CTC-membrane interaction on protein structure. CTC forms noncovalent trimers and supratrimeric oligomers. It contains two intrachain disulfide bridges, and its secondary structure is significantly affected by urea or heat only after disulfide reduction. The postulated Brichos domain of CTC, with homologs found in proteins associated with amyloid and proliferative disease, is up to 1000-fold more protected from limited proteolysis than the rest of CTC. The protein exposes hydrophobic surfaces, as determined by CTC binding to the environment-sensitive fluorescent probe 1,1'-bis(4-anilino-5,5'-naphthalenesulfonate). Fluorescence energy transfer experiments further reveal close proximity between bound 1,1'-bis(4-anilino-5,5'-naphthalenesulfonate) and tyrosine residues in CTC, some of which are conserved in all Brichos domains. CTC binds to unilamellar phospholipid vesicles with low micromolar dissociation constants, and differential scanning calorimetry and CD analyses indicate that membrane-bound CTC is less structurally ordered than the unbound protein. The exposed hydrophobic surfaces and the structural disordering that result from interactions with phospholipid membranes suggest a mechanism whereby CTC binds to misfolded SP-C in the endoplasmic reticulum membrane.     

Re-formation of disulphide bonds in reduced antithrombin III.

The separation of pepsin isoenzymes 1, 2, 3 and 5 (gastricsin) in human gastric juice was effected by chromatography on Mono Q ion-exchanger, and slow-moving proteinase was purified to homogeneity by using a modified procedure incorporating a novel affinity-chromatography step. The pH-activity profiles of these enzymes with mucus glycoprotein and basement-membrane substrates were determined; the profiles for pepsin 2 were noticeably different, and, in general, the pH optima for the hydrolysis of basement membrane were more acidic. Pepsin 1 expressed larger specificity constants (kcat./Km) than pepsin 3 with a series of synthetic peptide substrates, reflecting greater binding (smaller Km) by pepsin 1. Inhibitor studies at pH 1.7 and 4.5 with a series of P2-substituted lactoyl-pepstatins implied that valine at position P2 was optimal for inhibiting pepsins 1, 2 and 3 but detrimental for pepsin 5, whereas lysine at position P2 was tolerated well by pepsin 5 but not by pepsins 1, 2 and 3. The potency of lactoyl-pepstatin with lysine at position P2 did not increase as a function of pH. P2-substituted lactoyl-pepstatins failed to show any inhibitory selectivity among pepsins 1, 2 and 3.     

Identification of the disulphide bonds in human platelet glycocalicin

The glycoprotein Ib/IX complex on platelets is responsible for the first stage of haemostasis as an essential component in the primary adhesion of platelets to damaged vessel walls. Glycocalicin is the extracellular part of platelet glycoprotein Ib alpha and contains the von Willebrand factor and thrombin binding sites. Disulphide bonds are implicated in the von Willebrand binding site and studies with peptides point towards a region of glycocalicin with four cysteines as containing the binding sites for both von Willebrand factor and thrombin. The position and linkage of these two disulphide bonds are now determined to be 209-248 and 211-264 and the relevance of this double-loop structure for glycoprotein Ib/IX function is discussed.     

The ProDom database of protein domain families.

The ProDom database contains protein domain families generated from the SWISS-PROT database by automated sequence comparisons. It can be searched on the World Wide Web (http://protein.toulouse.inra. fr/prodom.html ) or by E-mail (prodom@toulouse.inra.fr) to study domain arrangements within known families or new proteins. Strong emphasis has been put on the graphical user interface which allows for interactive analysis of protein homology relationships. Recent improvements to the server include: ProDom search by keyword; links to PROSITE and PDB entries; more sensitive ProDom similarity search with BLAST or WU-BLAST; alignments of query sequences with homologous ProDom domain families; and links to the SWISS-MODEL server (http: //www.expasy.ch/swissmod/SWISS-MODEL.html ) for homology based 3-D domain modelling where possible.     

An integrated approach to the analysis and modeling of protein sequences and structures. III. A comparative study of sequence conservation in protein structural families using multiple structural alignments

The information required to generate a protein structure is contained in its amino acid sequence, but how three-dimensional information is mapped onto a linear sequence is still incompletely understood. Multiple structure alignments of similar protein structures have been used to investigate conserved sequence features but contradictory results have been obtained, due, in large part, to the absence of subjective criteria to be used in the construction of sequence profiles and in the quantitative comparison of alignment results. Here, we report a new procedure for multiple structure alignment and use it to construct structure-based sequence profiles for similar proteins. The definition of "similar" is based on the structural alignment procedure and on the protein structural distance (PSD) described in paper I of this series, which offers an objective measure for protein structure relationships. Our approach is tested in two well-studied groups of proteins; serine proteases and Ig-like proteins. It is demonstrated that the quality of a sequence profile generated by a multiple structure alignment is quite sensitive to the PSD used as a threshold for the inclusion of proteins in the alignment. Specifically, if the proteins included in the aligned set are too distant in structure from one another, there will be a dilution of information and patterns that are relevant to a subset of the proteins are likely to be lost.In order to understand better how the same three-dimensional information can be encoded in seemingly unrelated sequences, structure-based sequence profiles are constructed for subsets of proteins belonging to nine superfolds. We identify patterns of relatively conserved residues in each subset of proteins. It is demonstrated that the most conserved residues are generally located in the regions where tertiary interactions occur and that are relatively conserved in structure. Nevertheless, the conservation patterns are relatively weak in all cases studied, indicating that structure-determining factors that do not require a particular sequential arrangement of amino acids, such as secondary structure propensities and hydrophobic interactions, are important in encoding protein fold information. In general, we find that similar structures can fold without having a set of highly conserved residue clusters or a well-conserved sequence profile; indeed, in some cases there is no apparent conservation pattern common to structures with the same fold. Thus, when a group of proteins exhibits a common and well-defined sequence pattern, it is more likely that these sequences have a close evolutionary relationship rather than the similarities having arisen from the structural requirements of a given fold.     

Reduction of disulphide bonds in human haptoglobin 2-1

Gradual reduction of disulphide bonds in human haptoglobin, type 2-1, was carried out either by the use of sodium borohydride or 2-mercaptoethanol, newly exposed sulphydryl groups as determined by the Ellman's reagent and by the incorporation of [(14) C] acetamide, respectively. Cleavage of disulphide bonds resulted in the formation of a number of intermediates, separated in polyacrylamide gel electrophoresis, with sulphydryl groups blocked by the radioactive label. On the basis of molecular mass determinations, subunit composition of intermediates, was postulated. The ability of haptoglobin to form an active peroxidase-like complex with hemoglobin depended to a considerable extent on the presence of intact disulphide bonds. On the contrary,throughout the course of reduction of inter- and intrachain disulphides, antigenic reactivity was found to remain unchanged.     

Deterministic features of side-chain main-chain hydrogen bonds in globular protein structures

A total of 19 835 polar residues from a data set of 250 non-homologous and highly resolved protein crystal structures were used to identify side-chain main-chain (SC-MC) hydrogen bonds. The ratio of the number of SC-MC hydrogen bonds to the total number of polar residues is close to 1:2, indicating the ubiquitous nature of such hydrogen bonds. Close to 56% of the SC-MC hydrogen bonds are local involving side-chain acceptor/donor ('i') and a main-chain donor/acceptor within the window i-5 to i+5. These short-range hydrogen bonds form well defined conformational motifs characterized by specific combinations of backbone and side-chain torsion angles. (a) The Ser/Thr residues show the greatest preference in forming intra-helical hydrogen bonds between the atoms O(gamma)(i) and O(i-4). More than half the examples of such hydrogen bonds are found at the middle of alpha-helices rather than at their ends. The most favoured motif of these examples is alpha(R)alpha(R)alpha(R)alpha(R)(g(-)). (b) These residues also show great preference to form hydrogen bonds between O(gamma)(i) and O(i-3), which are closely related to the previous type and though intra-helical, these hydrogen bonds are more often found at the C-termini of helices than at the middle. The motif represented by alpha(R)alpha(R)alpha(R)alpha(R)(g(+)) is most preferred in these cases. (c) The Ser, Thr and Glu are the most frequently found residues participating in intra-residue hydrogen bonds (between the side-chain and main-chain of the same residue) which are characterized by specific motifs of the form beta(g(+)) for Ser/Thr residues and alpha(R)(g(-)g(+)t) for Glu/Gln. (d) The side-chain acceptor atoms of Asn/Asp and Ser/Thr residues show high preference to form hydrogen bonds with acceptors two residues ahead in the chain, which are characterized by the motifs beta (tt')alphaR and beta(t)alpha(R), respectively. These hydrogen bonded segments, referred to as Asx turns, are known to provide stability to type I and type I' beta-turns. (e) Ser/Thr residues often form a combination of SC-MC hydrogen bonds, with the side-chain donor hydrogen bonded to the carbonyl oxygen of its own peptide backbone and the side-chain acceptor hydrogen bonded to an amide hydrogen three residues ahead in the sequence. Such motifs are quite often seen at the beginning of alpha-helices, which are characterized by the beta(g(+))alpha(R)alpha(R) motif. A remarkable majority of all these hydrogen bonds are buried from the protein surface, away from the surrounding solvent. This strongly indicates the possibility of side-chains playing the role of the backbone, in the protein interiors, to satisfy the potential hydrogen bonding sites and maintaining the network of hydrogen bonds which is crucial to the structure of the protein.