Mammalian and yeast 14-3-3 isoforms form distinct patterns of dimers in vivo

The 14-3-3 protein family associates with many proteins involved in intracellular signalling. In many cases, there is a distinct preference for a particular isoform(s) of 14-3-3. A specific repertoire of 14-3-3 dimer formation may therefore influence which of the interacting proteins could be brought together. We have analysed the pattern of dimer formation for two of the most abundant isoforms of 14-3-3, epsilon ( epsilon ) and gamma (gamma), following their stable expression. This revealed a distinct preference for particular dimer combinations that is largely independent of cellular conditions. gamma 14-3-3 occurred as homodimers and also formed heterodimers, mainly with epsilon 14-3-3 (In PC12 and Cos cells). The epsilon isoform formed heterodimers with 14-3-3 beta, gamma, zeta, and eta, but no homodimers were detected. The two 14-3-3 homologues, BMH1 and BMH2 from Saccharomyces cerevisiae, were mainly heterodimers.     

Evolutionary and structural analyses of heterodimeric proteins composed of subunits with same fold

Heterodimeric proteins with homologous subunits of same fold are involved in various biological processes. The objective of this study is to understand the evolution of structural and functional features of such heterodimers. Using a non‐redundant dataset of 70 such heterodimers of known 3D structure and an independent dataset of 173 heterodimers from yeast, we note that the mean sequence identity between interacting homologous subunits is only 23–24% suggesting that, generally, highly diverged paralogues assemble to form such a heterodimer. We also note that the functional roles of interacting subunits/domains are generally quite different. This suggests that, though the interacting subunits/domains are homologous, the high evolutionary divergence characterize their high functional divergence which contributes to a gross function for the heterodimer considered as a whole. The inverse relationship between sequence identity and RMSD of interacting homologues in heterodimers is not followed. We also addressed the question of formation of homodimers of the subunits of heterodimers by generating models of fictitious homodimers on the basis of the 3D structures of the heterodimers. Interaction energies associated with these homodimers suggests that, in overwhelming majority of the cases, such homodimers are unlikely to be stable. Majority of the homologues of heterodimers of known structures form heterodimers (51.8%) and a small proportion (14.6%) form homodimers. Comparison of 3D structures of heterodimers with homologous homodimers suggests that interfacial nature of residues is not well conserved. In over 90% of the cases we note that the interacting subunits of heterodimers are co‐localized in the cell. Proteins 2015; 83:1766–1786. © 2015 Wiley Periodicals, Inc.     

Structural basis of diversity and homodimerization specificity of zinc-finger-associated domains in Drosophila

In arthropods, zinc finger-associated domains (ZADs) are found at the N-termini of many DNA-binding proteins with tandem arrays of Cys2-His2 zinc fingers (ZAD-C2H2 proteins). ZAD-C2H2 proteins undergo fast evolutionary lineage-specific expansion and functional diversification. Here, we show that all ZADs from Drosophila melanogaster form homodimers, but only certain ZADs with high homology can also heterodimerize. CG2712, for example, is unable to heterodimerize with its paralog, the previously characterized insulator protein Zw5, with which it shares 46% homology. We obtained a crystal structure of CG2712 protein''s ZAD domain that, in spite of a low sequence homology, has similar spatial organization with the only known ZAD structure (from Grauzone protein). Steric clashes prevented the formation of heterodimers between Grauzone and CG2712 ZADs. Using detailed structural analysis, site-directed mutagenesis, and molecular dynamics simulations, we demonstrated that rapid evolutionary acquisition of interaction specificity was mediated by the more energy-favorable formation of homodimers in comparison to heterodimers, and that this specificity was achieved by multiple amino acid substitutions resulting in the formation or breaking of stabilizing interactions. We speculate that specific homodimerization of ZAD-C2H2 proteins is important for their architectural role in genome organization.     

Binding properties and evolution of homodimers in protein-protein interaction networks

We demonstrate that protein–protein interaction networks in several eukaryotic organisms contain significantly more self-interacting proteins than expected if such homodimers randomly appeared in the course of the evolution. We also show that on average homodimers have twice as many interaction partners than non-self-interacting proteins. More specifically, the likelihood of a protein to physically interact with itself was found to be proportional to the total number of its binding partners. These properties of dimers are in agreement with a phenomenological model, in which individual proteins differ from each other by the degree of their ‘stickiness’ or general propensity toward interaction with other proteins including oneself. A duplication of self-interacting proteins creates a pair of paralogous proteins interacting with each other. We show that such pairs occur more frequently than could be explained by pure chance alone. Similar to homodimers, proteins involved in heterodimers with their paralogs on average have twice as many interacting partners than the rest of the network. The likelihood of a pair of paralogous proteins to interact with each other was also shown to decrease with their sequence similarity. This points to the conclusion that most of interactions between paralogs are inherited from ancestral homodimeric proteins, rather than established de novo after duplication. We finally discuss possible implications of our empirical observations from functional and evolutionary standpoints.     

Why Have Small Multidrug Resistance Proteins Not Evolved into Fused,Internally Duplicated Structures?

The increasing number of solved membrane protein structures has led to the recognition of a common feature in a large fraction of the small-molecule transporters: inverted repeat structures, formed by two fused homologous membrane domains with opposite orientation in the membrane. An evolutionary pathway in which the ancestral state is a single gene encoding a dual-topology membrane protein capable of forming antiparallel homodimers has been posited. A gene duplication event enables the evolution of two oppositely orientated proteins that form antiparallel heterodimers. Finally, fusion of the two genes generates an internally duplicated transporter with two oppositely orientated membrane domains. Strikingly, however, in the small multidrug resistance (SMR) family of transporters, no fused, internally duplicated proteins have been found to date. Here, we have analyzed fused versions of the dual-topology transporter EmrE, a member of the SMR family, by blue-native PAGE and in vivo activity measurements. We find that fused constructs give rise to both intramolecular inverted repeat structures and competing intermolecular dimers of varying activity. The formation of several intramolecularly and intermolecularly paired species indicates that a gene fusion event may lower the overall amount of active protein, possibly explaining the apparent absence of fused SMR proteins in nature.     

Molecular Basis for Bcl-2 Homology 3 Domain Recognition in the Bcl-2 Protein Family: IDENTIFICATION OF CONSERVED HOT SPOT INTERACTIONS*

The proteins of the Bcl-2 family are important regulators of apoptosis, or programmed cell death. These proteins regulate this fundamental biological process via the formation of heterodimers involving both pro- and anti-apoptotic family members. Disruption of the balance between anti- and pro-apoptotic Bcl-2 proteins is the cause of numerous pathologies. Bcl-xl, an anti-apoptotic protein of this family, is known to form heterodimers with multiple pro-apoptotic proteins, such as Bad, Bim, Bak, and Bid. To elucidate the molecular basis of this recognition process, we used molecular dynamics simulations coupled with the Molecular Mechanics/Poisson-Boltzmann Surface Area approach to identify the amino acids that make significant energetic contributions to the binding free energy of four complexes formed between Bcl-xl and pro-apoptotic Bcl-2 homology 3 peptides. A fifth protein-peptide complex composed of another anti-apoptotic protein, Bcl-w, in complex with the peptide from Bim was also studied. The results identified amino acids of both the anti-apoptotic proteins as well as the Bcl-2 homology 3 (BH3) domains of the pro-apoptotic proteins that make strong, recurrent interactions in the protein complexes. The calculations show that the two anti-apoptotic proteins, Bcl-xl and Bcl-w, share a similar recognition mechanism. Our results provide insight into the molecular basis for the promiscuous nature of this molecular recognition process by members of the Bcl-2 protein family. These amino acids could be targeted in the design of new mimetics that serve as scaffolds for new antitumoral molecules.     

Live cell FRET microscopy: homo- and heterodimerization of two human peroxisomal ABC transporters, the adrenoleukodystrophy protein (ALDP, ABCD1) and PMP70 (ABCD3)

The adrenoleukodystrophy protein (ALDP) and the 70-kDa peroxisomal membrane protein (PMP70) are half-ATP-binding cassette (ABC) transporters in the mammalian peroxisome membrane. Mutations in the gene encoding ALDP result in a devastating neurodegenerative disorder, X-linked adrenoleukodystrophy (X-ALD) that is associated with elevated levels of very long chain fatty acids because of impaired peroxisomal beta-oxidation. The interactions of peroxisomal ABC transporters, their role in the peroxisomal membrane, and their functions in disease pathogenesis are poorly understood. Studies on ABC transporters revealed that half-transporters have to dimerize to gain functionality. So far, conflicting observations are described for ALDP. By the use of in vitro methods (yeast two-hybrid and immunoprecipitation assays) on the one hand, it was shown that ALDP can form homodimers as well as heterodimers with PMP70 and ALDR, while on the other hand, it was demonstrated that ALDP and PMP70 exclusively homodimerize. To circumvent the problems of artificial interactions due to biochemical sample preparation in vitro, we investigated protein-protein interaction of ALDP in its physiological environment by FRET microscopy in intact living cells. The statistical relevance of FRET data was determined in two different ways using probability distribution shift analysis and Kolmogorov-Smirnov statistics. We demonstrate in vivo that ALDP and PMP70 form homodimers as well as ALDP/PMP70 heterodimers where ALDP homodimers predominate. Using C-terminal deletion constructs of ALDP, we demonstrate that the last 87 C-terminal amino acids harbor the most important protein domain mediating these interactions, and that the N-terminal transmembrane region of ALDP has an additional stabilization effect on ALDP homodimers. Loss of ALDP homo- or heterodimerization is highly relevant for understanding the disease mechanisms of X-ALD.     

Free leucine dissociates homo- and heterodimers formed between proteins containing leucine heptad repeats

A highly specific method for the dissociation of protein dimers has been developed. The method involves exposure of the dimers to free leucine at a concentration ranging between 3 and 10 mM. Using this method it has been possible to dissociate goat uterine oestrogen receptor homodimers, heterodimers formed between the non-activated oestrogen receptor (naER) and the oestrogen receptor activation factor (E-RAF) of the goat uterus, c-jun homodimers derived from bovine bone marrow and also glucocorticoid receptor homodimers isolated from rat liver cytosol. The pattern of dimer dissociation by leucine clearly differentiates two classes of proteins. The first is represented by steroid hormone receptors where dimerization is apparently contributed by both coiled-coil dimerization interfaces and the conserved heptad repeats of leucine. The second is represented by oncoproteins like c-fos and c-jun which dimerize through the exclusive involvement of leucine zippers. The patterns of dissociation of these two groups of proteins from the concerned affinity columns are distinctly different. This indicates a possibility that the elution pattern may be used as a yardstick to determine whether two proteins dimerize through the exclusive involvement of leucine zippers or whether coiled-coil interfaces are also involved in the dimerization process.     

Modes of interaction among yeast Nej1, Lif1 and Dnl4 proteins and comparison to human XLF, XRCC4 and Lig4

The nonhomologous end joining (NHEJ) pathway of double-strand break repair depends on DNA ligase IV and its interacting partner protein XRCC4 (Lif1 in yeast). A third yeast protein, Nej1, interacts with Lif1 and supports NHEJ, similar to the distantly related mammalian Nej1 orthologue XLF (also known as Cernunnos). XRCC4/Lif1 and XLF/Nej1 are themselves related and likely fold into similar coiled-coil structures, which suggests many possible modes of interaction between these proteins. Using yeast two-hybrid and co-precipitation methods we examined these interactions and the protein domains required to support them. Results suggest that stable coiled-coil homodimers are a predominant form of XLF/Nej1, just as for XRCC4/Lif1, but that similar heterodimers are not. XLF-XRCC4 and Nej1-Lif1 interactions were instead mediated independently of the coiled coil, and by different regions of XLF and Nej1. Specifically, the globular head of XRCC4/Lif1 interacted with N- and C-terminal domains of XLF and Nej1, respectively. Direct interactions between XLF/Nej1 and DNA ligase IV were also observed, but again appeared qualitatively different than the stable coiled-coil-mediated interaction between XRCC4/Lif1 and DNA ligase IV. The implications of these findings for DNA ligase IV function are considered in light of the evolutionary pattern in the XLF/XRCC4 and XLF/Nej1 family.     

The Subunit Interfaces of Weakly Associated Homodimeric Proteins

We analyzed subunit interfaces in 315 homodimers with an X-ray structure in the Protein Data Bank, validated by checking the literature for data that indicate that the proteins are dimeric in solution and that, in the case of the “weak” dimers, the homodimer is in equilibrium with the monomer. The interfaces of the 42 weak dimers, which are smaller by a factor of 2.4 on average than in the remainder of the set, are comparable in size with antibody-antigen or protease-inhibitor interfaces. Nevertheless, they are more hydrophobic than in the average transient protein-protein complex and similar in amino acid composition to the other homodimer interfaces. The mean numbers of interface hydrogen bonds and hydration water molecules per unit area are also similar in homodimers and transient complexes. Parameters related to the atomic packing suggest that many of the weak dimer interfaces are loosely packed, and we suggest that this contributes to their low stability. To evaluate the evolutionary selection pressure on interface residues, we calculated the Shannon entropy of homologous amino acid sequences at 60% sequence identity. In 93% of the homodimers, the interface residues are better conserved than the residues on the protein surface. The weak dimers display the same high degree of interface conservation as other homodimers, but their homologs may be heterodimers as well as homodimers. Their interfaces may be good models in terms of their size, composition, and evolutionary conservation for the labile subunit contacts that allow protein assemblies to share and exchange components, allosteric proteins to undergo quaternary structure transitions, and molecular machines to operate in the cell.     

Evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase family of the enolase superfamily

Understanding how proteins evolve to provide both exquisite specificity and proficient activity is a fundamental problem in biology that has implications for protein function prediction and protein engineering. To study this problem, we analyzed the evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase (OSBS/NAAAR) family, part of the mechanistically diverse enolase superfamily. Although all characterized members of the family catalyze the OSBS reaction, this family is extraordinarily divergent, with some members sharing <15% identity. In addition, a member of this family, Amycolatopsis OSBS/NAAAR, is promiscuous, catalyzing both dehydration and racemization. Although the OSBS/NAAAR family appears to have a single evolutionary origin, no sequence or structural motifs unique to this family could be identified; all residues conserved in the family are also found in enolase superfamily members that have different functions. Based on their species distribution, several uncharacterized proteins similar to Amycolatopsis OSBS/NAAAR appear to have been transmitted by lateral gene transfer. Like Amycolatopsis OSBS/NAAAR, these might have additional or alternative functions to OSBS because many are from organisms lacking the pathway in which OSBS is an intermediate. In addition to functional differences, the OSBS/NAAAR family exhibits surprising structural variations, including large differences in orientation between the two domains. These results offer several insights into protein evolution. First, orthologous proteins can exhibit significant structural variation, and specificity can be maintained with little conservation of ligand-contacting residues. Second, the discovery of a set of proteins similar to Amycolatopsis OSBS/NAAAR supports the hypothesis that new protein functions evolve through promiscuous intermediates. Finally, a combination of evolutionary, structural, and sequence analyses identified characteristics that might prime proteins, such as Amycolatopsis OSBS/NAAAR, for the evolution of new activities.     

Purification, cloning, and RXR identity of the HeLa cell factor with which RAR or TR heterodimerizes to bind target sequences efficiently.

We have purified and cloned a HeLa cell nuclear protein that strongly stimulates binding of retinoic acid and thyroid hormone receptors (RARs and TRs) to response elements. The purified protein is a human retinoid X receptor beta (hRXR beta). Three murine members of the RXR family (mRXR alpha, beta, and gamma) have also been cloned, and their interactions with RARs and TRs have been investigated. Under conditions where RAR, RXR, and TR bound poorly as homodimers to various response elements, strongly cooperative RAR-RXR and TR-RXR binding was observed. The binding efficiency was dependent on the sequence, relative orientation, and spacing of the repeated motifs of response elements. We show also that unstable RAR-RXR heterodimers were formed in solution, and that C-terminal sequences and the DNA-binding domains of both receptors were required for efficient formation of stable heterodimers on response elements. These findings suggest a convergence of the signaling pathways of some members of the nuclear receptor superfamily.     

Investigation at residue level of the early steps during the assembly of two proteins into supramolecular objects

Understanding the driving forces governing protein assembly requires the characterization of interactions at molecular level. We focus on two homologous oppositely charged proteins, lysozyme and α-lactalbumin, which can assemble into microspheres. The assembly early steps were characterized through the identification of interacting surfaces monitored at residue level by NMR chemical shift perturbations by titrating one (15)N-labeled protein with its unlabeled partner. While α-lactalbumin has a narrow interacting site, lysozyme has interacting sites scattered on a broad surface. The further assembly of these rather unspecific heterodimers into tetramers leads to the establishment of well-defined interaction sites. Within the tetramers, most of the electrostatic charge patches on the protein surfaces are shielded. Then, hydrophobic interactions, which are possible because α-lactalbumin is in a partially folded state, become preponderant, leading to the formation of larger oligomers. This approach will be particularly useful for rationalizing the design of protein assemblies as nanoscale devices.     

Evolution of antiparallel two-domain membrane proteins: tracing multiple gene duplication events in the DUF606 family

X-ray crystallography has revealed that many integral membrane proteins consist of two domains with a similar fold but opposite (antiparallel) orientation in the membrane. The proteins are believed to have evolved by gene duplication and gene fusion events from a dual topology ancestral membrane protein, that adapted both orientations in the membrane and formed antiparallel homodimers. Here, we present a detailed analysis of the DUF606 family of bacterial membrane proteins that contains the entire collection of intermediate states of such an evolutionary pathway: single genes that would code for dual topology homodimeric proteins, paired genes coding for homologous proteins with a fixed but opposite orientation in the membrane that would form heterodimers, and fused genes that encode antiparallel two-domain fusion proteins. Two types of paired genes can be discriminated corresponding to the order in which the genes coding for the two oppositely oriented proteins occur in the operon. On the protein level, the heterodimers resulting from the two types of gene pairs are indistinguishable. In contrast, two types of fused genes corresponding to the two possible orders in which the oppositely oriented domains are present in the encoded proteins, do result in discernible types of proteins. The large number of genetic and protein states in the DUF606 family allowed for a detailed phylogenic analysis that revealed a total of nine independent duplication events in the DUF606 family, five of which resulted in paired genes, and four resulted in fused genes. Noticeably, there was no evidence for a sequential mechanism in which fusions evolve from a pair of genes. Rather, an evolutionary mechanism is proposed by which antiparallel two-domain proteins are the direct result of a gene duplication event. Combining the phylogeny of proteins and hosting microorganisms allowed for a reconstruction of the evolutionary pathway.     

Clustering of adhesion receptors following exposure of insect blood cells to foreign surfaces

Cell-mediated immune responses of insects involve interactions of two main classes of blood cells (hemocytes) known as granular cells and plasmatocytes. In response to a foreign surface, these hemocytes suddenly transform from circulating, non-adherent cells to cells that interact and adhere to each other and the foreign surface. This report presents evidence that during this adhesive transformation the extracellular matrix (ECM) proteins lacunin and a ligand for peanut agglutinin (PNA) lectin are released by granular cells and bind to surfaces of both granular cells and plasmatocytes. ECM protein co-localizes on cell surfaces with the adhesive receptors integrin and neuroglian, a member of the immunoglobulin superfamily. The ECM protein(s) secreted by granular cells are hypothesized to interact with adhesion receptors such as neuroglian and integrin by cross linking and clustering them on hemocyte surfaces. This clustering of receptors is known to enhance the adhesiveness (avidity) of interacting mammalian immune cells. The formation of ring-shaped clusters of these adhesion receptors on surfaces of insect immune cells represents an evolutionary antecedent of the mammalian immunological synapse.     

The Analysis of Participation of Individual Proteins in the Protein Interactome Formation

It becomes increasingly clear that most proteins of living systems exist as components of various protein complexes rather than individual molecules. The use of various proteomic techniques significantly extended our knowledge not only about functioning of individual complexes but also formed a basis for systemic analysis of protein-protein interactions. In this study gel-filtration chromatography accompanied by mass spectrometry was used for the interactome analysis of human liver proteins. In six fractions (with average molecular masses of 45 kDa, 60 kDa, 85 kDa, 150 kDa, 250 kDa, and 440 kDa) 797 proteins were identified. In dependence of their distribution profiles in the fractions, these proteins could be subdivided into four groups: (1) single monomeric proteins that are not involved in formation of stable protein complexes; (2) proteins existing as homodimers or heterodimers with comparable partners; (3) proteins that are partially exist as monomers and partially as components of protein complexes; (4) proteins that do not exist in the monomolecular state, but also exist within protein complexes containing three or more subunits. Application of this approach to known isatin-binding proteins resulted in identification of proteins involved in formation of the homo- and heterodimers and mixed protein complexes.     

Different Epidermal Growth Factor (EGF) Receptor Ligands Show Distinct Kinetics and Biased or Partial Agonism for Homodimer and Heterodimer Formation

The EGF receptor has seven different cognate ligands. Previous work has shown that these different ligands are capable of inducing different biological effects, even in the same cell. To begin to understand the molecular basis for this variation, we used luciferase fragment complementation to measure ligand-induced dimer formation and radioligand binding to study the effect of the ligands on subunit-subunit interactions in EGF receptor (EGFR) homodimers and EGFR/ErbB2 heterodimers. In luciferase fragment complementation imaging studies, amphiregulin (AREG) functioned as a partial agonist, inducing only about half as much total dimerization as the other three ligands. However, unlike the other ligands, AREG showed biphasic kinetics for dimer formation, suggesting that its path for EGF receptor activation involves binding to both monomers and preformed dimers. EGF, TGFα, and betacellulin (BTC) appear to mainly stimulate receptor activation through binding to and dimerization of receptor monomers. In radioligand binding assays, EGF and TGFα exhibited increased affinity for EGFR/ErbB2 heterodimers compared with EGFR homodimers. By contrast, BTC and AREG showed a similar affinity for both dimers. Thus, EGF and TGFα are biased agonists, whereas BTC and AREG are balanced agonists with respect to selectivity of dimer formation. These data suggest that the differences in biological response to different EGF receptor ligands may result from partial agonism for dimer formation, differences in the kinetic pathway utilized to generate activated receptor dimers, and biases in the formation of heterodimers versus homodimers.