Design of an expression system for detecting folded protein domains and mapping macromolecular interactions by NMR.

Two protein expression vectors have been designed for the preparation of NMR samples. The vectors encode the immunoglobulin-binding domain of streptococcal protein G (GB1 domain) linked to the N-terminus of the desired proteins. This fusion strategy takes advantage of the small size, stable fold, and high bacterial expression capability of the GB1 domain to allow direct NMR spectroscopic analysis of the fusion protein by 1H-15N correlation spectroscopy. Using this system accelerates the initial assessment of protein NMR projects such that, in a matter of days, the solubility and stability of a protein can be determined. In addition, 15N-labeling of peptides and their testing for DNA binding are facilitated. Several examples are presented that demonstrate the usefulness of this technique for screening protein/DNA complexes, as well as for probing ligand-receptor interactions, using 15N-labeled GB1-peptide fusions and unlabeled target.     

Conservation of inter-residue interactions and prediction of folding rates of domain repeats

Domains are the main structural and functional units of larger proteins. They tend to be contiguous in primary structure and can fold and function independently. It has been observed that 10–20% of all encoded proteins contain duplicated domains and the average pairwise sequence identity between them is usually low. In the present study, we have analyzed the structural similarity between domain repeats of proteins with known structures available in the Protein Data Bank using structure-based inter-residue interaction measures such as the number of long-range contacts, surrounding hydrophobicity, and pairwise interaction energy. We used RADAR program for detecting the repeats in a protein sequence which were further validated using Pfam domain assignments. The sequence identity between the repeats in domains ranges from 20 to 40% and their secondary structural elements are well conserved. The number of long-range contacts, surrounding hydrophobicity calculations and pairwise interaction energy of the domain repeats clearly reveal the conservation of 3-D structure environment in the repeats of domains. The proportions of mainchain–mainchain hydrogen bonds and hydrophobic interactions are also highly conserved between the repeats. The present study has suggested that the computation of these structure-based parameters will give better clues about the tertiary environment of the repeats in domains. The folding rates of individual domains in the repeats predicted using the long-range order parameter indicate that the predicted folding rates correlate well with most of the experimentally observed folding rates for the analyzed independently folded domains.     

Finding biologically relevant protein domain interactions: conserved binding mode analysis

Proteins evolved through the shuffling of functional domains, and therefore, the same domain can be found in different proteins and species. Interactions between such conserved domains often involve specific, well-determined binding surfaces reflecting their important biological role in a cell. To find biologically relevant interactions we developed a method of systematically comparing and classifying protein domain interactions from the structural data. As a result, a set of conserved binding modes (CBMs) was created using the atomic detail of structure alignment data and the protein domain classification of the Conserved Domain Database. A conserved binding mode is inferred when different members of interacting domain families dock in the same way, such that their structural complexes superimpose well. Such domain interactions with recurring structural themes have greater significance to be biologically relevant, unlike spurious crystal packing interactions. Consequently, this study gives lower and upper bounds on the number of different types of interacting domain pairs in the structure database on the order of 1000-2000. We use CBMs to create domain interaction networks, which highlight functionally significant connections by avoiding many infrequent links between highly connected nodes. The CBMs also constitute a library of docking templates that may be used in molecular modeling to infer the characteristics of an unknown binding surface, just as conserved domains may be used to infer the structure of an unknown protein. The method's ability to sort through and classify large numbers of putative interacting domain pairs is demonstrated on the oligomeric interactions of globins.     

Prediction of VH–VL domain orientation for antibody variable domain modeling

The antigen‐binding site of antibodies forms at the interface of their two variable domains, VH and VL, making VH–VL domain orientation a factor that codetermines antibody specificity and affinity. Preserving VH–VL domain orientation in the process of antibody engineering is important in order to retain the original antibody properties, and predicting the correct VH–VL orientation has also been recognized as an important factor in antibody homology modeling. In this article, we present a fast sequence‐based predictor that predicts VH–VL domain orientation with Q² values ranging from 0.54 to 0.73 on the evaluation set. We describe VH–VL orientation in terms of the six absolute ABangle parameters that have recently been proposed as a means to separate the different degrees of freedom of VH–VL domain orientation. In order to assess the impact of adjusting VH–VL orientation according to our predictions, we use the set of antibody structures of the recently published Antibody Modeling Assessment (AMA) II study. In comparison to the original AMAII homology models, we find an improvement in the accuracy of VH–VL orientation modeling, which also translates into an improvement in the average root‐mean‐square deviation with regard to the crystal structures. Proteins 2015; 83:681–695. © 2015 Wiley Periodicals, Inc.     

Comparison of sequence and structure alignments for protein domains

Profile search methods based on protein domain alignments have proven to be useful tools in comparative sequence analysis. Domain alignments used by currently available search methods have been computed by sequence comparison. With the growth of the protein structure database, however, alignments of many domain pairs have also been computed by structure comparison. Here, we examine the extent to which information from these two sources agrees. We measure agreement with respect to identification of homologous regions in each protein, that is, with respect to the location of domain boundaries. We also measure agreement with respect to identification of homologous residue sites by comparing alignments and assessing the accuracy of the molecular models they predict. We find that domain alignments in publicly available collections based on sequence and structure comparison are largely consistent. However, the homologous regions identified by sequence comparison are often shorter than those identified by 3D structure comparison. In addition, when overall sequence similarity is low alignments from sequence comparison produce less accurate molecular models, suggesting that they less accurately identify homologous sites. These observations suggest that structure comparison results might be used to improve the overall accuracy of domain alignment collections and the performance of profile search methods based on them.     

RIP death domain structural interactions implicated in TNF-mediated proliferation and survival

Death domain (DD)-containing proteins are involved in both apoptosis and survival/proliferation signaling induced by activated death receptors. Here, a phylogenetic and structural analysis was performed to highlight differences in DD domains and their key regulatory interaction sites. The phylogenetic analysis shows that receptor DDs are more conserved than DDs in adaptors. Adaptor DDs can be subdivided into those that activate or inhibit apoptosis. Modeling of six homotypic DD interactions involved in the TNF signaling pathway implicates that the DD of RIP (Receptor interacting protein kinase 1) is capable of interacting with the DD of TRADD (TNFR1-associated death domain protein) in two different, exclusive ways: one that subsequently recruits CRADD (apoptosis/inflammation) and another that recruits NFkappaB (survival/proliferation).     

Increasing protein stability using a rational approach combining sequence homology and structural alignment: Stabilizing the WW domain

This study shows that a combination of sequence homology and structural information can be used to increase the stability of the WW domain by 2.5 kcal mol(-1) and increase the T(m) by 28 degrees C. Previous homology-based protein design efforts typically investigate positions with low sequence identity, whereas this study focuses on semi-conserved core residues and proximal residues, exploring their role(s) in mediating stabilizing interactions on the basis of structural considerations. The A20R and L30Y mutations allow increased hydrophobic interactions because of complimentary surfaces and an electrostatic interaction with a third residue adjacent to the ligand-binding hydrophobic cluster, increasing stability significantly beyond what additivity would predict for the single mutations. The D34T mutation situated in a pi-turn possibly disengages Asn31, allowing it to make up to three hydrogen bonds with the backbone in strand 1 and loop 2. The synergistic mutations A20R/L30Y in combination with the remotely located mutation D34T add together to create a hYap WW domain that is significantly more stable than any of the protein structures on which the design was based (Pin and FBP28 WW domains).     

Conservation of species interactions to achieve self-sustaining ecosystems

A desirable goal of nature management is to re-establish self-sustaining ecosystems that ensure the persistence of natural habitats and species without requiring active management. Such self-sustainability relies on functional species interactions; yet, species interactions are often overlooked in the conservation literature, and when designing species-specific management efforts. Some interactions may not be restored under general management (e.g. land protection), and may require additional specific management interventions. Interventions targeting these specific interactions fall in a gap between general and species-specific management, effectively bridging community- and population-level approaches to conservation management. We propose that managers should explicitly identify cases where active management of specific interaction partners is required to achieve population self-sustainability. In addition, they should ensure that general management interventions do not inadvertently conflict with naturally occurring interactions, potentially thwarting conservation targets. Interaction functionality may be restored by relying on native species and by identifying the spatial context in which interactions are most likely to re-establish, considering distributional range overlap of interaction partners, local variation in individual encounter rate or even spatial variation in the expected success rate (efficiency) of the focal interaction.     

Prediction of functional sites by analysis of sequence and structure conservation

We present a method for prediction of functional sites in a set of aligned protein sequences. The method selects sites which are both well conserved and clustered together in space, as inferred from the 3D structures of proteins included in the alignment. We tested the method using 86 alignments from the NCBI CDD database, where the sites of experimentally determined ligand and/or macromolecular interactions are annotated. In agreement with earlier investigations, we found that functional site predictions are most successful when overall background sequence conservation is low, such that sites under evolutionary constraint become apparent. In addition, we found that averaging of conservation values across spatially clustered sites improves predictions under certain conditions: that is, when overall conservation is relatively high and when the site in question involves a large macromolecular binding interface. Under these conditions it is better to look for clusters of conserved sites than to look for particular conserved sites.     

Unconventional interactions between water and heterocyclic nitrogens in protein structures

We report an unusual interaction in which a water molecule approaches the heterocyclic nitrogen of tryptophan and histidine along an axis that is roughly perpendicular to the aromatic plane of the side chain. The interaction is distinct from the well-known conventional aromatic hydrogen-bond, and it occurs at roughly the same frequency in protein structures. Calculations indicate that the water-indole interaction is favorable energetically, and we find several cases in which such contacts are conserved among structural orthologs. The indole-water interaction links side chains and peptide backbone in turn regions, connects the side chains in beta-sheets, and bridges secondary elements from different domains. We suggest that the water-indole interaction can be indirectly responsible for the quenching of tryptophan fluorescence that is observed in the folding of homeodomains and, possibly, many other proteins. We also observe a similar interaction between water and the imidazole nitrogens of the histidine side chain. Taken together, these observations suggest that the unconventional water-indole and water-imidazole interactions provide a small but favorable contribution to protein structures.     

Predicting protein domain interactions from coevolution of conserved regions

The knowledge of protein and domain interactions provide crucial insights into their function within a cell. Several computational methods have been proposed to detect interactions between proteins and their constitutive domains. In this work, we focus on approaches based on correlated evolution (coevolution) of sequences of interacting proteins. In this type of approach, often referred to as the mirrortree method, a high correlation of evolutionary histories of two proteins is used as an indicator to predict protein interactions. Recently, it has been observed that subtracting the underlying speciation process by separating coevolution due to common speciation divergence from that due to common function of interacting pairs greatly improves the predictive power of the mirrortree approach. In this article, we investigate possible improvements and limitations of this method. In particular, we demonstrate that the performance of the mirrortree method that can be further improved by restricting the coevolution analysis to the relatively conserved regions in the protein domain sequences (disregarding highly divergent regions). We provide a theoretical validation of our results leading to new insights into the interplay between coevolution and speciation of interacting proteins.     

Patterns of protein protein interactions in salt solutions and implications for protein crystallization

The second osmotic virial coefficients of seven proteins-ovalbumin, ribonuclease A, bovine serum albumin, alpha-lactalbumin, myoglobin, cytochrome c, and catalase-were measured in salt solutions. Comparison of the interaction trends in terms of the dimensionless second virial coefficient b(2) shows that, at low salt concentrations, protein-protein interactions can be either attractive or repulsive, possibly due to the anisotropy of the protein charge distribution. At high salt concentrations, the behavior depends on the salt: In sodium chloride, protein interactions generally show little salt dependence up to very high salt concentrations, whereas in ammonium sulfate, proteins show a sharp drop in b(2) with increasing salt concentration beyond a particular threshold. The experimental phase behavior of the proteins corroborates these observations in that precipitation always follows the drop in b(2). When the proteins crystallize, they do so at slightly lower salt concentrations than seen for precipitation. The b(2) measurements were extended to other salts for ovalbumin and catalase. The trends follow the Hofmeister series, and the effect of the salt can be interpreted as a water-mediated effect between the protein and salt molecules. The b(2) trends quantify protein-protein interactions and provide some understanding of the corresponding phase behavior. The results explain both why ammonium sulfate is among the best crystallization agents, as well as some of the difficulties that can be encountered in protein crystallization.     

Continuous and discontinuous domains: an algorithm for the automatic generation of reliable protein domain definitions.

An algorithm is presented for the fast and accurate definition of protein structural domains from coordinate data without prior knowledge of the number or type of domains. The algorithm explicitly locates domains that comprise one or two continuous segments of protein chain. Domains that include more than two segments are also located. The algorithm was applied to a nonredundant database of 230 protein structures and the results compared to domain definitions obtained from the literature, or by inspection of the coordinates on molecular graphics. For 70% of the proteins, the derived domains agree with the reference definitions, 18% show minor differences and only 12% (28 proteins) show very different definitions. Three screens were applied to identify the derived domains least likely to agree with the subjective definition set. These screens revealed a set of 173 proteins, 97% of which agree well with the subjective definitions. The algorithm represents a practical domain identification tool that can be run routinely on the entire structural database. Adjustment of parameters also allows smaller compact units to be identified in proteins.     

Long- and short-range interactions in native protein structures are consistent/minimally frustrated in sequence space

We show that long- and short-range interactions in almost all protein native structures are actually consistent with each other for coarse-grained energy scales; specifically we mean the long-range inter-residue contact energies and the short-range secondary structure energies based on peptide dihedral angles, which are potentials of mean force evaluated from residue distributions observed in protein native structures. This consistency is observed at equilibrium in sequence space rather than in conformational space. Statistical ensembles of sequences are generated by exchanging residues for each of 797 protein native structures with the Metropolis method. It is shown that adding the other category of interaction to either the short- or long-range interactions decreases the means and variances of those energies for essentially all protein native structures, indicating that both interactions consistently work by more-or-less restricting sequence spaces available to one of the interactions. In addition to this consistency, independence by these interaction classes is also indicated by the fact that there are almost no correlations between them when equilibrated using both interactions and significant but small, positive correlations at equilibrium using only one of the interactions. Evidence is provided that protein native sequences can be regarded approximately as samples from the statistical ensembles of sequences with these energy scales and that all proteins have the same effective conformational temperature. Designing protein structures and sequences to be consistent and minimally frustrated among the various interactions is a most effective way to increase protein stability and foldability.     

Three-dimensional domain swapping in homooligomeric proteins and its functional significance

In the process of oligomeric structure formation through a mechanism of three-dimensional domain swapping, one domain of a monomeric protein is replaced by the same domain from an identical monomer. The swapped domain can represent an entire tertiary globular domain or an element of secondary protein structure, such as an -helix or a -strand. Different examples of three-dimensional domain swapping are reviewed; the functional importance of this phenomenon and its role in the development of new properties by some proteins in the process of evolution are considered. The contribution of three-dimensional domain swapping to the formation of linear protein polymers and amyloids is discussed.     

On the relationship between the sequence conservation and the packing density profiles of the protein complexes

We have recently showed that the weighted contact number profiles (or the packing density profiles) of proteins are well correlated with those of the corresponding sequence conservation profiles. The results suggest that a protein structure may contain sufficient information about sequence conservation comparable to that derived from multiple homologous sequences. However, there are ambiguities concerning how to compute the packing density of the subunit of a protein complex. For the subunits of a complex, there are different ways to compute its packing density – one including the packing contributions of the other subunits and the other one excluding their contributions. Here we selected two sets of enzyme complexes. Set A contains complexes with the active sites comprising residues from multiple subunits, while set B contains those with the active sites residing on single subunits. In Set A, if the packing density profile of a subunit is computed considering the contributions of the other subunits of the complex, it will agree better with the sequence conservation profile. But in Set B the situations are reversed. The results may be due to the stronger functional and structural constraints on the evolution processes on the complexes of Set A than those of Set B to maintain the enzymatic functions of the complexes. The comparison of the packing density and the sequence conservation profiles may provide a simple yet potentially useful way to understanding the structural and evolutionary couplings between the subunits of protein complexes. Proteins 2013; 81:1192–1199. © 2013 Wiley Periodicals, Inc.     

Persistently conserved positions in structurally similar, sequence dissimilar proteins: Roles in preserving protein fold and function

Many protein pairs that share the same fold do not have any detectable sequence similarity, providing a valuable source of information for studying sequence-structure relationship. In this study, we use a stringent data set of structurally similar, sequence-dissimilar protein pairs to characterize residues that may play a role in the determination of protein structure and/or function. For each protein in the database, we identify amino-acid positions that show residue conservation within both close and distant family members. These positions are termed "persistently conserved". We then proceed to determine the "mutually" persistently conserved (MPC) positions: those structurally aligned positions in a protein pair that are persistently conserved in both pair mates. Because of their intra- and interfamily conservation, these positions are good candidates for determining protein fold and function. We find that 45% of the persistently conserved positions are mutually conserved. A significant fraction of them are located in critical positions for secondary structure determination, they are mostly buried, and many of them form spatial clusters within their protein structures. A substitution matrix based on the subset of MPC positions shows two distinct characteristics: (i) it is different from other available matrices, even those that are derived from structural alignments; (ii) its relative entropy is high, emphasizing the special residue restrictions imposed on these positions. Such a substitution matrix should be valuable for protein design experiments.     

Conservation of structural fluctuations in homologous protein kinases and its implications on functional sites

Our aim is to explore the similarities in structural fluctuations of homologous kinases. Gaussian Network Model based Normal Mode Analysis was performed on 73 active conformation structures in Ser/Thr/Tyr kinase superfamily. Categories of kinases with progressive evolutionary divergence, viz. (i) Same kinase with many crystal structures, (ii) Within‐Subfamily, (iii) Within‐Family, (iv) Within‐Group, and (v) Across‐Group, were analyzed. We identified a flexibility signature conserved in all kinases involving residues in and around the catalytic loop with consistent low‐magnitude fluctuations. However, the overall structural fluctuation profiles are conserved better in closely related kinases (Within‐Subfamily and Within‐family) than in distant ones (Within‐Group and Across‐Group). A substantial 65.4% of variation in flexibility was not accounted by variation in sequences or structures. Interestingly, we identified substructural residue‐wise fluctuation patterns characteristic of kinases of different categories. Specifically, we recognized statistically significant fluctuations unique to families of protein kinase A, cyclin‐dependent kinases, and nonreceptor tyrosine kinases. These fluctuation signatures localized to sites known to participate in protein‐protein interactions typical of these kinase families. We report for the first time that residues characterized by fluctuations unique to the group/family are involved in interactions specific to the group/family. As highlighted for Src family, local regions with differential fluctuations are proposed as attractive targets for drug design. Overall, our study underscores the importance of consideration of fluctuations, over and above sequence and structural features, in understanding the roles of sites characteristic of kinases. Proteins 2016; 84:957–978. © 2016 Wiley Periodicals, Inc.     

Interresidue interactions in protein classes

The free energy difference between folded and unfolded state is about the same for most proteins and it is not more than the energy of a few noncovalent interactions. In addition to the numerous noncovalent interactions, some proteins contain one or more disulfide bonds, which, as covalent crosslinks, significantly stabilize their tertiary structure. Correlation between the presence of disulfide bond(s), and the number noncovalent interresidue interactions of various kinds is analyzed here. The number of interactions per residue is almost the same for all protein. Also the number of long-range interactions per residue is the same in all proteins. Proteins with S(SINGLE BOND)S bond(s) (extracellular proteins) have more medium-range and fewer short-range interactions than those without S(SINGLE BOND)S bonds. However, the difference is independent of the number of these covalent crosslinks. We concluded that the different distributions of the various kinds of noncovalent interaction reflect the needs of proteins in the different environments, the extracellular and the intracellular ones, rather than the presence of the disulfide bond(s). We also pointed out that the observed differences in the distributions of short- and medium-range interactions are in good agreement with different secondary structure compositions of extracellular and intracellular proteins. Proteins 27:360–366, 1997. © 1997 Wiley-Liss, Inc.