Conservation analysis and structure prediction of the SH2 family of phosphotyrosine binding domains.

Src homology 2 (SH2) regions are short (approximately 100 amino acids), non-catalytic domains conserved among a wide variety of proteins involved in cytoplasmic signaling induced by growth factors. It is thought that SH2 domains play an important role in the intracellular response to growth factor stimulation by binding to phosphotyrosine containing proteins. In this paper we apply the techniques of multiple sequence alignment, secondary structure prediction and conservation analysis to 67 SH2 domain amino acid sequences. This combined approach predicts seven core secondary structure regions with the pattern beta-alpha-beta-beta-beta-beta-alpha, identifies those residues most likely to be buried in the hydrophobic core of the native SH2 domain, and highlights patterns of conservation indicative of secondary structural elements. Residues likely to be involved in phosphotyrosine binding are shown and orientations of the predicted secondary structures suggested which could enable such residues to cooperate in phosphate binding. We propose a consensus pattern that encapsulates the principal conserved features of the SH2 domains. Comparison of the proposed SH2 domain of akt to this pattern shows only 12/40 matches, suggesting that this domain may not exhibit SH2-like properties.     

Ligand-induced strain in hydrogen bonds of the c-Src SH3 domain detected by NMR

Changes in the molecular conformation of proteins can result from a variety of perturbations, and can play crucial roles in the regulation of biological activity. A new solution NMR method has been applied to monitor ligand-induced changes in hydrogen bond geometry in the chicken c-Src SH3 domain. The structural response of this domain to ligand binding has been investigated by measuring trans-hydrogen bond (15)N-(13)C' scalar couplings in the free state and when bound to the high affinity class I ligand RLP2, containing residues RALPPLPRY. A comparison between hydrogen bonds in high resolution X-ray structures of this domain and those observed via (h3)J(NC') couplings in solution shows remarkable agreement. Two backbone-to-side-chain hydrogen bonds are observed in solution, and each appears to play a role in stabilization of loop structure. Reproducible ligand-induced changes in trans-hydrogen bond scalar couplings are observed across the domain that translate into changes in hydrogen bond length ranging between 0.02 to 0.12 A. The observed changes can be rationalized by an induced fit mechanism in which hydrogen bonds across the protein participate in a compensatory response to forces imparted at the protein-ligand interface. Upon ligand binding, mutual intercalation of the two Leu-Pro segments of the ligand between three aromatic side-chains protruding from the SH3 surface wedges apart secondary structural elements within the SH3 domain. This disruption is transmitted in a domino-like effect across the domain through networks of hydrogen bonded peptide planes. The unprecedented resolution obtained demonstrates the ability to characterize subtle structural rearrangements within a protein upon perturbation, and represents a new step in the endeavor to understand how hydrogen bonds contribute to the stabilization and function of biological macromolecules.     

Sequence and structural analysis of cellular retinoic acid-binding proteins reveals a network of conserved hydrophobic interactions

Proteins in the intracellular lipid-binding protein (iLBP) family show remarkably high structural conservation despite their low-sequence identity. A multiple-sequence alignment using 52 sequences of iLBP family members revealed 15 fully conserved positions, with a disproportionately high number of these (n=7) located in the relatively small helical region. The conserved positions displayed high structural conservation based on comparisons of known iLBP crystal structures. It is striking that the beta-sheet domain had few conserved positions, despite its high structural conservation. This observation prompted us to analyze pair-wise interactions within the beta-sheet region to ask whether structural information was encoded in interacting amino acid pairs. We conducted this analysis on the iLBP family member, cellular retinoic acid-binding protein I (CRABP I), whose folding mechanism is under study in our laboratory. Indeed, an analysis based on a simple classification of hydrophobic and polar amino acids revealed a network of conserved interactions in CRABP I that cluster spatially, suggesting a possible nucleation site for folding. Significantly, a small number of residues participated in multiple conserved interactions, suggesting a key role for these sites in the structure and folding of CRABP I. The results presented here correlate well with available experimental evidence on folding of CRABPs and their family members and suggest future experiments. The analysis also shows the usefulness of considering pair-wise conservation based on a simple classification of amino acids, in analyzing sequences and structures to find common core regions among homologues.     

Conserved water molecules in MHC class-I molecules and their putative structural and functional roles

A set of conserved water positions making direct contacts with the alpha1 and alpha2 domains of the MHC class-I protein was identified by a cluster analysis in 12 high-resolution crystal structures of proteins from different allele types and different species, comprising human, mouse and rat. The analysis revealed a total of 63 clusters, corresponding to water molecules, whose positions are conserved in half or more of the analyzed structures. Analysis of these clusters shows that the most conserved water positions-those appearing in the largest fraction of the structures-were also the most accurately defined, as measured by their normalized crystallographic B-factor. Not too surprisingly, these positions displayed better overlap and formed more H-bonds with the protein. In a second part of this work, a detailed analysis is presented of three of the most conserved water positions and their putative structural and functional roles are discussed. The most highly conserved of the three appears to play an important role in stabilizing the conformation of a twisted beta-turn between residues 118 and 122 (numbering of HLA-B3501, PDB code 1A1N). An equivalent water molecule was found to be associated with a similar beta-turn in 43 unrelated structures surveyed in the PDB, leading to the suggestion that this water molecule plays an important structural role in this type of turn. The second water molecule makes hydrogen bonds with residues lining pocket B in the peptide-binding groove and is suggested to play a role in modulating peptide recognition. The third highly conserved water molecule is located at the first kink of the alpha2 helix, possibly playing a role in determining the position of the N-terminal segment of that helix, which also carries side chains in contact with the bound peptide. This information on conserved water positions in MHC class-I molecules should be helpful in modeling interactions with bound peptide antigens and in designing new peptides with tailor-made affinities.     

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.     

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long‐range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.     

The guanidino-group modifying enzymes: structural basis for their diversity and commonality

The guanidino-group modifying enzyme (GME) superfamily contains many drug targets, including metabolic enzymes from pathogenic microorganisms as well as key regulatory proteins from higher eukaryotes. These enzymes, despite their diverse sequences, adopt the common alpha/beta propeller fold and catalyze the modification of (methylated) guanidino groups. Our structural superposition and structure-based alignment for the GMEs have identified key residues that are involved in the catalysis and substrate binding. We have shown that conserved guanidino-carboxyl interactions are utilized in two different ways; the acidic residues in the catalytic site form hydrogen bonds to the substrate guanidino group, and the enzyme Arg residues at several key positions recognize the carboxyl group of the substrate and fix its orientation. Based on this observation, we have proposed rules for classifying the GME sequences and predicting their molecular function from the conservation of the key acidic and Arg residues. Other novel motifs have been identified, which involve residues that are not in direct contact with the substrate but are likely to stabilize the active-site conformation through hydrogen-bonding networks. In addition, we have examined the domain architecture of the GMEs. Although most members consist of a single catalytic domain, fold recognition analysis has identified a likely bifunctional enzyme from a cyanobacterium. It has also revealed common immunoglobulin-like beta-sandwich domains found in the enzymes that recognize protein substrates. These findings will be useful for predicting the precise mechanism of action for potential novel targets and designing therapeutic compounds against them.     

The CD8alpha from sea bass (Dicentrarchus labrax L.): Cloning, expression and 3D modelling

In this paper we describe the cloning, expression and structural study by modelling techniques of the CD8alpha from sea bass (Dicentrarchus labrax L.). The sea bass CD8alpha cDNA is comprised of 1490 bp and is translated in one reading frame to give a protein of 217 amino acids, with a predicted 26 amino acids signal peptide, a 88 bp 5'-UTR and a 748 bp 3'-UTR. A multiple alignment of CD8alpha from sea bass with other known CD8alpha sequences shows the conservation of most amino acid residues involved in the peculiar structural domains found within CD8alpha's. Cysteine residues that are involved in disulfide bonding to form the V domain are conserved. In contrast, an extra cysteine residue found in most mammals in this region is not present in sea bass. The transmembrane and cytoplasmic regions are the most conserved regions within the molecule in the alignment analysis. However, the motif (CXCP) that is thought to be responsible for binding p56lck is missing in the sea bass sequence. Phylogenetic analysis conducted using amino acid sequences showed that sea bass CD8alpha grouped with other known teleost sequences and that three different clusters were formed by the mammalian, avian and fish CD8alpha sequences. The thymus was the tissue with the highest CD8alpha expression, followed by gut, gills, peripheral blood leukocytes and spleen. Lower CD8alpha mRNA levels were found in head kidney, liver and brain. It was possible to create a partial 3D model using the human and mouse structures as template. The CD8alpha 11-120 amino acid region was taken into consideration and the best obtained 3D model shows the presence of ten beta-strands, involving about 50% of the sequence. The global structure was defined as an immunoglobulin-like beta-sandwich made of two anti-parallel sheets. Two cysteines were present in this region and they were at a suitable distance to form an S-S bond as seen in the template human and mouse structures.     

Relationship between sequence conservation and three-dimensional structure in a large family of esterases, lipases, and related proteins.

Based on the recently determined X-ray structures of Torpedo californica acetylcholinesterase and Geotrichum candidum lipase and on their three-dimensional superposition, an improved alignment of a collection of 32 related amino acid sequences of other esterases, lipases, and related proteins was obtained. On the basis of this alignment, 24 residues are found to be invariant in 29 sequences of hydrolytic enzymes, and an additional 49 are well conserved. The conservation in the three remaining sequences is somewhat lower. The conserved residues include the active site, disulfide bridges, salt bridges, and residues in the core of the proteins. Most invariant residues are located at the edges of secondary structural elements. A clear structural basis for the preservation of many of these residues can be determined from comparison of the two X-ray structures.     

Conservation of structure and function among tyrosine recombinases: homology-based modeling of the lambda integrase core-binding domain

Tyrosine recombinases participate in diverse biological processes by catalyzing recombination between specific DNA sites. Although a conserved protein fold has been described for the catalytic (CAT) domains of five recombinases, structural relationships between their core-binding (CB) domains remain unclear. Despite differences in the specificity and affinity of core-type DNA recognition, a conserved binding mechanism is suggested by the shared two-domain motif in crystal structure models of the recombinases Cre, XerD and Flp. We have found additional evidence for conservation of the CB domain fold. Comparison of XerD and Cre crystal structures showed that their CB domains are closely related; the three central α-helices of these domains are superposable to within 1.44 Å. A structure-based multiple sequence alignment containing 25 diverse CB domain sequences provided evidence for widespread conservation of both structural and functional elements in this fold. Based upon the Cre and XerD crystal structures, we employed homology modeling to construct a three-dimensional structure for the λ integrase CB domain. The model provides a conceptual framework within which many previously identified, functionally important amino acid residues were investigated. In addition, the model predicts new residues that may participate in core-type DNA binding or dimerization, thereby providing hypotheses for future genetic and biochemical experiments.     

Determinants of the SRC homology domain 3-like fold

In eukaryotes, the Src homology domain 3 (SH3) is a very important motif in signal transduction. SH3 domains recognize poly-proline-rich peptides and are involved in protein-protein interactions. Until now, the existence of SH3 domains has not been demonstrated in prokaryotes. However, the structure of the C-terminal domain of DtxR clearly shows that the fold of this domain is very similar to that of the SH3 domain. In addition, there is evidence that the C-terminal domain of DtxR binds to poly-proline-rich regions. Other bacterial proteins have domains that are structurally similar to the SH3 domain but whose functions are unknown or differ from that of the SH3 domain. The observed similarities between the structures of the C-terminal domain of DtxR and the SH3 domain constitute a perfect system to gain insight into their function and information about their evolution. Our results show that the C-terminal domain of DtxR shares a number of conserved key hydrophobic positions not recognizable from sequence comparison that might be responsible for the integrity of the SH3-like fold. Structural alignment of an ensemble of such domains from unrelated proteins shows a common structural core that seems to be conserved despite the lack of sequence similarity. This core constitutes the minimal requirements of protein architecture for the SH3-like fold.     

Allosteric loss-of-function mutations in HIV-1 Nef from a long-term non-progressor

Activation of Src family kinases by human immunodeficiency virus type 1 (HIV-1) Nef may play an important role in the pathogenesis of HIV/AIDS. Here we investigated whether diverse Nef sequences universally activate Hck, a Src family member expressed in macrophages and other HIV-1 target cells. In general, we observed that Hck activation is a highly conserved Nef function. However, we identified an unusual Nef variant from an HIV-positive individual that did not develop AIDS which failed to activate Hck despite the presence of conserved residues linked to Hck SH3 domain binding and kinase activation. Amino acid sequence alignment with active Nef proteins revealed differences in regions not previously implicated in Hck activation, including a large internal flexible loop absent from available Nef structures. Substitution of these residues in active Nef compromised Hck activation without affecting SH3 domain binding. These findings show that residues at a distance from the SH3 domain binding site influence Nef interactions allosterically with a key effector protein linked to AIDS progression.     

Altered dynamics in Lck SH3 upon binding to the LBD1 domain of Herpesvirus saimiri Tip

The Tip protein from Herpesvirus saimiri interacts with the SH3 domain from the Src-family kinase Lck via a proline-containing sequence termed LBD1. Src-family kinase SH3 domains related to Lck have been shown to be dynamic in solution and partially unfold under physiological conditions. The rate of such partial unfolding is reduced by viral protein binding. To determine if the Lck SH3 domain displayed similar behavior, the domain was investigated with hydrogen exchange and mass spectrometry. Lck SH3 was found to be highly dynamic in solution. While other SH3 domains require as much as 10,000 sec to become totally deuterated, Lck SH3 became almost completely labeled within 200 sec. A partial unfolding event involving 8-10 residues was observed with a half-life of approximately 10 sec. Tip LBD1 binding did not cause gross structural changes in Lck SH3 but globally stabilized the domain and reduced the rate of partial unfolding by a factor of five. The region of partial unfolding in Lck SH3 was found to be similar to that identified for other SH3 domains that partially unfold. Although the sequence conservation between Lck SH3 and other closely related SH3 domains is high, the dynamics do not appear to be conserved.     

Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: The case of the fold-type I, pyridoxal-5'-phosphate-dependent enzymes

The wealth of biological information provided by structural and genomic projects opens new prospects of understanding life and evolution at the molecular level. In this work, it is shown how computational approaches can be exploited to pinpoint protein structural features that remain invariant upon long evolutionary periods in the fold-type I, PLP-dependent enzymes. A nonredundant set of 23 superposed crystallographic structures belonging to this superfamily was built. Members of this family typically display high-structural conservation despite low-sequence identity. For each structure, a multiple-sequence alignment of orthologous sequences was obtained, and the 23 alignments were merged using the structural information to obtain a comprehensive multiple alignment of 921 sequences of fold-type I enzymes. The structurally conserved regions (SCRs), the evolutionarily conserved residues, and the conserved hydrophobic contacts (CHCs) were extracted from this data set, using both sequence and structural information. The results of this study identified a structural pattern of hydrophobic contacts shared by all of the superfamily members of fold-type I enzymes and involved in native interactions. This profile highlights the presence of a nucleus for this fold, in which residues participating in the most conserved native interactions exhibit preferential evolutionary conservation, that correlates significantly (r = 0.70) with the extent of mean hydrophobic contact value of their apolar fraction.     

Protein folds propelled by diversity

Many proteins involved in key biological processes are modular in nature. A group of these, the beta-propeller proteins, fold by packing 4-stranded beta-sheets in a circular array. The members of this group are increasingly numerous and, although their modular building blocks all preserve the same basic conformation, they do not have similar sequences. These proteins have extreme functional and phylogenetic diversity. Here, features of the beta-propeller fold are reviewed through comparisons of available structural coordinates. Structure-based sequence alignments combined with analyses of superpositions of individual modular units reveal conserved general features such as hydrogen bonds, beta-turns and positions of hydrophobic contacts. The lack of significant sequence identity is compensated by sets of interactions which stabilise the fold differently in distinct structures. Re-occurring aspartates make contacts to exposed backbone amides in turns or peptide connections within the same sheet. The sole factor responsible for the number of sheets that assemble in the array is the size of the hydrophobic residues that pack into the cores between the sheets. Whilst there is no overall sequence conservation, it may be possible to detect new members of this fold through sequence searches that take into account the repeated nature of the modular assembly as well as the positions of hydrophobic residues and H-bonding side chains.     

Folding of a highly conserved diverging turn motif from the SH3 domain

Recent NMR structural characterization studies showed that a seven-residue segment (FKKGERL) from the src SH3 domain adopts the nativelike diverging type II beta-turn in aqueous solution in support of the prediction based on the I-sites library of sequence structural motifs. We study the conformational variability and folding/unfolding thermodynamics of this peptide in explicit solvent using replica-exchange molecular dynamics simulations, which greatly enhances the sampling of the conformational space. This peptide samples three main free energy basins (nativelike, intermediate, and unfolded) separated by small barriers. The nativelike basin is fractionally populated (DeltaG(300K) = 0.4 kcal/mol) with structures that satisfy a subset of the NMR-derived constraints. The intrinsic stability of the diverging turn is examined in relationship to the nature of three specific contacts: a turn-hydrogen bond, a mainchain-to-sidechain hydrogen bond, and an end-to-end hydrophobic contact. We have carried out simulations of mutants at the highly conserved GE positions in the sequence. The mutation E5D destabilizes the isolated diverging turn motif, contrary to the observation that this mutation stabilizes the fyn SH3 domain. The G4T mutation also destabilizes the isolated diverging turn; however, the extent of destabilization is smaller than that of the reverse mutation in the drk SH3.     

Molecular Determinants of PAM2 Recognition by the MLLE Domain of Poly(A)-Binding Protein

MLLE (previously known as PABC) is a peptide-binding domain that is found in poly(A)-binding protein (PABP) and EDD (E3 isolated by differential display), a HECT E3 ubiquitin ligase also known as HYD (hyperplastic discs tumor suppressor) or UBR5. The MLLE domain from PABP recruits various regulatory proteins and translation factors to poly(A) mRNAs through binding of a conserved 12 amino acid peptide motif called PAM2 (for PABP-interacting motif 2). Here, we determined crystal structures of the MLLE domain from PABP alone and in complex with PAM2 peptides from PABP-interacting protein 2. The structures provide a detailed view of hydrophobic determinants of the MLLE binding coded by PAM2 positions 3, 5, 7, 10, and 12 and reveal novel intermolecular polar contacts. In particular, the side chain of the invariant MLLE residue K580 forms hydrogen bonds with the backbone of PAM2 residues 5 and 7. The structures also show that peptide residues outside of the conserved PAM2 motif contribute to binding. Altogether, the structures provide a significant advance in understanding the molecular basis for the binding of PABP by PAM2-containing proteins involved in translational control, mRNA deadenylation, and other cellular processes.