The origins of the evolutionary signal used to predict protein-protein interactions

Background: The correlation of genetic distances between pairs of protein sequence alignments has been used to infer protein-protein interactions. It has been suggested that these correlations are based on the signal of co-evolution between interacting proteins. However, although mutations in different proteins associated with maintaining an interaction clearly occur (particularly in binding interfaces and neighbourhoods), many other factors contribute to correlated rates of sequence evolution. Proteins in the same genome are usually linked by shared evolutionary history and so it would be expected that there would be topological similarities in their phylogenetic trees, whether they are interacting or not. For this reason the underlying species tree is often corrected for. Moreover processes such as expression level, are known to effect evolutionary rates. However, it has been argued that the correlated rates of evolution used to predict protein interaction explicitly includes shared evolutionary history; here we test this hypothesis. Results: In order to identify the evolutionary mechanisms giving rise to the correlations between interaction proteins, we use phylogenetic methods to distinguish similarities in tree topologies from similarities in genetic distances. We use a range of datasets of interacting and non-interacting proteins from Saccharomyces cerevisiae. We find that the signal of correlated evolution between interacting proteins is predominantly a result of shared evolutionary rates, rather than similarities in tree topology, independent of evolutionary divergence. Conclusions: Since interacting proteins do not have tree topologies that are more similar than the control group of non-interacting proteins, it is likely that coevolution does not contribute much to, if any, of the observed correlations.

Improved Detection of Remote Homologues Using Cascade PSI-BLAST: Influence of Neighbouring Protein Families on Sequence Coverage

Background: Development of sensitive sequence search procedures for the detection of distant relationships between proteins at superfamily/fold level is still a big challenge. The intermediate sequence search approach is the most frequently employed manner of identifying remote homologues effectively. In this study, examination of serine proteases of prolyl oligopeptidase, rhomboid and subtilisin protein families were carried out using plant serine proteases as queries from two genomes including A. thaliana and O. sativa and 13 other families of unrelated folds to identify the distant homologues which could not be obtained using PSI-BLAST. Methodology/Principal Findings: We have proposed to start with multiple queries of classical serine protease members to identify remote homologues in families, using a rigorous approach like Cascade PSI-BLAST. We found that classical sequence based approaches, like PSI-BLAST, showed very low sequence coverage in identifying plant serine proteases. The algorithm was applied on enriched sequence database of homologous domains and we obtained overall average coverage of 88% at family, 77% at superfamily or fold level along with specificity of similar to 100% and Mathew's correlation coefficient of 0.91. Similar approach was also implemented on 13 other protein families representing every structural class in SCOP database. Further investigation with statistical tests, like jackknifing, helped us to better understand the influence of neighbouring protein families. Conclusions/Significance: Our study suggests that employment of multiple queries of a family for the Cascade PSI-BLAST searches is useful for predicting distant relationships effectively even at superfamily level. We have proposed a generalized strategy to cover all the distant members of a particular family using multiple query sequences. Our findings reveal that prior selection of sequences as query and the presence of neighbouring families can be important for covering the search space effectively in minimal computational time. This study also provides an understanding of the `bridging' role of related families.

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 naive Bayes classifier model to predict the interfacial residues. Our predictions are highly accurate (approximate to 85%) and specific (approximate to 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. (c) 2013 Wiley Periodicals, Inc.

The Glycosylated Rv1860 Protein of Mycobacterium tuberculosis Inhibits Dendritic Cell Mediated TH1 and TH17 Polarization of T Cells and Abrogates Protective Immunity Conferred by BCG

We previously reported interferon gamma secretion by human CD4(+) and CD8(+) T cells in response to recombinant E. coli-expressed Rv1860 protein of Mycobacterium tuberculosis (MTB) as well as protection of guinea pigs against a challenge with virulent MTB following prime-boost immunization with DNA vaccine and poxvirus expressing Rv1860. In contrast, a Statens Serum Institute Mycobacterium bovis BCG (BCG-SSI) recombinant expressing MTB Rv1860 (BCG-TB1860) showed loss of protective ability compared to the parent BCG strain expressing the control GFP protein (BCG-GFP). Since Rv1860 is a secreted mannosylated protein of MTB and BCG, we investigated the effect of BCG-TB1860 on innate immunity. Relative to BCG-GFP, BCG-TB1860 effected a significant near total reduction both in secretion of cytokines IL-2, IL-12p40, IL-12p70, TNF-alpha, IL-6 and IL-10, and up regulation of co-stimulatory molecules MHC-II, CD40, CD54, CD80 and CD86 by infected bone marrow derived dendritic cells (BMDC), while leaving secreted levels of TGF-beta unchanged. These effects were mimicked by BCG-TB1860His which carried a 6-Histidine tag at the C-terminus of Rv1860, killed sonicated preparations of BCG-TB1860 and purified H37Rv-derived Rv1860 glycoprotein added to BCG-GFP, but not by E. coli-expressed recombinant Rv1860. Most importantly, BMDC exposed to BCG-TB1860 failed to polarize allogeneic as well as syngeneic T cells to secrete IFN-gamma and IL-17 relative to BCG-GFP. Splenocytes from mice infected with BCG-SSI showed significantly less proliferation and secretion of IL-2, IFN-gamma and IL-17, but secreted higher levels of IL-10 in response to in vitro restimulation with BCG-TB1860 compared to BCG-GFP. Spleens from mice infected with BCG-TB1860 also harboured significantly fewer DC expressing MHC-II, IL-12, IL-2 and TNF-alpha compared to mice infected with BCG-GFP. Glycoproteins of MTB, through their deleterious effects on DC may thus contribute to suppress the generation of a TH1- and TH17-dominated adaptive immune response that is vital for protection against tuberculosis.

CANDO and the infinite drug discovery frontier

The Computational Analysis of Novel Drug Opportunities (CANDO) platform (http://protinfo.org/cando) uses similarity of compound-proteome interaction signatures to infer homology of compound/drug behavior. We constructed interaction signatures for 3733 human ingestible compounds covering 48,278 protein structures mapping to 2030 indications based on basic science methodologies to predict and analyze protein structure, function, and interactions developed by us and others. Our signature comparison and ranking approach yielded benchmarking accuracies of 12-25% for 1439 indications with at least two approved compounds. We prospectively validated 49/82 `high value' predictions from nine studies covering seven indications, with comparable or better activity to existing drugs, which serve as novel repurposed therapeutics. Our approach may be generalized to compounds beyond those approved by the FDA, and can also consider mutations in protein structures to enable personalization. Our platform provides a holistic multiscale modeling framework of complex atomic, molecular, and physiological systems with broader applications in medicine and engineering.

CINCINNATA in Antirrhinum majus directly modulates genes involved in cytokinin and auxin signaling

Mutations in the CINCINNATA (CIN) gene in Antirrhinum majus and its orthologs in Arabidopsis result in crinkly leaves as a result of excess growth towards the leaf margin. CIN homologs code for TCP (TEOSINTE-BRANCHED 1, CYCLOIDEA, PROLIFERATING CELL FACTOR 1 AND 2) transcription factors and are expressed in a broad zone in a growing leaf distal to the proliferation zone where they accelerate cell maturation. Although a few TCP targets are known, the functional basis of CIN-mediated leaf morphogenesis remains unclear. We compared the global transcription profiles of wild-type and the cin mutant of A. majus to identify the targets of CIN. We cloned and studied the direct targets using RNA in situ hybridization, DNA-protein interaction, chromatin immunoprecipitation and reporter gene analysis. Many of the genes involved in the auxin and cytokinin signaling pathways showed altered expression in the cin mutant. Further, we showed that CIN binds to genomic regions and directly promotes the transcription of a cytokinin receptor homolog HISTIDINE KINASE 4 (AmHK4) and an IAA3/SHY2 (INDOLE-3-ACETIC ACID INDUCIBLE 3/SHORT HYPOCOTYL 2) homolog in A. majus. Our results suggest that CIN limits excess cell proliferation and maintains the flatness of the leaf surface by directly modulating the hormone pathways involved in patterning cell proliferation and differentiation during leaf growth. 10.1111/(ISSN)1469-8137

Metallophosphoesterases: structural fidelity with functional promiscuity

Calcineurin-like metallophosphoesterases (MPEs) form a large superfamily of binuclear metal-ion-centre-containing enzymes that hydrolyse phosphomono-, phosphodi-or phosphotri-esters in a metal-dependent manner. The MPE domain is found in Mre11/SbcD DNA-repair enzymes, mammalian phosphoprotein phosphatases, acid sphingomyelinases, purple acid phosphatases, nucleotidases and bacterial cyclic nucleotide phosphodiesterases. Despite this functional diversity, MPEs show a remarkably similar structural fold and active-site architecture. In the present review, we summarize the available structural, biochemical and functional information on these proteins. We also describe how diversification and specialization of the core MPE fold in various MPEs is achieved by amino acid substitution in their active sites, metal ions and regulatory effects of accessory domains. Finally, we discuss emerging roles of these proteins as non-catalytic protein-interaction scaffolds. Thus we view the MPE superfamily as a set of proteins with a highly conserved structural core that allows embellishment to result in dramatic and niche-specific diversification of function.

Cyclic nucleotide binding and structural changes in the isolated GAF domain of Anabaena adenylyl cyclase, CyaB2

GAF domains are a large family of regulatory domains, and a subset are found associated with enzymes involved in cyclic nucleotide (cNMP) metabolism such as adenylyl cyclases and phosphodiesterases. CyaB2, an adenylyl cyclase from Anabaena, contains two GAF domains in tandem at the N-terminus and an adenylyl cyclase domain at the C-terminus. Cyclic AMP, but not cGMP, binding to the GAF domains of CyaB2 increases the activity of the cyclase domain leading to enhanced synthesis of cAMP. Here we show that the isolated GAFb domain of CyaB2 can bind both cAMP and cGMP, and enhanced specificity for cAMP is observed only when both the GAFa and the GAFb domains are present in tandem(GAFab domain). In silico docking and mutational analysis identified distinct residues important for interaction with either cAMP or cGMP in the GAFb domain. Structural changes associated with ligand binding to the GAF domains could not be detected by bioluminescence resonance energy transfer (BRET) experiments. However, amide hydrogen-deuterium exchange mass spectrometry (HDXMS) experiments provided insights into the structural basis for cAMP-induced allosteric regulation of the GAF domains, and differences in the changes induced by cAMP and cGMP binding to the GAF domain. Thus, our findings could allow the development of molecules that modulate the allosteric regulation by GAF domains present in pharmacologically relevant proteins.

The Structure of the Holo-Acyl Carrier Protein of Leishmania major Displays a Remarkably Different Phosphopantetheinyl Transferase Binding Interface

The genome of Leishmania major encodes a type II fatty acid biosynthesis pathway for which no structural or biochemical information exists. Here, for the first time, we have characterized the central player of the pathway, the acyl carrier protein (LmACP), using nuclear magnetic resonance (NMR). Structurally, the LmACP molecule is similar to other type II ACPs, comprising a four-helix bundle, enclosing a hydrophobic core. Dissimilarities in sequence, however, exist in helix II (recognition helix) of the protein. The enzymatic conversion of apo-LmACP into the holo form using type I (Escherichia coli AcpS) and type II (Sfp type) phosphopantetheinyl transferases (PPTs) is relatively slow. Mutagenesis studies underscore the importance of the residues present at the protein protein interaction interface of LmACP in modulating the activity of PPTs. Interestingly, the cognate PPT for this ACP, the L. major 4'-phosphopantetheinyl transferase (LmPPT), does not show any enzymatic activity toward it, though it readily converts other type I and type II ACPs into their holo forms. NMR chemical shift perturbation studies suggest a moderately tight complex between LmACP and its cognate PPT, suggesting inhibition. We surmise that the unique surface of LmACP might have evolved to complement its cognate enzyme (LmPPT), possibly for the purpose of regulation.

Mononuclear five-coordinate cobalt(II) complexes with N-4-coordinate pyrazole based ligand and pseudohalogens: synthesis, structures, DNA and protein binding study

A series of mononuclear five-coordinate cobalt(II) complexes, Co(dbdmp)(X)]Y, where dbdmp=N,N-diethyl-N,N-bis((3,5-dimethyl-1H-pyrazol-1-yl)methyl)ethane-1, 2-diamine, X=N-3(-)/NCO-/NCS- and Y=PF6-/BF4-/ClO4-, have been synthesized and characterized by microanalyses and spectroscopic techniques. Crystal structures of Co(N-3)(dbdmp)]PF6 (1), Co(N-3)(dbdmp)]ClO4 (3), Co(NCO)(dbdmp)]PF6 (4), Co(NCO)(dbdmp)]ClO4 (6), and Co(NCS)(dbdmp)]ClO4 (9) have been solved by single-crystal X-ray diffraction studies and showed that all the complexes have distorted trigonal bipyramidal geometry; PF6- counter anion containing complexes Co(N-3)(dbdmp)]PF6 and Co(NCO)(dbdmp)]PF6 have chiral space groups. The binding ability of synthesized complexes with CT-DNA and bovine serum albumin (BSA) has been studied by spectroscopic methods and viscosity measurements. The experimental results of absorption titration of cobalt(II) complexes with CT-DNA indicate that the complexes have ability to form adducts and they can stabilize the DNA helix. The cobalt(II) complexes exhibit good binding propensity to BSA protein.

Protein Structure and Function: Looking through the Network of Side-Chain Interactions

Network theory has become an excellent method of choice through which biological data are smoothly integrated to gain insights into complex biological problems. Understanding protein structure, folding, and function has been an important problem, which is being extensively investigated by the network approach. Since the sequence uniquely determines the structure, this review focuses on the networks of non-covalently connected amino acid side chains in proteins. Questions in structural biology are addressed within the framework of such a formalism. While general applications are mentioned in this review, challenging problems which have demanded the attention of scientific community for a long time, such as allostery and protein folding, are considered in greater detail. Our aim has been to explore these important problems through the eyes of networks. Various methods of constructing protein structure networks (PSN) are consolidated. They include the methods based on geometry, edges weighted by different schemes, and also bipartite network of protein-nucleic acid complexes. A number of network metrics that elegantly capture the general features as well as specific features related to phenomena, such as allostery and protein model validation, are described. Additionally, an integration of network theory with ensembles of equilibrium structures of a single protein or that of a large number of structures from the data bank has been presented to perceive complex phenomena from network perspective. Finally, we discuss briefly the capabilities, limitations, and the scope for further explorations of protein structure networks.

Comparative sequence analyses of genome and transcriptome reveal novel transcripts and variants in the Asian elephant Elephas maximus

The Asian elephant Elephas maximus and the African elephant Loxodonta africana that diverged 5-7 million years ago exhibit differences in their physiology, behaviour and morphology. A comparative genomics approach would be useful and necessary for evolutionary and functional genetic studies of elephants. We performed sequencing of E. maximus and map to L. africana at similar to 15X coverage. Through comparative sequence analyses, we have identified Asian elephant specific homozygous, non-synonymous single nucleotide variants (SNVs) that map to 1514 protein coding genes, many of which are involved in olfaction. We also present the first report of a high-coverage transcriptome sequence in E. maximus from peripheral blood lymphocytes. We have identified 103 novel protein coding transcripts and 66-long non-coding (lnc)RNAs. We also report the presence of 181 protein domains unique to elephants when compared to other Afrotheria species. Each of these findings can be further investigated to gain a better understanding of functional differences unique to elephant species, as well as those unique to elephantids in comparison with other mammals. This work therefore provides a valuable resource to explore the immense research potential of comparative analyses of transcriptome and genome sequences in the Asian elephant.