Chirality-independent protein-protein recognition between transmembrane domains in vivo

Stereospecificity in protein-protein recognition and docking is an unchallenged dogma. Soluble proteins provide the main source of evidence for stereospecificity. In contrast, within the membrane little is known about the role of stereospecificity in the recognition process. Here, we have reassessed the stereospecificity of protein-protein recognition by testing whether it holds true for the well-defined glycophorin A (GPA) transmembrane domain in vivo. We found that the all-D amino acid GPA transmembrane domain and two all-D mutants specifically associated with an all-L GPA transmembrane domain, within the membrane milieu of Escherichia coli. Molecular dynamics techniques reveal a possible structural explanation to the observed interaction between all-D and all-L transmembrane domains. A very strong correlation was found between amino acid residues at the interface of both the all-L homodimer structure and the mixed L/D heterodimer structure, suggesting that the original interactions are conserved. The results suggest that GPA helix-helix recognition within the membrane is chirality-independent.     

Both transmembrane domains of SecG contribute to signal sequence recognition by the Escherichia coli protein export machinery

A chimeric protein containing the uncleaved signal sequence of plasminogen activators inhibitor-2 (PAI2) fused to alkaline phosphatase (AP) interferes with Escherichia coli protein export and arrests growth. Suppressors of this toxicity include secG mutations that define the Thr-41-Leu-42-Phe-43 (TLF) domain of SecG. These mutations slow down the export of PAI2-AP. Another construct encoding a truncated PAI2 signal sequence (hB-AP) is also toxic. Most suppressors exert their effect on both chimeric proteins. We describe here five secG suppressors that only suppress the toxicity of hB-AP and selectively slow down its export. These mutations do not alter the TLF domain: three encode truncated SecG, whereas two introduce Arg residues in the transmembrane domains of SecG. The shortest truncated protein only contains 13 residues of SecG, suggesting that the mutation is equivalent to a null allele. Indeed, a secG disruption selectively suppresses the toxicity of hB-AP. However, the missense mutations are not null alleles. They allow SecG binding to SecYE, although with reduced affinity. Furthermore, these mutated SecG are functional, as they facilitate the export of endogenous proteins. Thus, SecG participates in signal sequence recognition, and both transmembrane domains of SecG contribute to ensure normal signal sequence recognition by the translocase.     

Phosphorylation by cAMP-dependent protein kinase modulates the structural coupling between the transmembrane and cytosolic domains of phospholamban

We have used frequency-domain fluorescence spectroscopy to investigate the structural linkage between the transmembrane and cytosolic domains of the regulatory protein phospholamban (PLB). Using an engineered PLB having a single cysteine (Cys(24)) derivatized with the fluorophore N-(1-pyrenyl)maleimide (PMal), we have used fluorescence resonance energy transfer (FRET) to measure the average spatial separation and conformational heterogeneity between PMal bound to Cys(24) in the transmembrane domain and Tyr(6) in the cytosolic domain near the amino terminus of PLB. In these measurements, PMal serves as a FRET donor, and Tyr(6) serves as a FRET acceptor following its nitration by tetranitromethane. The native structure of PLB is retained following site-directed mutagenesis and chemical modification, as indicated by the ability of the derivatized PLB to fully regulate the Ca-ATPase following their co-reconstitution. To assess how phosphorylation modulates the structure of PLB itself, FRET measurements were made following reconstitution of PLB in membrane vesicles made from extracted sarcoplasmic reticulum membrane lipids. We find that the cytosolic domain of PLB assumes a wide range of conformations relative to the transmembrane sequence, consistent with other structural data indicating the presence of a flexible hinge region between the transmembrane and cytosolic domains of PLB. Phosphorylation of Ser(16) by PKA results in a 3 A decrease in the spatial separation between PMal at Cys(24) and nitroTyr(6) and an almost 2-fold decrease in conformational heterogeneity, suggesting a stabilization of the hinge region of PLB possibly through an electrostatic linkage between phosphoSer(16) and Arg(13) that promotes a coil-to-helix transition. This structural transition has the potential to function as a conformational switch, since inhibition of the Ca-ATPase requires disruption of the secondary structure of PLB in the vicinity of the hinge element to permit association with the nucleotide binding domain at a site located approximately 50 A above the membrane surface. Following phosphorylation, the stabilization of the helical content in the hinge domain will disrupt this inhibitory interaction by reducing the maximal dimension of the cytosolic domain of PLB. Thus, stabilization of the structure of PLB following phosphorylation of Ser(16) is part of a switching mechanism, which functions to alter binding interactions between PLB and the nucleotide binding domain of the Ca-ATPase that modulates enzyme inhibition.     

Factors limiting the performance of prediction-based fold recognition methods

In the past few years, a new generation of fold recognition methods has been developed, in which the classical sequence information is combined with information obtained from secondary structure and, sometimes, accessibility predictions. The results are promising, indicating that this approach may compete with potential-based methods (Rost B et al., 1997, J Mol Biol 270:471-480). Here we present a systematic study of the different factors contributing to the performance of these methods, in particular when applied to the problem of fold recognition of remote homologues. Our results indicate that secondary structure and accessibility prediction methods have reached an accuracy level where they are not the major factor limiting the accuracy of fold recognition. The pattern degeneracy problem is confirmed as the major source of error of these methods. On the basis of these results, we study three different options to overcome these limitations: normalization schemes, mapping of the coil state into the different zones of the Ramachandran plot, and post-threading graphical analysis.     

Diversity in protein recognition by PTB domains

Phosphotyrosine-binding (PTB) domains were originally identified as modular domains that recognize phosphorylated Asn-Pro-Xxx-p Tyr-containing proteins. Recent binding and structural studies of PTB domain complexes with target peptides have revealed a number of deviations from the previously described mode of interaction, with respect to both the sequences of possible targets and their structures within the complexes. This diversity of recognition by PTB domains extends and strengthens our general understanding of modular binding domain recognition.     

A single aldehyde group can serve as a structural element for recognition by transmembrane protein CD36

Transmembrane protein CD36 is considered to bind its distinct ligands such as long-chain fatty acids primarily by recognizing their terminal carboxyl moiety. In this study, we provide evidence that long-chain fatty aldehydes, such as oleic aldehyde, can be recognized by CD36. We suggest that a single aldehyde group may also serve as one of the structural elements recognizable by CD36.     

Conserved structural determinants in three-fingered protein domains

The three-dimensional structures of some components of snake venoms forming so-called 'three-fingered protein' domains (TFPDs) are similar to those of the ectodomains of activin, bone morphogenetic protein and transforming growth factor-beta receptors, and to a variety of proteins encoded by the Ly6 and Plaur genes. The analysis of sequences of diverse snake toxins, various ectodomains of the receptors that bind activin and other cytokines, and numerous gene products encoded by the Ly6 and Plaur families of genes has revealed that they differ considerably from each other. The sequences of TFPDs may consist of up to six disulfide bonds, three of which have the same highly conserved topology. These three disulfide bridges and an asparagine residue in the C-terminal part of TFPDs are essential for the TFPD-like fold. Analyses of the three-dimensional structures of diverse TFPDs have revealed that the three highly conserved disulfides impose a major stabilizing contribution to the TFPD-like fold, in both TFPDs contained in some snake venoms and ectodomains of several cellular receptors, whereas the three remaining disulfide bonds impose specific geometrical constraints in the three fingers of some TFPDs.     

The use of functional domains to improve transmembrane protein topology prediction

Transmembrane proteins affect vital cellular functions and pathogenesis, and are a focus of drug design. It is difficult to obtain diffraction quality crystals to study transmembrane protein structure. Computational tools for transmembrane protein topology prediction fill in the gap between the abundance of transmembrane proteins and the scarcity of known membrane protein structures. Their prediction accuracy is still inadequate: TMHMM, the current state-of-the-art method, has less than 52% accuracy in topology prediction on one set of transmembrane proteins of known topology. Based on the observation that there are functional domains that occur preferentially internal or external to the membrane, we have extended the model of TMHMM to incorporate functional domains, using a probabilistic approach originally developed for computational gene finding. Our extension is better than TMHMM in predicting the topology of transmembrane proteins. As prediction of functional domain improves, our system's prediction accuracy will likely improve as well.     

Perfringolysin O association with ordered lipid domains: implications for transmembrane protein raft affinity

Upon interaction with cholesterol, perfringolysin O (PFO) inserts into membranes and forms a rigid transmembrane (TM) β-barrel. PFO is believed to interact with liquid ordered lipid domains (lipid rafts). Because the origin of TM protein affinity for rafts is poorly understood, we investigated PFO raft affinity in vesicles having coexisting ordered and disordered lipid domains. Fluorescence resonance energy transfer (FRET) from PFO Trp to domain-localized acceptors indicated that PFO generally has a raft affinity between that of LW peptide (low raft affinity) and cholera toxin B (high raft affinity) in vesicles containing ordered domains rich in brain sphingomyelin or distearoylphosphatidylcholine. FRET also showed that ceramide, which increases exposure of cholesterol to water and thus displaces it from rafts, does not displace PFO from ordered domains. This can be explained by shielding of PFO-bound cholesterol from water. Finally, FRET showed that PFO affinity for ordered domains was higher in its non-TM (prepore) form than in its TM form, demonstrating that the TM portion of PFO interacts unfavorably with rafts. Microscopy studies in giant unilamellar vesicles confirmed that PFO exhibits intermediate raft affinity, and showed that TM PFO (but not non-TM PFO) concentrated at the edges of liquid ordered domains. These studies suggest that a combination of binding to raft-associating molecules and having a rigid TM structure that is unable to pack well in a highly ordered lipid environment can control TM protein domain localization. To accommodate these constraints, raft-associated TM proteins in cells may tend to locate within liquid disordered shells encapsulated within ordered domains.     

Guiding probabilistic search of the protein conformational space with structural profiles

The roughness of the protein energy surface poses a significant challenge to search algorithms that seek to obtain a structural characterization of the native state. Recent research seeks to bias search toward near-native conformations through one-dimensional structural profiles of the protein native state. Here we investigate the effectiveness of such profiles in a structure prediction setting for proteins of various sizes and folds. We pursue two directions. We first investigate the contribution of structural profiles in comparison to or in conjunction with physics-based energy functions in providing an effective energy bias. We conduct this investigation in the context of Metropolis Monte Carlo with fragment-based assembly. Second, we explore the effectiveness of structural profiles in providing projection coordinates through which to organize the conformational space. We do so in the context of a robotics-inspired search framework proposed in our lab that employs projections of the conformational space to guide search. Our findings indicate that structural profiles are most effective in obtaining physically realistic near-native conformations when employed in conjunction with physics-based energy functions. Our findings also show that these profiles are very effective when employed instead as projection coordinates to guide probabilistic search toward undersampled regions of the conformational space.     

An ankyrin-like protein with transmembrane domains is specifically lost after oncogenic transformation of human fibroblasts

We have identified a novel transformation-sensitive mRNA, which is present in cultured fibroblasts but is lacking in SV40 transformed cells as well as in many mesenchymal tumor cell lines. The corresponding gene is located on human chromosome 8 in band 8q13. The open reading frame of the mRNA encodes a protein of 1119 amino acids forming two distinct domains. The N-terminal domain consists of 18 repeats that are related to the cytoskeletal protein ankyrin. The C-terminal domain contains six putative transmembrane segments that resemble many ion channels. This overall structure is reminiscent of TRP-like proteins that function as store-operated calcium channels. The novel protein with an Mr of 130 kDa is expressed at a very low level in human fibroblasts and at a moderate level in liposarcoma cells. Overexpression in eukaryotic cells appears to interfere with normal growth, suggesting that it might play a direct or indirect role in signal transduction and growth control.     

Towards an automatic classification of protein structural domains based on structural similarity

Background  

Formal classification of a large collection of protein structures aids the understanding of evolutionary relationships among them. Classifications involving manual steps, such as SCOP and CATH, face the challenge of increasing volume of available structures. Automatic methods such as FSSP or Dali Domain Dictionary, yield divergent classifications, for reasons not yet fully investigated. One possible reason is that the pairwise similarity scores used in automatic classification do not adequately reflect the judgments made in manual classification. Another possibility is the difference between manual and automatic classification procedures. We explore the degree to which these two factors might affect the final classification.     

3Dee: a database of protein structural domains

The 3Dee database is a repository of protein structural domains. It stores alternative domain definitions for the same protein, organises domains into sequence and structural hierarchies, contains non-redundant set(s) of sequences and structures, multiple structure alignments for families of domains, and allows previous versions of the database to be regenerated. AVAILABILITY: 3Dee is accessible on the World Wide Web at the URL http://barton.ebi.ac.uk/servers/3Dee.html.     

Enhanced genome annotation using structural profiles in the program 3D-PSSM

A method (three-dimensional position-specific scoring matrix, 3D-PSSM) to recognise remote protein sequence homologues is described. The method combines the power of multiple sequence profiles with knowledge of protein structure to provide enhanced recognition and thus functional assignment of newly sequenced genomes. The method uses structural alignments of homologous proteins of similar three-dimensional structure in the structural classification of proteins (SCOP) database to obtain a structural equivalence of residues. These equivalences are used to extend multiply aligned sequences obtained by standard sequence searches. The resulting large superfamily-based multiple alignment is converted into a PSSM. Combined with secondary structure matching and solvation potentials, 3D-PSSM can recognise structural and functional relationships beyond state-of-the-art sequence methods. In a cross-validated benchmark on 136 homologous relationships unambiguously undetectable by position-specific iterated basic local alignment search tool (PSI-Blast), 3D-PSSM can confidently assign 18 %. The method was applied to the remaining unassigned regions of the Mycoplasma genitalium genome and an additional 13 regions were assigned with 95 % confidence. 3D-PSSM is available to the community as a web server: http://www.bmm.icnet.uk/servers/3dpssm Copyright 2000 Academic Press.     

Different transmembrane domains associate with distinct endoplasmic reticulum components during membrane integration of a polytopic protein

We have been studying the insertion of the seven transmembrane domain (TM) protein opsin to gain insights into how the multiple TMs of polytopic proteins are integrated at the endoplasmic reticulum (ER). We find that the ER components associated with the first and second TMs of the nascent opsin polypeptide chain are clearly distinct. The first TM (TM1) is adjacent to the alpha and beta subunits of the Sec61 complex, and a novel component, a protein associated with the ER translocon of 10 kDa (PAT-10). The most striking characteristic of PAT-10 is that it remains adjacent to TM1 throughout the biogenesis and membrane integration of the full-length opsin polypeptide. TM2 is also found to be adjacent to Sec61alpha and Sec61beta during its membrane integration. However, TM2 does not form any adducts with PAT-10; rather, a transient association with the TRAM protein is observed. We show that the association of PAT-10 with opsin TM1 does not require the N-glycosylation of the nascent chain and occurs irrespective of the amino acid sequence and transmembrane topology of TM1. We conclude that the precise makeup of the ER membrane insertion site can be distinct for the different transmembrane domains of a polytopic protein. We find that the environment of a particular TM can be influenced by both the "stage" of nascent chain biosynthesis reached, and the TM's relative location within the polypeptide.     

Using protein structural information in evolutionary inference: transmembrane proteins

We present a model of amino acid sequence evolution based on a hidden Markov model that extends to transmembrane proteins previous methods that incorporate protein structural information into phylogenetics. Our model aims to give a better understanding of processes of molecular evolution and to extract structural information from multiple alignments of transmembrane sequences and use such information to improve phylogenetic analyses. This should be of value in phylogenetic studies of transmembrane proteins: for example, mitochondrial proteins have acquired a special importance in phylogenetics and are mostly transmembrane proteins. The improvement in fit to example data sets of our new model relative to less complex models of amino acid sequence evolution is statistically tested. To further illustrate the potential utility of our method, phylogeny estimation is performed on primate CCR5 receptor sequences, sequences of l and m subunits of the light reaction center in purple bacteria, guinea pig sequences with respect to lagomorph and rodent sequences of calcitonin receptor and K-substance receptor, and cetacean sequences of cytochrome b.