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.     

Ligand and pathogen specificity of the Atlantic salmon serum C-type lectin

Background

An Atlantic salmon (Salmo salar) C-type lectin (SSL) binds to mannose and related sugars as well as to the surface of Aeromonas salmonicida. To characterize this lectin as a pathogen recognition receptor in salmon, aspects of its interaction with molecules and with intact pathogens were investigated.

Methods

SSL was isolated using whole-yeast-affinity and mannan-affinity chromatography. The binding of SSL to the two major surface molecules of A. salmonicida, lipopolysaccharide (LPS) and A-layer protein was investigated by western blotting and enzyme-linked immunosorbent assays. Microbial binding specificity of SSL was examined by whole cell binding assays using a range of species. Carbohydrate ligand specificity of SSL was examined using glycan array analysis and frontal affinity chromatography.

Results

SSL showed binding to bacteria and yeast including, Pseudomonas fluorescens, A. salmonicida, A. hydrophila, Pichia pastoris, and Saccharomyces cerevisiae, but there was no detectable binding to Yersinia ruckeri. In antimicrobial assays, SSL showed no activity against Escherichia coli, Bacillus subtilis, S. cerevisiae, or A. salmonicida, but it was found to agglutinate E. coli. The major surface molecule of A. salmonicida recognized by SSL was shown to be LPS and not the A-layer protein. LPS binding was mannose-inhibitable. Glycans containing N-acetylglucosamine were shown to be predominant ligands.

Conclusion

SSL has a distinct ligand preference while allowing recognition of a wide variety of related carbohydrate structures.

General Significance

SSL is likely to function as a wide-spectrum pattern recognition protein.     

Comparative electrostatic analysis of adenylyl cyclase for isoform dependent regulation properties

The enzyme adenylyl cyclase (AC) plays a pivotal role in a variety of signal transduction pathways inside the cell, where it catalyzes the cyclization of adenosine triphosphate (ATP) into the second‐messenger cyclic adenosine monophosphate (cAMP). Among other roles, AC regulates processes involved in neural plasticity, innervation of smooth muscles of the heart and the endocrine system of the pancreas. The functional diversity of AC is manifested in its different isoforms, each having a specific regulation pattern. There is an increasing amount of data available concerning the regulatory properties of AC isoforms, however little is known about the interactions on a structural level. Here, we conducted a comparative electrostatic analysis of the catalytic domains of all nine transmembrane AC isoforms with the aim of detecting, verifying and predicting the binding sites of molecular regulators on AC. The results provide support for the positioning of the binding site of the inhibitory protein G_iα at a pseudo‐symmetric position to the stimulatory G_sα binding site. They also provide a structural interpretation of the Gβγ interaction with ACs 2, 4, and 7 and suggest a new binding site for RGS2. Comparison of the small molecule binding sites on AC shows that overall they have high electrostatic similarity, but regions of electrostatic differences are identified. These could provide a basis for the development of novel compounds with isoform‐specific modulatory effects on AC. Proteins 2016; 84:1844–1858. © 2016 Wiley Periodicals, Inc.     

C-type lectin from red swamp crayfish Procambarus clarkii participates in cellular immune response

Lectins are potential immune recognition proteins. In this study, a novel C-type lectin (Pc-Lec1) is reported in freshwater crayfish Procambarus clarkii. Pc-Lec1 encodes a protein of 163 amino acids with a putative signal peptide and a single carbohydrate recognition domain. It was constitutively expressed in various tissues of a normal crayfish, especially in the hepatopancreas and gills. Expressions of Pc-Lec1 were up-regulated in the hepatopancreas and gills of crayfish challenged with Vibrio anguillarum, Staphylococcus aureus, or the white spot syndrome virus. Recombinant mature Pc-Lec1 bound bacteria and polysaccharides (peptidoglycan, lipoteichoic acid, and lipopolysaccharide) but did not agglutinate bacteria. Pc-Lec1 enhanced hemocyte encapsulation of the sepharose beads in vitro, and the blocking of beads by a polyclonal antibody inhibited encapsulation. Pc-Lec1 promoted clearance of V. anguillarum in vivo. These results suggest that Pc-Lec1 is a pattern recognition receptor and participates in cellular immune response. Pc-Lec1 performs its function as an opsonin by enhancing the encapsulation or clearance of pathogenic bacteria.     

The PYRIN domain: a member of the death domain-fold superfamily

PYRIN domains were identified recently as putative protein-protein interaction domains at the N-termini of several proteins thought to function in apoptotic and inflammatory signaling pathways. The approximately 95 residue PYRIN domains have no statistically significant sequence homology to proteins with known three-dimensional structure. Using secondary structure prediction and potential-based fold recognition methods, however, the PYRIN domain is predicted to be a member of the six-helix bundle death domain-fold superfamily that includes death domains (DDs), death effector domains (DEDs), and caspase recruitment domains (CARDs). Members of the death domain-fold superfamily are well established mediators of protein-protein interactions found in many proteins involved in apoptosis and inflammation, indicating further that the PYRIN domains serve a similar function. An homology model of the PYRIN domain of CARD7/DEFCAP/NAC/NALP1, a member of the Apaf-1/Ced-4 family of proteins, was constructed using the three-dimensional structures of the FADD and p75 neurotrophin receptor DDs, and of the Apaf-1 and caspase-9 CARDs, as templates. Validation of the model using a variety of computational techniques indicates that the fold prediction is consistent with the sequence. Comparison of a circular dichroism spectrum of the PYRIN domain of CARD7/DEFCAP/NAC/NALP1 with spectra of several proteins known to adopt the death domain-fold provides experimental support for the structure prediction.     

Evolutionary history and stress regulation of the lectin superfamily in higher plants

Background  

Lectins are a class of carbohydrate-binding proteins. They play roles in various biological processes. However, little is known about their evolutionary history and their functions in plant stress regulation. The availability of full genome sequences from various plant species makes it possible to perform a whole-genome exploration for further understanding their biological functions.     

Structure,function and evolution of the hemerythrin‐like domain superfamily

Hemerythrin‐like proteins have generally been studied for their ability to reversibly bind oxygen through their binuclear nonheme iron centers. However, in recent years, it has become increasingly evident that some members of the hemerythrin‐like superfamily also participate in many other biological processes. For instance, the binuclear nonheme iron site of YtfE, a hemerythrin‐like protein involved in the repair of iron centers in Escherichia coli, catalyzes the reduction of nitric oxide to nitrous oxide, and the human F‐box/LRR‐repeat protein 5, which contains a hemerythrin‐like domain, is involved in intracellular iron homeostasis. Furthermore, structural data on hemerythrin‐like domains from two proteins of unknown function, PF0695 from Pyrococcus furiosus and NMB1532 from Neisseria meningitidis, show that the cation‐binding sites, typical of hemerythrin, can be absent or be occupied by metal ions other than iron. To systematically investigate this functional and structural diversity of the hemerythrin‐like superfamily, we have collected hemerythrin‐like sequences from a database comprising fully sequenced proteomes and generated a cluster map based on their all‐against‐all pairwise sequence similarity. Our results show that the hemerythrin‐like superfamily comprises a large number of protein families which can be classified into three broad groups on the basis of their cation‐coordinating residues: (a) signal‐transduction and oxygen‐carrier hemerythrins (H‐HxxxE‐HxxxH‐HxxxxD); (b) hemerythrin‐like (H‐HxxxE‐H‐HxxxE); and, (c) metazoan F‐box proteins (H‐HExxE‐H‐HxxxE). Interestingly, all but two hemerythrin‐like families exhibit internal sequence and structural symmetry, suggesting that a duplication event may have led to the origin of the hemerythrin domain.     

Origins of the specificity of inhibitor P218 toward wild-type and mutant PfDHFR: a molecular dynamics analysis

Molecular dynamics simulations were performed to evaluate the origin of the antimalarial effect of the lead compound P218. The simulations of the ligand in the cavities of wild-type, mutant Plasmodium falciparum Dihydrofolate Reductase (PfDHFR) and the human DHFR revealed the differences in the atomic-level interactions and also provided explanation for the specificity of this ligand toward PfDHFR. The binding free energy estimation using Molecular Mechanics Poisson-Boltzmann Surface Area method revealed that P218 has higher binding affinity (~ ?30 to ?35 kcal/mol) toward PfDHFR (both in wild-type and mutant forms) than human DHFR (~ ?22 kcal/mol), corroborating the experimental observations. Intermolecular hydrogen bonding analysis of the trajectories showed that P218 formed two stable hydrogen bonds with human DHFR (Ile7 and Glu30), wild-type and double-mutant PfDHFR’s (Asp54 and Arg122), while it formed three stable hydrogen bonds with quadruple-mutant PfDHFR (Asp54, Arg59, and Arg122). Additionally, P218 binding in PfDHFR is stabilized by hydrogen bonds with residues Ile14 and Ile164. It was found that mutant residues do not reduce the binding affinity of P218 to PfDHFR, in contrast, Cys59Arg mutation strongly favors inhibitor binding to quadruple-mutant PfDHFR. The atomistic-level details explored in this work will be highly useful for the design of non-resistant novel PfDHFR inhibitors as antimalarial agents.     

Integrating phylogenetic and ecological distances reveals new insights into parasite host specificity

The range of hosts a pathogen infects (host specificity) is a key element of disease risk that may be influenced by both shared phylogenetic history and shared ecological attributes of prospective hosts. Phylospecificity indices quantify host specificity in terms of host relatedness, but can fail to capture ecological attributes that increase susceptibility. For instance, similarity in habitat niche may expose phylogenetically unrelated host species to similar pathogen assemblages. Using a recently proposed method that integrates multiple distances, we assess the relative contributions of host phylogenetic and functional distances to pathogen host specificity (functional–phylogenetic host specificity). We apply this index to a data set of avian malaria parasite (Plasmodium and Haemoproteus spp.) infections from Melanesian birds to show that multihost parasites generally use hosts that are closely related, not hosts with similar habitat niches. We also show that host community phylogenetic ß‐diversity (Pßd) predicts parasite Pßd and that individual host species carry phylogenetically clustered Haemoproteus parasite assemblages. Our findings were robust to phylogenetic uncertainty, and suggest that phylogenetic ancestry of both hosts and parasites plays important roles in driving avian malaria host specificity and community assembly. However, restricting host specificity analyses to either recent or historical timescales identified notable exceptions, including a ‘habitat specialist’ parasite that infects a diversity of unrelated host species with similar habitat niches. This work highlights that integrating ecological and phylogenetic distances provides a powerful approach to better understand drivers of pathogen host specificity and community assembly.     

Computational models explain the oligosaccharide specificity of cyanovirin-N

The prokaryotic lectin cyanovirin-N (CV-N) is a potent inhibitor of HIV envelope-mediated cell entry, and thus is a leading candidate among a new class of potential anti-HIV microbicides. The activity of CV-N is a result of interactions with the D1 arm of high-mannose oligosaccharides on the viral glycoprotein gp120. Here, we present computationally refined models of CV-N recognition of the di- and trisaccharides that represent the terminal three sugars of the D1 arm by each CV-N binding site. These models complement existing structural data, both from NMR spectroscopy and X-ray crystallography. When used with a molecular dynamics/continuum electrostatic (MD/PBSA) approach to compute binding free energies, these models explain the relative affinity of each site for the two saccharides. This work presents the first validation of the application of continuum electrostatic models to carbohydrate-protein association. Taken as a whole, the results both provide models of CV-N sugar recognition and demonstrate the utility of these computational methods for the study of carbohydrate-binding proteins.     

Multivariate sequence analysis reveals additional function impacting residues in the SDR superfamily

The “extended” type of short chain dehydrogenases/reductases (SDR), share a remarkable similarity in their tertiary structures inspite of being highly divergent in their functions and sequences. We have carried out principal component analysis (PCA) on structurally equivalent residue positions of 10 SDR families using information theoretic measures like Jensen–Shannon divergence and average shannon entropy as variables. The results classify residue positions in the SDR fold into six groups, one of which is characterized by low Shannon entropies but high Jensen–Shannon divergence against the reference family SDR1E, suggesting that these positions are responsible for the specific functional identities of individual SDR families, distinguishing them from the reference family SDR1E. Site directed mutagenesis of three residues from this group in the enzyme UDP‐Galactose 4‐epimerase belonging to SDR1E shows that the mutants promote the formation of NADH containing abortive complexes. Finally, molecular dynamics simulations have been used to suggest a mechanism by which the mutants interfere with the re‐oxidation of NADH leading to the formation of abortive complexes. Proteins 2014; 82:2842–2856. © 2014 Wiley Periodicals, Inc.     

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.     

Conformation of the Ras-binding domain of Raf studied by molecular dynamics and free energy simulations

Recognition of Ras by its downstream target Raf is mediated by a Ras-recognition region in the Ras-binding domain (RBD) of Raf. Residues 78–89 in this region occupy two different conformations in the ensemble of NMR solution structures of the RBD: a fully α-helical one, and one where 87–90 form a type IV β-turn. Molecular dynamics simulations of the RBD in solution were performed to explore the stability of these and other possible conformations of both the wild-type RBD and the R89K mutant, which does not bind Ras. The simulations sample a fully helical conformation for residues 78–89 similar to the NMR helical structures, a conformation where 85–89 form a 3₁₀-helical turn, and a conformation where 87–90 form a type I |iB-turn, whose free energies are all within 0.3 kcal/mol of each other. NOE patterns and H^α chemical shifts from the simulations are in reasonable agreement with experiment. The NMR turn structure is calculated to be 3 kcal/mol higher than the three above conformations. In a simulation with the same implicit solvent model used in the NMR structure generation, the turn conformation relaxes into the fully helical conformation, illustrating possible structural artifacts introduced by the implicit solvent model. With the Raf R89K mutant, simulations sample a fully helical and a turn conformation, the turn being 0.9 kcal/mol more stable. Thus, the mutation affects the population of RBD conformations, and this is expected to affect Ras binding. For example, if the fully helical conformation of residues 78–89 is required for binding, its free energy increase in R89K will increase the binding free energy by about 0.6 kcal/mol. Proteins 31:186–200, 1998. © 1998 Wiley-Liss, Inc.     

Kinetic consequences of native state optimization of surface-exposed electrostatic interactions in the Fyn SH3 domain

Optimization of surface exposed charge-charge interactions in the native state has emerged as an effective means to enhance protein stability; but the effect of electrostatic interactions on the kinetics of protein folding is not well understood. To investigate the kinetic consequences of surface charge optimization, we characterized the folding kinetics of a Fyn SH3 domain variant containing five amino acid substitutions that was computationally designed to optimize surface charge-charge interactions. Our results demonstrate that this optimized Fyn SH3 domain is stabilized primarily through an eight-fold acceleration in the folding rate. Analyses of the constituent single amino acid substitutions indicate that the effects of optimization of charge-charge interactions on folding rate are additive. This is in contrast to the trend seen in folded state stability, and suggests that electrostatic interactions are less specific in the transition state compared to the folded state. Simulations of the transition state using a coarse-grained chain model show that native electrostatic contacts are weakly formed, thereby making the transition state conducive to nonspecific, or even nonnative, electrostatic interactions. Because folding from the unfolded state to the folding transition state for small proteins is accompanied by an increase in charge density, nonspecific electrostatic interactions, that is, generic charge density effects can have a significant contribution to the kinetics of protein folding. Thus, the interpretation of the effects of amino acid substitutions at surface charged positions may be complicated and consideration of only native-state interactions may fail to provide an adequate picture.     

Homology modeling of the Abl-SH3 domain

A tertiary structure model of the Abl-SH3 domain is predicted by using homology modeling techniques coupled to molecular dynamics simulations. Two template proteins were used, Fyn-SH3 and Spc-SH3. The refined model was extensively checked for errors using criteria based on stereochemistry, packing, solvation free-energy, accessible surface areas, and contact analyses. The different checking methods do not totally agree, as each one evaluates a different characteristic of protein structures. Several zones of the protein are more susceptible to incorporating errors. These include residues 13, 15, 35, 39, 45, 46, 50, and 60. An interesting finding is that the measurement of the Cα chirality correlated well with the rest of the criteria, suggesting that this parameter might be a good indicator of correct local conformation. Deviations of more than 4 degrees may be indicative of poor local structure. © 1994 Wiley-Liss, Inc.     

Beyond PKA: Evolutionary and structural insights that define a docking and dimerization domain superfamily

Protein-interaction domains can create unique macromolecular complexes that drive evolutionary innovation. By combining bioinformatic and phylogenetic analyses with structural approaches, we have discovered that the docking and dimerization (D/D) domain of the PKA regulatory subunit is an ancient and conserved protein fold. An archetypal function of this module is to interact with A-kinase-anchoring proteins (AKAPs) that facilitate compartmentalization of this key cell-signaling enzyme. Homology searching reveals that D/D domain proteins comprise a superfamily with 18 members that function in a variety of molecular and cellular contexts. Further in silico analyses indicate that D/D domains segregate into subgroups on the basis of their similarity to type I or type II PKA regulatory subunits. The sperm autoantigenic protein 17 (SPA17) is a prototype of the type II or R2D2 subgroup that is conserved across metazoan phyla. We determined the crystal structure of an extended D/D domain from SPA17 (amino acids 1–75) at 1.72 Å resolution. This revealed a four-helix bundle-like configuration featuring terminal β-strands that can mediate higher order oligomerization. In solution, SPA17 forms both homodimers and tetramers and displays a weak affinity for AKAP18. Quantitative approaches reveal that AKAP18 binding occurs at nanomolar affinity when SPA17 heterodimerizes with the ropporin-1-like D/D protein. These findings expand the role of the D/D fold as a versatile protein-interaction element that maintains the integrity of macromolecular architectures within organelles such as motile cilia.     

Phosphorylation-induced transient intrinsic structure in the kinase-inducible domain of CREB facilitates its recognition by the KIX domain of CBP

Phosphorylation at Ser-133 of the kinase inducible domain of CREB (KID) triggers its binding to the KIX domain of CBP via a concomitant coil-to-helix transition. The exact role of this key event is still puzzling: it does not switch between disordered and ordered states, nor its direct interactions fully account for selectivity. Hence, we reasoned that phosphorylation may shift the conformational preferences of KID towards a binding-competent state. To this end we investigated the intrinsic structural properties of the unbound KID in phosphorylated and unphosphorylated forms by simulated annealing and molecular dynamics simulations. Although helical populations show subtle differences, phosphorylation reduces the flexibility of the turn segment connecting the two helices in the complexed structure and induces a transient structural element that corresponds to its bound conformation. It is stabilized by the pSer-133-Arg-131 interaction, which is absent from the unphosphorylated KID. Diminishing this coupling decreases the 3.1 kcal/mol contribution of pSer-133 to the binding free energy (DeltaGbind) of the phosphorylated KID to KIX by 1.1 kcal/mol, as computed in reference to Ser-133. In a binding competent form of the S133E KID mutant, the contribution of Glu-133 to DeltaGbind is by 1.5 kcal/mol smaller than that of pSer, suggesting that altered structural properties due to pSer --> Glu replacement impair the binding affinity. Thus, we propose that phoshorylation contributes to selectivity not merely by the direct interactions of the phosphate group with KIX, but also by promoting the formation of a transient structural element in the highly conserved turn segment.     

Homology modelling of the core domain of the endogenous lectin comitin: structural basis for its mannose-binding specificity

The N-terminal core domain of comitin from the slime mold Dictyostelium discoideum has been modelled from the X-ray coordinates of the monocot mannose-binding lectin from snowdrop (Galanthus nivalis). Docking experiments performed on the three-dimensional model showed that two of the three mannose-binding sites of the comitin monomer are functional. They are located at both ends of the comitin dimer whereas the actin-interacting region occurs in the central hinge region where both monomers are non covalently associated. This distribution is fully consistent with the bifunctional character of comitin which is believed to link the Golgi vesicles exhibiting mannosylated membrane glycans to the actin cytoskeleton in the cell.     

Crystal structure analysis reveals a spring‐loaded latch as molecular mechanism for GDF‐5–type I receptor specificity

Dysregulation of growth and differentiation factor 5 (GDF‐5) signalling, a member of the TGF‐β superfamily, is strongly linked to skeletal malformation. GDF‐5‐mediated signal transduction involves both BMP type I receptors, BMPR‐IA and BMPR‐IB. However, mutations in either GDF‐5 or BMPR‐IB lead to similar phenotypes, indicating that in chondrogenesis GDF‐5 signalling seems to be exclusively mediated through BMPR‐IB. Here, we present structural insights into the GDF‐5:BMPR‐IB complex revealing how binding specificity for BMPR‐IB is generated on a molecular level. In BMPR‐IB, a loop within the ligand‐binding epitope functions similar to a latch allowing high‐affinity binding of GDF‐5. In BMPR‐IA, this latch is in a closed conformation leading to steric repulsion. The new structural data now provide also a molecular basis of how phenotypically relevant missense mutations in GDF‐5 might impair receptor binding and activation.