Integration of protein motions with molecular networks reveals different mechanisms for permanent and transient interactions

The integration of molecular networks with other types of data, such as changing levels of gene expression or protein-structural features, can provide richer information about interactions than the simple node-and-edge representations commonly used in the network community. For example, the mapping of 3D-structural data onto networks enables classification of proteins into singlish- or multi-interface hubs (depending on whether they have >2 interfaces). Similarly, interactions can be classified as permanent or transient, depending on whether their interface is used by only one or by multiple partners. Here, we incorporate an additional dimension into molecular networks: dynamic conformational changes. We parse the entire PDB structural databank for alternate conformations of proteins and map these onto the protein interaction network, to compile a first version of the Dynamic Structural Interaction Network (DynaSIN). We make this network available as a readily downloadable resource file, and we then use it to address a variety of downstream questions. In particular, we show that multi-interface hubs display a greater degree of conformational change than do singlish-interface ones; thus, they show more plasticity which perhaps enables them to utilize more interfaces for interactions. We also find that transient associations involve smaller conformational changes than permanent ones. Although this may appear counterintuitive, it is understandable in the following framework: as proteins involved in transient interactions shuttle between interchangeable associations, they interact with domains that are similar to each other and so do not require drastic structural changes for their activity. We provide evidence for this hypothesis through showing that interfaces involved in transient interactions bind fewer classes of domains than those in a control set.     

Multiple protein–protein interactions converging on the Prp38 protein during activation of the human spliceosome

Spliceosomal Prp38 proteins contain a conserved amino-terminal domain, but only higher eukaryotic orthologs also harbor a carboxy-terminal RS domain, a hallmark of splicing regulatory SR proteins. We show by crystal structure analysis that the amino-terminal domain of human Prp38 is organized around three pairs of antiparallel α-helices and lacks similarities to RNA-binding domains found in canonical SR proteins. Instead, yeast two-hybrid analyses suggest that the amino-terminal domain is a versatile protein–protein interaction hub that possibly binds 12 other spliceosomal proteins, most of which are recruited at the same stage as Prp38. By quantitative, alanine surface-scanning two-hybrid screens and biochemical analyses we delineated four distinct interfaces on the Prp38 amino-terminal domain. In vitro interaction assays using recombinant proteins showed that Prp38 can bind at least two proteins simultaneously via two different interfaces. Addition of excess Prp38 amino-terminal domain to in vitro splicing assays, but not of an interaction-deficient mutant, stalled splicing at a precatalytic stage. Our results show that human Prp38 is an unusual SR protein, whose amino-terminal domain is a multi-interface protein–protein interaction platform that might organize the relative positioning of other proteins during splicing.     

Interaction interfaces of protein domains are not topologically equivalent across families within superfamilies: Implications for metabolic and signaling pathways

Using a data set of aligned protein domain superfamilies of known three-dimensional structure, we compared the location of interdomain interfaces on the tertiary folds between members of distantly related protein domain superfamilies. The data set analyzed is comprised of interdomain interfaces, with domains occurring within a polypeptide chain and those between two polypeptide chains. We observe that, in general, the interfaces between protein domains are formed entirely in different locations on the tertiary folds in such pairs. This variation in the location of interface happens in protein domains involved in a wide range of functions, such as enzymes, adapters, and domains that bind protein ligands, or cofactors. While basic biochemical functionality is preserved at the domain superfamily level, the effect of biochemical function on protein assemblies is different in these protein domains related by superfamily. The divergence between proteins, in most cases, is coupled with domain recruitment, with different modes of interaction with the recruited domain. This is in complete contrast to the observation that in closely related homologous protein domains, almost always the interaction interfaces are topologically equivalent. In a small subset of interacting domains within proteins related by remote homology, we observe that the relative positioning of domains with respect to one another is preserved. Based on the analysis of multidomain proteins of known or unknown structure, we suggest that variation in protein-protein interactions in members within a superfamily could serve as diverging points in otherwise parallel metabolic or signaling pathways. We discuss a few representative cases of diverging pathways involving domains in a superfamily.     

Accurate prediction of interfacial residues in two‐domain proteins using evolutionary information: Implications for three‐dimensional modeling

With the preponderance of multidomain proteins in eukaryotic genomes, it is essential to recognize the constituent domains and their functions. Often function involves communications across the domain interfaces, and the knowledge of the interacting sites is essential to our understanding of the structure–function relationship. Using evolutionary information extracted from homologous domains in at least two diverse domain architectures (single and multidomain), we predict the interface residues corresponding to domains from the two‐domain proteins. We also use information from the three‐dimensional structures of individual domains of two‐domain proteins to train naïve Bayes classifier model to predict the interfacial residues. Our predictions are highly accurate (～85%) and specific (～95%) to the domain–domain interfaces. This method is specific to multidomain proteins which contain domains in at least more than one protein architectural context. Using predicted residues to constrain domain–domain interaction, rigid‐body docking was able to provide us with accurate full‐length protein structures with correct orientation of domains. We believe that these results can be of considerable interest toward rational protein and interaction design, apart from providing us with valuable information on the nature of interactions. Proteins 2014; 82:1219–1234. © 2013 Wiley Periodicals, Inc.     

Structure of the Hemoglobin-IsdH Complex Reveals the Molecular Basis of Iron Capture by Staphylococcus aureus

Staphylococcus aureus causes life-threatening disease in humans. The S. aureus surface protein iron-regulated surface determinant H (IsdH) binds to mammalian hemoglobin (Hb) and extracts heme as a source of iron, which is an essential nutrient for the bacteria. However, the process of heme transfer from Hb is poorly understood. We have determined the structure of IsdH bound to human Hb by x-ray crystallography at 4.2 Å resolution, revealing the structural basis for heme transfer. One IsdH molecule is bound to each α and β Hb subunit, suggesting that the receptor acquires iron from both chains by a similar mechanism. Remarkably, two near iron transporter (NEAT) domains in IsdH perform very different functions. An N-terminal NEAT domain binds α/β globin through a site distant from the globin heme pocket and, via an intervening structural domain, positions the C-terminal heme-binding NEAT domain perfectly for heme transfer. These data, together with a 2.3 Å resolution crystal structure of the isolated N-terminal domain bound to Hb and small-angle x-ray scattering of free IsdH, reveal how multiple domains of IsdH cooperate to strip heme from Hb. Many bacterial pathogens obtain iron from human hemoglobin using proteins that contain multiple NEAT domains and other domains whose functions are poorly understood. Our results suggest that, rather than acting as isolated units, NEAT domains may be integrated into higher order architectures that employ multiple interaction interfaces to efficiently extract heme from host proteins.     

Light-induced Changes in the Dimerization Interface of Bacteriophytochromes

Phytochromes are dimeric photoreceptor proteins that sense red light levels in plants, fungi, and bacteria. The proteins are structurally divided into a light-sensing photosensory module consisting of PAS, GAF, and PHY domains and a signaling output module, which in bacteriophytochromes typically is a histidine kinase (HK) domain. Existing structural data suggest that two dimerization interfaces exist between the GAF and HK domains, but their functional roles remain unclear. Using mutational, biochemical, and computational analyses of the Deinococcus radiodurans phytochrome, we demonstrate that two dimerization interfaces between sister GAF and HK domains stabilize the dimer with approximately equal contributions. The existence of both dimerization interfaces is critical for thermal reversion back to the resting state. We also find that a mutant in which the interactions between the GAF domains were removed monomerizes under red light. This implies that the interactions between the HK domains are significantly altered by photoconversion. The results suggest functional importance of the dimerization interfaces in bacteriophytochromes.     

Regulation of Response Regulator Autophosphorylation through Interdomain Contacts

DNA-binding response regulators (RRs) of the OmpR/PhoB subfamily alternate between inactive and active conformational states, with the latter having enhanced DNA-binding affinity. Phosphorylation of an aspartate residue in the receiver domain, usually via phosphotransfer from a cognate histidine kinase, stabilizes the active conformation. Many of the available structures of inactive OmpR/PhoB family proteins exhibit extensive interfaces between the N-terminal receiver and C-terminal DNA-binding domains. These interfaces invariably involve the α4-β5-α5 face of the receiver domain, the locus of the largest differences between inactive and active conformations and the surface that mediates dimerization of receiver domains in the active state. Structures of receiver domain dimers of DrrB, DrrD, and MtrA have been determined, and phosphorylation kinetics were analyzed. Analysis of phosphotransfer from small molecule phosphodonors has revealed large differences in autophosphorylation rates among OmpR/PhoB RRs. RRs with substantial domain interfaces exhibit slow rates of phosphorylation. Rates are greatly increased in isolated receiver domain constructs. Such differences are not observed between autophosphorylation rates of full-length and isolated receiver domains of a RR that lacks interdomain interfaces, and they are not observed in histidine kinase-mediated phosphotransfer. These findings suggest that domain interfaces restrict receiver domain conformational dynamics, stabilizing an inactive conformation that is catalytically incompetent for phosphotransfer from small molecule phosphodonors. Inhibition of phosphotransfer by domain interfaces provides an explanation for the observation that some RRs cannot be phosphorylated by small molecule phosphodonors in vitro and provides a potential mechanism for insulating some RRs from small molecule-mediated phosphorylation in vivo.     

The many faces of protein-protein interactions: A compendium of interface geometry

A systematic classification of protein-protein interfaces is a valuable resource for understanding the principles of molecular recognition and for modelling protein complexes. Here, we present a classification of domain interfaces according to their geometry. Our new algorithm uses a hybrid approach of both sequential and structural features. The accuracy is evaluated on a hand-curated dataset of 416 interfaces. Our hybrid procedure achieves 83% precision and 95% recall, which improves the earlier sequence-based method by 5% on both terms. We classify virtually all domain interfaces of known structure, which results in nearly 6,000 distinct types of interfaces. In 40% of the cases, the interacting domain families associate in multiple orientations, suggesting that all the possible binding orientations need to be explored for modelling multidomain proteins and protein complexes. In general, hub proteins are shown to use distinct surface regions (multiple faces) for interactions with different partners. Our classification provides a convenient framework to query genuine gene fusion, which conserves binding orientation in both fused and separate forms. The result suggests that the binding orientations are not conserved in at least one-third of the gene fusion cases detected by a conventional sequence similarity search. We show that any evolutionary analysis on interfaces can be skewed by multiple binding orientations and multiple interaction partners. The taxonomic distribution of interface types suggests that ancient interfaces common to the three major kingdoms of life are enriched by symmetric homodimers. The classification results are online at http://www.scoppi.org.     

In the Multi-domain Protein Adenylate Kinase,Domain Insertion Facilitates Cooperative Folding while Accommodating Function at Domain Interfaces

Having multiple domains in proteins can lead to partial folding and increased aggregation. Folding cooperativity, the all or nothing folding of a protein, can reduce this aggregation propensity. In agreement with bulk experiments, a coarse-grained structure-based model of the three-domain protein, E. coli Adenylate kinase (AKE), folds cooperatively. Domain interfaces have previously been implicated in the cooperative folding of multi-domain proteins. To understand their role in AKE folding, we computationally create mutants with deleted inter-domain interfaces and simulate their folding. We find that inter-domain interfaces play a minor role in the folding cooperativity of AKE. On further analysis, we find that unlike other multi-domain proteins whose folding has been studied, the domains of AKE are not singly-linked. Two of its domains have two linkers to the third one, i.e., they are inserted into the third one. We use circular permutation to modify AKE chain-connectivity and convert inserted-domains into singly-linked domains. We find that domain insertion in AKE achieves the following: (1) It facilitates folding cooperativity even when domains have different stabilities. Insertion constrains the N- and C-termini of inserted domains and stabilizes their folded states. Therefore, domains that perform conformational transitions can be smaller with fewer stabilizing interactions. (2) Inter-domain interactions are not needed to promote folding cooperativity and can be tuned for function. In AKE, these interactions help promote conformational dynamics limited catalysis. Finally, using structural bioinformatics, we suggest that domain insertion may also facilitate the cooperative folding of other multi-domain proteins.     

Sequence and structural determinants of strand swapping in cadherin domains: do all cadherins bind through the same adhesive interface?

Cadherins are cell surface adhesion proteins important for tissue development and integrity. Type I and type II, or classical, cadherins form adhesive dimers via an interface formed through the exchange, or “swapping”, of the N-terminal β-strands from their membrane-distal EC1 domains. Here, we ask which sequence and structural features in EC1 domains are responsible for β-strand swapping and whether members of other cadherin families form similar strand-swapped binding interfaces. We created a comprehensive database of multiple alignments of each type of cadherin domain. We used the known three-dimensional structures of classical cadherins to identify conserved positions in multiple sequence alignments that appear to be crucial determinants of the cadherin domain structure. We identified features that are unique to EC1 domains. On the basis of our analysis, we conclude that all cadherin domains have very similar overall folds but, with the exception of classical and desmosomal cadherin EC1 domains, most of them do not appear to bind through a strand-swapping mechanism. Thus, non-classical cadherins that function in adhesion are likely to use different protein-protein interaction interfaces. Our results have implications for the evolution of molecular mechanisms of cadherin-mediated adhesion in vertebrates.     

The structural basis of mode of activation and functional diversity: a case study with HtrA family of serine proteases

HtrA (High temperature requirement protease A) proteins that are primarily involved in protein quality control belong to a family of serine proteases conserved from bacteria to humans. HtrAs are oligomeric proteins that share a common trimeric pyramidal architecture where each monomer comprises a serine protease domain and one or two PDZ domains. Although the overall structural integrity is well maintained and they exhibit similar mechanism of activation, subtle conformational changes and structural plasticity especially in the flexible loop regions and domain interfaces lead to differences in their active site conformation and hence in their specificity and functions.     

Multi-domain proteins in the three kingdoms of life: orphan domains and other unassigned regions

Comparative studies of the proteomes from different organisms have provided valuable information about protein domain distribution in the kingdoms of life. Earlier studies have been limited by the fact that only about 50% of the proteomes could be matched to a domain. Here, we have extended these studies by including less well-defined domain definitions, Pfam-B and clustered domains, MAS, in addition to Pfam-A and SCOP domains. It was found that a significant fraction of these domain families are homologous to Pfam-A or SCOP domains. Further, we show that all regions that do not match a Pfam-A or SCOP domain contain a significantly higher fraction of disordered structure. These unstructured regions may be contained within orphan domains or function as linkers between structured domains. Using several different definitions we have re-estimated the number of multi-domain proteins in different organisms and found that several methods all predict that eukaryotes have approximately 65% multi-domain proteins, while the prokaryotes consist of approximately 40% multi-domain proteins. However, these numbers are strongly dependent on the exact choice of cut-off for domains in unassigned regions. In conclusion, all eukaryotes have similar fractions of multi-domain proteins and disorder, whereas a high fraction of repeating domain is distinguished only in multicellular eukaryotes. This implies a role for repeats in cell-cell contacts while the other two features are important for intracellular functions.     

The generation of new protein functions by the combination of domains

During evolution, many new proteins have been formed by the process of gene duplication and combination. The genes involved in this process usually code for whole domains. Small proteins contain one domain; medium and large proteins contain two or more domains. We have compared homologous domains that occur in both one-domain proteins and multidomain proteins. We have determined (1) how the functions of the individual domains in the multidomain proteins combine to produce their overall functions and (2) the extent to which these functions are similar to those in the one-domain homologs. We describe how domain combinations increase the specificity of enzymes; act as links between domains that have functional roles; regulate activity; combine within one chain functions that can act either independently, in concert or in new contexts; and provide the structural framework for the evolution of entirely new functions.     

真菌中锌簇蛋白的结构和功能

锌簇家族蛋白即Zn₂Cys₆类锌指蛋白,是真菌中特有的一类蛋白,它们属于转录因子类,广泛参与真菌中初级和次级代谢、胁迫应答和细胞分裂等生命活动的调控。锌簇蛋白主要包括N端的DNA结合结构域、中间的调节结构域和C端的酸性区域,其中DNA结合结构域包含锌指基序并负责结合靶基因的启动子。目前已经解析了多个锌簇家族转录因子DNA结合结构域的三维结构,并发现该家族中一些蛋白能够参与调控多个基因的表达,但缺乏对其结构、动力学和功能关系的全面分析。本文综合分析了不同锌簇蛋白与DNA结合的结构特征,总结其结构域与功能的关系,指出锌簇蛋白研究的重要方向,旨在为锌簇家族蛋白的深入研究提供思路。     

Understanding the role of domain–domain linkers in the spatial orientation of domains in multi-domain proteins

Inter-domain linkers (IDLs)’ bridge flanking domains and support inter-domain communication in multi-domain proteins. Their sequence and conformational preferences enable them to carry out varied functions. They also provide sufficient flexibility to facilitate domain motions and, in conjunction with the interacting interfaces, they also regulate the inter-domain geometry (IDG). In spite of the basic intuitive understanding of the inter-domain orientations with respect to linker conformations and interfaces, we still do not entirely understand the precise relationship among the three. We show that IDG is evolutionarily well conserved and is constrained by the domain–domain interface interactions. The IDLs modulate the interactions by varying their lengths, conformations and local structure, thereby affecting the overall IDG. Results of our analysis provide guidelines in modelling of multi-domain proteins from the tertiary structures of constituent domain components.     

Crystal structures of Nova-1 and Nova-2 K-homology RNA-binding domains.

BACKGROUND: Nova-1 and Nova-2 are related neuronal proteins that were initially cloned using antisera obtained from patients with the autoimmune neurological disease paraneoplastic opsoclonus-myoclonus ataxia (POMA). Both of these disease gene products contain three RNA-binding motifs known as K-homology or KH domains, and their RNA ligands have been identified via binding-site selection experiments. The KH motif structure has been determined previously using NMR spectroscopy, but not using X-ray crystallography. Many proteins contain more than one KH domain, yet there is no published structural information regarding the behavior of such multimers. RESULTS: We have obtained the first X-ray crystallographic structures of KH-domain-containing proteins. Structures of the third KH domains (KH3) of Nova-1 and Nova-2 were determined by multiple isomorphous replacement and molecular replacement at 2.6 A and 2.0 A, respectively. These highly similar RNA-binding motifs form a compact protease-resistant domain resembling an open-faced sandwich, consisting of a three-stranded antiparallel beta sheet topped by three alpha helices. In both Nova crystals, the lattice is composed of symmetric tetramers of KH3 domains that are created by two dimer interfaces. CONCLUSIONS: The crystal structures of both Nova KH3 domains are similar to the previously determined NMR structures. The most significant differences among the KH domains involve changes in the positioning of one or more of the alpha helices with respect to the betasheet, particularly in the NMR structure of the KH1 domain of the Fragile X disease protein FMR-1. Loop regions in the KH domains are clearly visible in the crystal structure, unlike the NMR structures, revealing the conformation of the invariant Gly-X-X-Gly segment that is thought to participate in RNA-binding and of the variable region. The tetrameric arrangements of the Nova KH3 domains provide insights into how KH domains may interact with each other in proteins containing multiple KH motifs.     

SnapDRAGON: a method to delineate protein structural domains from sequence data

We describe a method to identify protein domain boundaries from sequence information alone based on the assumption that hydrophobic residues cluster together in space. SnapDRAGON is a suite of programs developed to predict domain boundaries based on the consistency observed in a set of alternative ab initio three-dimensional (3D) models generated for a given protein multiple sequence alignment. This is achieved by running a distance geometry-based folding technique in conjunction with a 3D-domain assignment algorithm. The overall accuracy of our method in predicting the number of domains for a non-redundant data set of 414 multiple alignments, representing 185 single and 231 multiple-domain proteins, is 72.4 %. Using domain linker regions observed in the tertiary structures associated with each query alignment as the standard of truth, inter-domain boundary positions are delineated with an accuracy of 63.9 % for proteins comprising continuous domains only, and 35.4 % for proteins with discontinuous domains. Overall, domain boundaries are delineated with an accuracy of 51.8 %. The prediction accuracy values are independent of the pair-wise sequence similarities within each of the alignments. These results demonstrate the capability of our method to delineate domains in protein sequences associated with a wide variety of structural domain organisation.     

Functional fingerprints of folds: evidence for correlated structure-function evolution

Using structural similarity clustering of protein domains: protein domain universe graph (PDUG), and a hierarchical functional annotation: gene ontology (GO) as two evolutionary lenses, we find that each structural cluster (domain fold) exhibits a distribution of functions that is unique to it. These functional distributions are functional fingerprints that are specific to characteristic structural clusters and vary from cluster to cluster. Furthermore, as structural similarity threshold for domain clustering in the PDUG is relaxed we observe an influx of earlier-diverged domains into clusters. These domains join clusters without destroying the functional fingerprint. These results can be understood in light of a divergent evolution scenario that posits correlated divergence of structural and functional traits in protein domains from one or few progenitors.     

Length constraints of multi-domain proteins in metazoans

The increasing number of annotated genome sequences in public databases has made it possible to study the length distributions and domain composition of proteins at unprecedented resolution. To identify factors that influence protein length in metazoans, we performed an analysis of all domain-annotated proteins from a total of 49 animal species from Ensembl (v.56) or EnsemblMetazoa (v.3). Our results indicate that protein length constraints are not fixed as a linear function of domain count and can vary based on domain content. The presence of repeating domains was associated with relaxation of the constraints that govern protein length. Conversely, for proteins with unique domains, length constraints were generally maintained with increased domain counts. It is clear that mean (and median) protein length and domain composition vary significantly between metazoans and other kingdoms; however, the connections between function, domain content, and length are unclear. We incorporated Gene Ontology (GO) annotation to identify biological processes, cellular components, or molecular functions that favor the incorporation of multi-domain proteins. Using this approach, we identified multiple GO terms that favor the incorporation of multi-domain proteins; interestingly, several of the GO terms with elevated domain counts were not restricted to a single gene family. The findings presented here represent an important step in resolving the complex relationship between protein length, function, and domain content. The comparison of the data presented in this work to data from other kingdoms is likely to reveal additional differences in the regulation of protein length.