1
|
Gao Q, Zhang C, Li M, Yu T. Protein-Protein Interaction Prediction Model Based on ProtBert-BiGRU-Attention. J Comput Biol 2024. [PMID: 39069885 DOI: 10.1089/cmb.2023.0297] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024] Open
Abstract
The physiological activities within cells are mainly regulated through protein-protein interactions (PPI). Therefore, studying protein interactions has become an essential part of researching protein function and mechanisms. Traditional biological experiments required for PPI prediction are expensive and time consuming. For this reason, many methods based on predicting PPI from protein sequences have been proposed in recent years. However, existing computational methods usually require the combination of evolutionary feature information of proteins to predict PPI docking situations. Because different relevant features of selected proteins are chosen, there may be differences in the predicted results for PPI. This article proposes a PPI prediction method based on the pretrained protein sequence model ProtBert, combined with the Bidirectional Gated Recurrent Unit (BiGRU) and attention mechanism. Only using protein sequence information and leveraging ProtBert's powerful ability to capture amino acid feature information, BiGRU is used for further feature extraction of the amino acid vectors output by ProtBert. The attention mechanism is then applied to enhance the focus on different amino acid features and improve the expression ability of protein sequence features, ultimately obtaining binary classification results for protein interactions. Experimental results show that our proposed ProtBert-BiGRU-Attention model has good predictive performance for PPI. Through relevant comparative experiments, it has been proven that our model performs well in protein binary prediction. Furthermore, through the ablation experiment of the model, different deep learning modules' contributions to the prediction have been demonstrated.
Collapse
Affiliation(s)
- Qian Gao
- College of Computer and Control Engineering, Qiqihar University, Qiqihar, China
| | - Chi Zhang
- College of Computer and Control Engineering, Qiqihar University, Qiqihar, China
| | - Ming Li
- College of Computer and Control Engineering, Qiqihar University, Qiqihar, China
| | - Tianfei Yu
- College of Life Science and Agriculture Forestry, Qiqihar University, Qiqihar, China
| |
Collapse
|
2
|
Little J, Chikina M, Clark NL. Evolutionary rate covariation is a reliable predictor of co-functional interactions but not necessarily physical interactions. eLife 2024; 12:RP93333. [PMID: 38415754 PMCID: PMC10942632 DOI: 10.7554/elife.93333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024] Open
Abstract
Co-functional proteins tend to have rates of evolution that covary over time. This correlation between evolutionary rates can be measured over the branches of a phylogenetic tree through methods such as evolutionary rate covariation (ERC), and then used to construct gene networks by the identification of proteins with functional interactions. The cause of this correlation has been hypothesized to result from both compensatory coevolution at physical interfaces and nonphysical forces such as shared changes in selective pressure. This study explores whether coevolution due to compensatory mutations has a measurable effect on the ERC signal. We examined the difference in ERC signal between physically interacting protein domains within complexes compared to domains of the same proteins that do not physically interact. We found no generalizable relationship between physical interaction and high ERC, although a few complexes ranked physical interactions higher than nonphysical interactions. Therefore, we conclude that coevolution due to physical interaction is weak, but present in the signal captured by ERC, and we hypothesize that the stronger signal instead comes from selective pressures on the protein as a whole and maintenance of the general function.
Collapse
Affiliation(s)
- Jordan Little
- Department of Human Genetics, University of UtahSalt Lake CityUnited States
| | - Maria Chikina
- Department of Computational Biology, University of PittsburghPittsburghUnited States
| | - Nathan L Clark
- Department of Human Genetics, University of UtahSalt Lake CityUnited States
- Department of Biological Sciences, University of PittsburghPittsburghUnited States
| |
Collapse
|
3
|
Kogay R, Zhaxybayeva O. Co-evolution of gene transfer agents and their alphaproteobacterial hosts. J Bacteriol 2024; 206:e0039823. [PMID: 38240570 PMCID: PMC10883770 DOI: 10.1128/jb.00398-23] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Accepted: 12/19/2023] [Indexed: 02/23/2024] Open
Abstract
Gene transfer agents (GTAs) are enigmatic elements that resemble small viruses and are known to be produced during nutritional stress by some bacteria and archaea. The production of GTAs is regulated by quorum sensing, under which a small fraction of the population acts as GTA producers, while the rest becomes GTA recipients. In contrast to canonical viruses, GTAs cannot propagate themselves because they package pieces of the producing cell's genome. In alphaproteobacteria, GTAs are mostly vertically inherited and reside in their hosts' genomes for hundreds of millions of years. While GTAs' ability to transfer genetic material within a population and their long-term preservation suggest an increased fitness of GTA-producing microbes, the associated benefits and type of selection that maintains GTAs are poorly understood. By comparing rates of evolutionary change in GTA genes to the rates in gene families abundantly present across 293 alphaproteobacterial genomes, we detected 59 gene families that likely co-evolve with GTA genes. These gene families are predominantly involved in stress response, DNA repair, and biofilm formation. We hypothesize that biofilm formation enables the physical proximity of GTA-producing cells, limiting GTA-derived benefits only to a group of closely related cells. We further conjecture that the population structure of biofilm-forming sub-populations ensures that the trait of GTA production is maintained despite the inevitable rise of "cheating" genotypes. Because release of GTA particles kills the producing cell, maintenance of GTAs is an exciting example of social evolution in a microbial population.IMPORTANCEGene transfer agents (GTAs) are viruses domesticated by some archaea and bacteria as vehicles for carrying pieces of the host genome. Produced under certain environmental conditions, GTA particles can deliver DNA to neighboring, closely related cells. The function of GTAs remains uncertain. While making GTAs is suicidal for a cell, GTA-encoding genes are widespread in genomes of alphaproteobacteria. Such GTA persistence implies functional benefits but raises questions about how selection maintains this lethal trait. By showing that GTA genes co-evolve with genes involved in stress response, DNA repair, and biofilm formation, we provide support for the hypothesis that GTAs facilitate DNA exchange during the stress conditions and present a model for how GTAs persist in biofilm-forming bacterial populations despite being lethal.
Collapse
Affiliation(s)
- Roman Kogay
- Department of Biological Sciences, Dartmouth College, Hanover, New Hampshire, USA
| | - Olga Zhaxybayeva
- Department of Biological Sciences, Dartmouth College, Hanover, New Hampshire, USA
- Department of Computer Science, Dartmouth College, Hanover, New Hampshire, USA
| |
Collapse
|
4
|
Chen J, Gu Z, Lai L, Pei J. In silico protein function prediction: the rise of machine learning-based approaches. MEDICAL REVIEW (2021) 2023; 3:487-510. [PMID: 38282798 PMCID: PMC10808870 DOI: 10.1515/mr-2023-0038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Accepted: 10/11/2023] [Indexed: 01/30/2024]
Abstract
Proteins function as integral actors in essential life processes, rendering the realm of protein research a fundamental domain that possesses the potential to propel advancements in pharmaceuticals and disease investigation. Within the context of protein research, an imperious demand arises to uncover protein functionalities and untangle intricate mechanistic underpinnings. Due to the exorbitant costs and limited throughput inherent in experimental investigations, computational models offer a promising alternative to accelerate protein function annotation. In recent years, protein pre-training models have exhibited noteworthy advancement across multiple prediction tasks. This advancement highlights a notable prospect for effectively tackling the intricate downstream task associated with protein function prediction. In this review, we elucidate the historical evolution and research paradigms of computational methods for predicting protein function. Subsequently, we summarize the progress in protein and molecule representation as well as feature extraction techniques. Furthermore, we assess the performance of machine learning-based algorithms across various objectives in protein function prediction, thereby offering a comprehensive perspective on the progress within this field.
Collapse
Affiliation(s)
- Jiaxiao Chen
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
| | - Zhonghui Gu
- Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
| | - Luhua Lai
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
- Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
- BNLMS, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
- Research Unit of Drug Design Method, Chinese Academy of Medical Sciences (2021RU014), Beijing, China
| | - Jianfeng Pei
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
- Research Unit of Drug Design Method, Chinese Academy of Medical Sciences (2021RU014), Beijing, China
| |
Collapse
|
5
|
Bhalla D, van Noort V. Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes. J Mol Evol 2023; 91:628-646. [PMID: 37392220 DOI: 10.1007/s00239-023-10124-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2023] [Accepted: 06/13/2023] [Indexed: 07/03/2023]
Abstract
The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.
Collapse
Affiliation(s)
- Diksha Bhalla
- Centre of Microbial and Plant Genetics, Faculty of Bioscience Engineering, KU Leuven, Leuven, Belgium.
| | - Vera van Noort
- Centre of Microbial and Plant Genetics, Faculty of Bioscience Engineering, KU Leuven, Leuven, Belgium
- Institute of Biology Leiden, Leiden University, Leiden, The Netherlands
| |
Collapse
|
6
|
Vukičević D, Matijević D. The Connection of the Generalized Robinson-Foulds Metric with Partial Wiener Indices. Acta Biotheor 2023; 71:5. [PMID: 36695929 DOI: 10.1007/s10441-023-09457-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Accepted: 01/11/2023] [Indexed: 01/26/2023]
Abstract
In this work we propose the partial Wiener index as one possible measure of branching in phylogenetic evolutionary trees. We establish the connection between the generalized Robinson-Foulds (RF) metric for measuring the similarity of phylogenetic trees and partial Wiener indices by expressing the number of conflicting pairs of edges in the generalized RF metric in terms of partial Wiener indices. To do so we compute the minimum and maximum value of the partial Wiener index [Formula: see text], where [Formula: see text] is a binary rooted tree with root [Formula: see text] and [Formula: see text] leaves. Moreover, under the Yule probabilistic model, we show how to compute the expected value of [Formula: see text]. As a direct consequence, we give exact formulas for the upper bound and the expected number of conflicting pairs. By doing so we provide a better theoretical understanding of the computational complexity of the generalized RF metric.
Collapse
Affiliation(s)
- Damir Vukičević
- Department of Mathematics, Faculty of Science, University of Split, Ruđera Boškovića 33, 21000, Split, Croatia.
| | - Domagoj Matijević
- Department of Mathematics, University of Osijek, Trg Lj. Gaja 6, 31000, Osijek, Croatia
| |
Collapse
|
7
|
Gerardos A, Dietler N, Bitbol AF. Correlations from structure and phylogeny combine constructively in the inference of protein partners from sequences. PLoS Comput Biol 2022; 18:e1010147. [PMID: 35576238 PMCID: PMC9135348 DOI: 10.1371/journal.pcbi.1010147] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2021] [Revised: 05/26/2022] [Accepted: 04/27/2022] [Indexed: 11/19/2022] Open
Abstract
Inferring protein-protein interactions from sequences is an important task in computational biology. Recent methods based on Direct Coupling Analysis (DCA) or Mutual Information (MI) allow to find interaction partners among paralogs of two protein families. Does successful inference mainly rely on correlations from structural contacts or from phylogeny, or both? Do these two types of signal combine constructively or hinder each other? To address these questions, we generate and analyze synthetic data produced using a minimal model that allows us to control the amounts of structural constraints and phylogeny. We show that correlations from these two sources combine constructively to increase the performance of partner inference by DCA or MI. Furthermore, signal from phylogeny can rescue partner inference when signal from contacts becomes less informative, including in the realistic case where inter-protein contacts are restricted to a small subset of sites. We also demonstrate that DCA-inferred couplings between non-contact pairs of sites improve partner inference in the presence of strong phylogeny, while deteriorating it otherwise. Moreover, restricting to non-contact pairs of sites preserves inference performance in the presence of strong phylogeny. In a natural data set, as well as in realistic synthetic data based on it, we find that non-contact pairs of sites contribute positively to partner inference performance, and that restricting to them preserves performance, evidencing an important role of phylogeny.
Collapse
Affiliation(s)
- Andonis Gerardos
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Nicola Dietler
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Anne-Florence Bitbol
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
- * E-mail:
| |
Collapse
|
8
|
Steenwyk JL, Phillips MA, Yang F, Date SS, Graham TR, Berman J, Hittinger CT, Rokas A. An orthologous gene coevolution network provides insight into eukaryotic cellular and genomic structure and function. SCIENCE ADVANCES 2022; 8:eabn0105. [PMID: 35507651 PMCID: PMC9067921 DOI: 10.1126/sciadv.abn0105] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Accepted: 03/16/2022] [Indexed: 06/14/2023]
Abstract
The evolutionary rates of functionally related genes often covary. We present a gene coevolution network inferred from examining nearly 3 million orthologous gene pairs from 332 budding yeast species spanning ~400 million years of evolution. Network modules provide insight into cellular and genomic structure and function. Examination of the phenotypic impact of network perturbation using deletion mutant data from the baker's yeast Saccharomyces cerevisiae, which were obtained from previously published studies, suggests that fitness in diverse environments is affected by orthologous gene neighborhood and connectivity. Mapping the network onto the chromosomes of S. cerevisiae and Candida albicans revealed that coevolving orthologous genes are not physically clustered in either species; rather, they are often located on different chromosomes or far apart on the same chromosome. The coevolution network captures the hierarchy of cellular structure and function, provides a roadmap for genotype-to-phenotype discovery, and portrays the genome as a linked ensemble of genes.
Collapse
Affiliation(s)
- Jacob L. Steenwyk
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
| | - Megan A. Phillips
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
| | - Feng Yang
- Shmunis School of Biomedical and Cancer Research, Tel Aviv University, Ramat Aviv, Israel
- Department of Pharmacology, Shanghai Tenth People’s Hospital, Tongji University School of Medicine, Shanghai, China
| | - Swapneeta S. Date
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
| | - Todd R. Graham
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
| | - Judith Berman
- Shmunis School of Biomedical and Cancer Research, Tel Aviv University, Ramat Aviv, Israel
| | - Chris Todd Hittinger
- Laboratory of Genetics, DOE Great Lakes Bioenergy Research Center, Wisconsin Energy Institute, Center for Genomic Science Innovation, J.F. Crow Institute for the Study of Evolution, University of Wisconsin-Madison, Madison, WI, USA
| | - Antonis Rokas
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
| |
Collapse
|
9
|
Abstract
Since the large-scale experimental characterization of protein–protein interactions (PPIs) is not possible for all species, several computational PPI prediction methods have been developed that harness existing data from other species. While PPI network prediction has been extensively used in eukaryotes, microbial network inference has lagged behind. However, bacterial interactomes can be built using the same principles and techniques; in fact, several methods are better suited to bacterial genomes. These predicted networks allow systems-level analyses in species that lack experimental interaction data. This review describes the current network inference and analysis techniques and summarizes the use of computationally-predicted microbial interactomes to date.
Collapse
|
10
|
Palenchar PM. The Influence of Codon Usage, Protein Abundance, and Protein Stability on Protein Evolution Vary by Evolutionary Distance and the Type of Protein. Protein J 2022; 41:216-229. [PMID: 35147896 DOI: 10.1007/s10930-022-10045-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/01/2022] [Indexed: 12/01/2022]
Abstract
In general, the evolutionary rate of proteins is not primarily related to protein and amino acid functions, and factors such as protein abundance, codon usage, and the protein's TM are more important. To better understand the factors that affect protein evolution, E. coli MG1655 orthologs were compared to those in closely related bacteria and to more distantly related prokaryotes, eukaryotes, and archaea. Also, the evolution of different types of proteins was studied. The analyses indicate that the amino acid conservation of enzymes that do not use macromolecules (e.g. DNA, RNA, and proteins) as substrates and that carry out metabolic processes involving small molecules (i.e. small molecule enzymes) is different than other enzymes. For example, the small molecule enzymes have a lower percent identity than other enzymes when sequences from closely related bacteria are compared. Analyses indicate the lower percent identity is not a result of the amino acid or codon usage of the small molecule enzymes. The small molecule enzymes also don't have a significantly lower protein abundance indicating that is also not likely an important factor driving differences in amino acid conservation. Analyses indicate different methods to measure the TM of proteins have different relationships between amino acid conservation over different evolutionary distances. In totality, the results demonstrate that the relationship between the factors thought to affect protein evolution (protein abundance, codon usage, and proteins TMs) and protein evolution are complex and depend on the factor, the organisms, and the type of proteins being analyzed.
Collapse
Affiliation(s)
- Peter M Palenchar
- Department of Chemistry, Villanova University, 800 E. Lancaster Ave, Villanova, PA, 19805, USA.
| |
Collapse
|
11
|
Wang X, Wong LM, McElvain ME, Martire S, Lee WH, Li CZ, Fisher FA, Maheshwari RL, Wu ML, Imun MC, Murad R, Warshaviak DT, Yin J, Kamb A, Xu H. A rational approach to assess off-target reactivity of a dual-signal integrator for T cell therapy. Toxicol Appl Pharmacol 2022; 437:115894. [DOI: 10.1016/j.taap.2022.115894] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Revised: 01/15/2022] [Accepted: 01/19/2022] [Indexed: 01/16/2023]
|
12
|
Wangchuk J, Chatterjee A, Patil S, Madugula SK, Kondabagil K. The coevolution of large and small terminases of bacteriophages is a result of purifying selection leading to phenotypic stabilization. Virology 2021; 564:13-25. [PMID: 34598064 DOI: 10.1016/j.virol.2021.09.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2021] [Revised: 09/14/2021] [Accepted: 09/14/2021] [Indexed: 10/20/2022]
Abstract
Genome packaging in many dsDNA phages requires a series of precisely coordinated actions of two phage-coded proteins, namely, large terminase (TerL) and small terminase (TerS) with DNA and ATP, and with each other. Despite the strict functional conservation, TerL and TerS homologs exhibit large sequence variations. We investigated the sequence variability across eight phage types and observed a coevolutionary framework wherein the genealogy of TerL homologs mirrored that of the corresponding TerS homologs. Furthermore, a high purifying selection observed (dN/dS«1) indicated strong structural constraints on both TerL and TerS, and identify coevolving residues in TerL and TerS of phage T4 and lambda. Using the highly coevolving (correlation coefficient of 0.99) TerL and TerS of phage N4, we show that their biochemical features are similar to the phylogenetically divergent phage λ terminases. We also demonstrate using the Surface Plasma Resonance (SPR) technique that phage N4 TerL transiently interacts with TerS.
Collapse
Affiliation(s)
- Jigme Wangchuk
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, India
| | - Anirvan Chatterjee
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, India
| | - Supriya Patil
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, India
| | - Santhosh Kumar Madugula
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, India
| | - Kiran Kondabagil
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, India.
| |
Collapse
|
13
|
Rangel LT, Soucy SM, Setubal JC, Gogarten JP, Fournier GP. An efficient, non-phylogenetic method for detecting genes sharing evolutionary signals in phylogenomic datasets. Genome Biol Evol 2021; 13:6352501. [PMID: 34390574 PMCID: PMC8483891 DOI: 10.1093/gbe/evab187] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/11/2021] [Indexed: 11/25/2022] Open
Abstract
Assessing the compatibility between gene family phylogenies is a crucial and often computationally demanding step in many phylogenomic analyses. Here, we describe the Evolutionary Similarity Index (IES), a means to assess shared evolution between gene families using a weighted orthogonal distance regression model applied to sequence distances. The utilization of pairwise distance matrices circumvents comparisons between gene tree topologies, which are inherently uncertain and sensitive to evolutionary model choice, phylogenetic reconstruction artifacts, and other sources of error. Furthermore, IES enables the many-to-many pairing of multiple copies between similarly evolving gene families. This is done by selecting non-overlapping pairs of copies, one from each assessed family, and yielding the least sum of squared residuals. Analyses of simulated gene family data sets show that IES’s accuracy is on par with popular tree-based methods while also less susceptible to noise introduced by sequence alignment and evolutionary model fitting. Applying IES to an empirical data set of 1,322 genes from 42 archaeal genomes identified eight major clusters of gene families with compatible evolutionary trends. The most cohesive cluster consisted of 62 genes with compatible evolutionary signal, which occur as both single-copy and multiple homologs per genome; phylogenetic analysis of concatenated alignments from this cluster produced a tree closely matching previously published species trees for Archaea. Four other clusters are mainly composed of accessory genes with limited distribution among Archaea and enriched toward specific metabolic functions. Pairwise evolutionary distances obtained from these accessory gene clusters suggest patterns of interphyla horizontal gene transfer. An IES implementation is available at https://github.com/lthiberiol/evolSimIndex.
Collapse
Affiliation(s)
- Luiz Thibério Rangel
- Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
- Corresponding author: E-mail:
| | - Shannon M Soucy
- Department of Biomedical Data Science, Geisel School of Medicine, Dartmouth College, Hanover, New Hampshire, USA
| | - João C Setubal
- Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, Brasil
| | - Johann Peter Gogarten
- Department of Molecular and Cell Biology, University of Connecticut, USA
- Institute for Systems Genomics, University of Connecticut, USA
| | - Gregory P Fournier
- Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| |
Collapse
|
14
|
Forsythe ES, Williams AM, Sloan DB. Genome-wide signatures of plastid-nuclear coevolution point to repeated perturbations of plastid proteostasis systems across angiosperms. THE PLANT CELL 2021; 33:980-997. [PMID: 33764472 PMCID: PMC8226287 DOI: 10.1093/plcell/koab021] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Accepted: 01/16/2021] [Indexed: 05/05/2023]
Abstract
Nuclear and plastid (chloroplast) genomes experience different mutation rates, levels of selection, and transmission modes, yet key cellular functions depend on their coordinated interactions. Functionally related proteins often show correlated changes in rates of sequence evolution across a phylogeny [evolutionary rate covariation (ERC)], offering a means to detect previously unidentified suites of coevolving and cofunctional genes. We performed phylogenomic analyses across angiosperm diversity, scanning the nuclear genome for genes that exhibit ERC with plastid genes. As expected, the strongest hits were highly enriched for genes encoding plastid-targeted proteins, providing evidence that cytonuclear interactions affect rates of molecular evolution at genome-wide scales. Many identified nuclear genes functioned in post-transcriptional regulation and the maintenance of protein homeostasis (proteostasis), including protein translation (in both the plastid and cytosol), import, quality control, and turnover. We also identified nuclear genes that exhibit strong signatures of coevolution with the plastid genome, but their encoded proteins lack organellar-targeting annotations, making them candidates for having previously undescribed roles in plastids. In sum, our genome-wide analyses reveal that plastid-nuclear coevolution extends beyond the intimate molecular interactions within chloroplast enzyme complexes and may be driven by frequent rewiring of the machinery responsible for maintenance of plastid proteostasis in angiosperms.
Collapse
Affiliation(s)
- Evan S Forsythe
- Department of Biology, Colorado State University, Fort Collins, Colorado 80523, USA
| | - Alissa M Williams
- Department of Biology, Colorado State University, Fort Collins, Colorado 80523, USA
| | - Daniel B Sloan
- Department of Biology, Colorado State University, Fort Collins, Colorado 80523, USA
| |
Collapse
|
15
|
Patil S, Kondabagil K. Coevolutionary and Phylogenetic Analysis of Mimiviral Replication Machinery Suggest the Cellular Origin of Mimiviruses. Mol Biol Evol 2021; 38:2014-2029. [PMID: 33570580 PMCID: PMC8097291 DOI: 10.1093/molbev/msab003] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
Mimivirus is one of the most complex and largest viruses known. The origin and evolution of Mimivirus and other giant viruses have been a subject of intense study in the last two decades. The two prevailing hypotheses on the origin of Mimivirus and other viruses are the reduction hypothesis, which posits that viruses emerged from modern unicellular organisms; whereas the virus-first hypothesis proposes viruses as relics of precellular forms of life. In this study, to gain insights into the origin of Mimivirus, we have carried out extensive phylogenetic, correlation, and multidimensional scaling analyses of the putative proteins involved in the replication of its 1.2-Mb large genome. Correlation analysis and multidimensional scaling methods were validated using bacteriophage, bacteria, archaea, and eukaryotic replication proteins before applying to Mimivirus. We show that a large fraction of mimiviral replication proteins, including polymerase B, clamp, and clamp loaders are of eukaryotic origin and are coevolving. Although phylogenetic analysis places some components along the lineages of phage and bacteria, we show that all the replication-related genes have been homogenized and are under purifying selection. Collectively our analysis supports the idea that Mimivirus originated from a complex cellular ancestor. We hypothesize that Mimivirus has largely retained complex replication machinery reminiscent of its progenitor while losing most of the other genes related to processes such as metabolism and translation.
Collapse
Affiliation(s)
- Supriya Patil
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra, India
| | - Kiran Kondabagil
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra, India
| |
Collapse
|
16
|
Trivial and nontrivial error sources account for misidentification of protein partners in mutual information approaches. Sci Rep 2021; 11:6902. [PMID: 33767294 PMCID: PMC7994710 DOI: 10.1038/s41598-021-86455-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Accepted: 03/15/2021] [Indexed: 12/01/2022] Open
Abstract
The problem of finding the correct set of partners for a given pair of interacting protein families based on multi-sequence alignments (MSAs) has received great attention over the years. Recently, the native contacts of two interacting proteins were shown to store the strongest mutual information (MI) signal to discriminate MSA concatenations with the largest fraction of correct pairings. Although that signal might be of practical relevance in the search for an effective heuristic to solve the problem, the number of MSA concatenations with near-native MI is large, imposing severe limitations. Here, a Genetic Algorithm that explores possible MSA concatenations according to a MI maximization criteria is shown to find degenerate solutions with two error sources, arising from mismatches among (i) similar and (ii) non-similar sequences. If mistakes made among similar sequences are disregarded, type-(i) solutions are found to resolve correct pairings at best true positive (TP) rates of 70%—far above the very same estimates in type-(ii) solutions. A machine learning classification algorithm helps to show further that differences between optimized solutions based on TP rates are not artificial and may have biological meaning associated with the three-dimensional distribution of the MI signal. Type-(i) solutions may therefore correspond to reliable results for predictive purposes, found here to be more likely obtained via MI maximization across protein systems having a minimum critical number of amino acid contacts on their interaction surfaces (N > 200).
Collapse
|
17
|
Salmanian S, Pezeshk H, Sadeghi M. Inter-protein residue covariation information unravels physically interacting protein dimers. BMC Bioinformatics 2020; 21:584. [PMID: 33334319 PMCID: PMC7745481 DOI: 10.1186/s12859-020-03930-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2020] [Accepted: 12/09/2020] [Indexed: 01/04/2023] Open
Abstract
BACKGROUND Predicting physical interaction between proteins is one of the greatest challenges in computational biology. There are considerable various protein interactions and a huge number of protein sequences and synthetic peptides with unknown interacting counterparts. Most of co-evolutionary methods discover a combination of physical interplays and functional associations. However, there are only a handful of approaches which specifically infer physical interactions. Hybrid co-evolutionary methods exploit inter-protein residue coevolution to unravel specific physical interacting proteins. In this study, we introduce a hybrid co-evolutionary-based approach to predict physical interplays between pairs of protein families, starting from protein sequences only. RESULTS In the present analysis, pairs of multiple sequence alignments are constructed for each dimer and the covariation between residues in those pairs are calculated by CCMpred (Contacts from Correlated Mutations predicted) and three mutual information based approaches for ten accessible surface area threshold groups. Then, whole residue couplings between proteins of each dimer are unified into a single Frobenius norm value. Norms of residue contact matrices of all dimers in different accessible surface area thresholds are fed into support vector machine as single or multiple feature models. The results of training the classifiers by single features show no apparent different accuracies in distinct methods for different accessible surface area thresholds. Nevertheless, mutual information product and context likelihood of relatedness procedures may roughly have an overall higher and lower performances than other two methods for different accessible surface area cut-offs, respectively. The results also demonstrate that training support vector machine with multiple norm features for several accessible surface area thresholds leads to a considerable improvement of prediction performance. In this context, CCMpred roughly achieves an overall better performance than mutual information based approaches. The best accuracy, sensitivity, specificity, precision and negative predictive value for that method are 0.98, 1, 0.962, 0.96, and 0.962, respectively. CONCLUSIONS In this paper, by feeding norm values of protein dimers into support vector machines in different accessible surface area thresholds, we demonstrate that even small number of proteins in pairs of multiple alignments could allow one to accurately discriminate between positive and negative dimers.
Collapse
Affiliation(s)
- Sara Salmanian
- Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Hamid Pezeshk
- School of Mathematics, Statistics and Computer Science, College of Science, University of Tehran, Tehran, Iran
- Present Address: Department of Mathematics and Statistics, Concordia University, Montreal, Canada
- School of Biological Sciences, Institute for Research in Fundamental Sciences, Tehran, Iran
| | - Mehdi Sadeghi
- National Institute of Genetic Engineering and Biotechnology, Tehran, Iran
| |
Collapse
|
18
|
Dornburg A, Wang Z, Wang J, Mo ES, López-Giráldez F, Townsend JP. Comparative Genomics within and across Bilaterians Illuminates the Evolutionary History of ALK and LTK Proto-Oncogene Origination and Diversification. Genome Biol Evol 2020; 13:5983394. [PMID: 33196781 PMCID: PMC7851593 DOI: 10.1093/gbe/evaa228] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/23/2020] [Indexed: 12/14/2022] Open
Abstract
Comparative genomic analyses have enormous potential for identifying key genes central to human health phenotypes, including those that promote cancers. In particular, the successful development of novel therapeutics using model species requires phylogenetic analyses to determine molecular homology. Accordingly, we investigate the evolutionary histories of anaplastic lymphoma kinase (ALK)—which can underlie tumorigenesis in neuroblastoma, nonsmall cell lung cancer, and anaplastic large-cell lymphoma—its close relative leukocyte tyrosine kinase (LTK) and their candidate ligands. Homology of ligands identified in model organisms to those functioning in humans remains unclear. Therefore, we searched for homologs of the human genes across metazoan genomes, finding that the candidate ligands Jeb and Hen-1 were restricted to nonvertebrate species. In contrast, the ligand augmentor (AUG) was only identified in vertebrates. We found two ALK-like and four AUG-like protein-coding genes in lamprey. Of these six genes, only one ALK-like and two AUG-like genes exhibited early embryonic expression that parallels model mammal systems. Two copies of AUG are present in nearly all jawed vertebrates. Our phylogenetic analysis strongly supports the presence of previously unrecognized functional convergences of ALK and LTK between actinopterygians and sarcopterygians—despite contemporaneous, highly conserved synteny of ALK and LTK. These findings provide critical guidance regarding the propriety of fish and mammal models with regard to model organism-based investigation of these medically important genes. In sum, our results provide the phylogenetic context necessary for effective investigations of the functional roles and biology of these critically important receptors.
Collapse
Affiliation(s)
- Alex Dornburg
- Department of Bioinformatics and Genomics, University of North Carolina Charlotte
| | - Zheng Wang
- Department of Ecology and Evolutionary Biology, Yale University, New Haven.,Department of Biostatistics, Yale School of Public Health, New Haven, Connecticut
| | - Junrui Wang
- Department of Ecology and Evolutionary Biology, Yale University, New Haven.,Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Elizabeth S Mo
- Yale Combined Program in the Biological and Biomedical Sciences, Yale School of Medicine, Yale University, New Haven
| | | | - Jeffrey P Townsend
- Department of Ecology and Evolutionary Biology, Yale University, New Haven.,Department of Biostatistics, Yale School of Public Health, New Haven, Connecticut.,Program in Microbiology, Yale University, New Haven
| |
Collapse
|
19
|
Casas-Pastor D, Diehl A, Fritz G. Coevolutionary Analysis Reveals a Conserved Dual Binding Interface between Extracytoplasmic Function σ Factors and Class I Anti-σ Factors. mSystems 2020; 5:e00310-20. [PMID: 32753504 PMCID: PMC7406223 DOI: 10.1128/msystems.00310-20] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Accepted: 07/17/2020] [Indexed: 11/30/2022] Open
Abstract
Extracytoplasmic function σ factors (ECFs) belong to the most abundant signal transduction mechanisms in bacteria. Among the diverse regulators of ECF activity, class I anti-σ factors are the most important signal transducers in response to internal and external stress conditions. Despite the conserved secondary structure of the class I anti-σ factor domain (ASDI) that binds and inhibits the ECF under noninducing conditions, the binding interface between ECFs and ASDIs is surprisingly variable between the published cocrystal structures. In this work, we provide a comprehensive computational analysis of the ASDI protein family and study the different contact themes between ECFs and ASDIs. To this end, we harness the coevolution of these diverse protein families and predict covarying amino acid residues as likely candidates of an interaction interface. As a result, we find two common binding interfaces linking the first alpha-helix of the ASDI to the DNA-binding region in the σ4 domain of the ECF, and the fourth alpha-helix of the ASDI to the RNA polymerase (RNAP)-binding region of the σ2 domain. The conservation of these two binding interfaces contrasts with the apparent quaternary structure diversity of the ECF/ASDI complexes, partially explaining the high specificity between cognate ECF and ASDI pairs. Furthermore, we suggest that the dual inhibition of RNAP- and DNA-binding interfaces is likely a universal feature of other ECF anti-σ factors, preventing the formation of nonfunctional trimeric complexes between σ/anti-σ factors and RNAP or DNA.IMPORTANCE In the bacterial world, extracytoplasmic function σ factors (ECFs) are the most widespread family of alternative σ factors, mediating many cellular responses to environmental cues, such as stress. This work uses a computational approach to investigate how these σ factors interact with class I anti-σ factors-the most abundant regulators of ECF activity. By comprehensively classifying the anti-σs into phylogenetic groups and by comparing this phylogeny to the one of the cognate ECFs, the study shows how these protein families have coevolved to maintain their interaction over evolutionary time. These results shed light on the common contact residues that link ECFs and anti-σs in different phylogenetic families and set the basis for the rational design of anti-σs to specifically target certain ECFs. This will help to prevent the cross talk between heterologous ECF/anti-σ pairs, allowing their use as orthogonal regulators for the construction of genetic circuits in synthetic biology.
Collapse
Affiliation(s)
- Delia Casas-Pastor
- Center for Synthetic Microbiology (SYNMIKRO), Philipps-University Marburg, Marburg, Germany
| | - Angelika Diehl
- Center for Synthetic Microbiology (SYNMIKRO), Philipps-University Marburg, Marburg, Germany
- School of Molecular Sciences, University of Western Australia, Perth, Australia
| | - Georg Fritz
- Center for Synthetic Microbiology (SYNMIKRO), Philipps-University Marburg, Marburg, Germany
- School of Molecular Sciences, University of Western Australia, Perth, Australia
| |
Collapse
|
20
|
Gandarilla-Pérez CA, Mergny P, Weigt M, Bitbol AF. Statistical physics of interacting proteins: Impact of dataset size and quality assessed in synthetic sequences. Phys Rev E 2020; 101:032413. [PMID: 32290011 DOI: 10.1103/physreve.101.032413] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Accepted: 03/04/2020] [Indexed: 11/07/2022]
Abstract
Identifying protein-protein interactions is crucial for a systems-level understanding of the cell. Recently, algorithms based on inverse statistical physics, e.g., direct coupling analysis (DCA), have allowed to use evolutionarily related sequences to address two conceptually related inference tasks: finding pairs of interacting proteins and identifying pairs of residues which form contacts between interacting proteins. Here we address two underlying questions: How are the performances of both inference tasks related? How does performance depend on dataset size and the quality? To this end, we formalize both tasks using Ising models defined over stochastic block models, with individual blocks representing single proteins and interblock couplings protein-protein interactions; controlled synthetic sequence data are generated by Monte Carlo simulations. We show that DCA is able to address both inference tasks accurately when sufficiently large training sets of known interaction partners are available and that an iterative pairing algorithm allows to make predictions even without a training set. Noise in the training data deteriorates performance. In both tasks we find a quadratic scaling relating dataset quality and size that is consistent with noise adding in square-root fashion and signal adding linearly when increasing the dataset. This implies that it is generally good to incorporate more data even if their quality are imperfect, thereby shedding light on the empirically observed performance of DCA applied to natural protein sequences.
Collapse
Affiliation(s)
- Carlos A Gandarilla-Pérez
- Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire de Biologie Computationnelle et Quantitative (LCQB, UMR 7238), F-75005 Paris, France.,Facultad de Física, Universidad de la Habana, San Lázaro y L, Vedado, Habana 4, CP-10400, Cuba
| | - Pierre Mergny
- Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire de Biologie Computationnelle et Quantitative (LCQB, UMR 7238), F-75005 Paris, France.,Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire Jean Perrin (LJP, UMR 8237), F-75005 Paris, France
| | - Martin Weigt
- Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire de Biologie Computationnelle et Quantitative (LCQB, UMR 7238), F-75005 Paris, France
| | - Anne-Florence Bitbol
- Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire Jean Perrin (LJP, UMR 8237), F-75005 Paris, France.,Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
21
|
Dench J, Hinz A, Aris‐Brosou S, Kassen R. Identifying the drivers of computationally detected correlated evolution among sites under antibiotic selection. Evol Appl 2020; 13:781-793. [PMID: 32211067 PMCID: PMC7086105 DOI: 10.1111/eva.12900] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2019] [Revised: 10/02/2019] [Accepted: 11/14/2019] [Indexed: 11/29/2022] Open
Abstract
The ultimate causes of correlated evolution among sites in a genome remain difficult to tease apart. To address this problem directly, we performed a high-throughput search for correlated evolution among sites associated with resistance to a fluoroquinolone antibiotic using whole-genome data from clinical strains of Pseudomonas aeruginosa, before validating our computational predictions experimentally. We show that for at least two sites, this correlation is underlain by epistasis. Our analysis also revealed eight additional pairs of synonymous substitutions displaying correlated evolution underlain by physical linkage, rather than selection associated with antibiotic resistance. Our results provide direct evidence that both epistasis and physical linkage among sites can drive the correlated evolution identified by high-throughput computational tools. In other words, the observation of correlated evolution is not by itself sufficient evidence to guarantee that the sites in question are epistatic; such a claim requires additional evidence, ideally coming from direct estimates of epistasis, based on experimental evidence.
Collapse
Affiliation(s)
- Jonathan Dench
- Department of BiologyUniversity of OttawaOttawaOntarioCanada
| | - Aaron Hinz
- Department of BiologyUniversity of OttawaOttawaOntarioCanada
| | - Stéphane Aris‐Brosou
- Department of BiologyUniversity of OttawaOttawaOntarioCanada
- Department of Mathematics and StatisticsUniversity of OttawaOttawaOntarioCanada
| | - Rees Kassen
- Department of BiologyUniversity of OttawaOttawaOntarioCanada
| |
Collapse
|
22
|
Coevolutive, evolutive and stochastic information in protein-protein interactions. Comput Struct Biotechnol J 2019; 17:1429-1435. [PMID: 31871588 PMCID: PMC6906720 DOI: 10.1016/j.csbj.2019.10.005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Revised: 10/19/2019] [Accepted: 10/22/2019] [Indexed: 11/24/2022] Open
Abstract
Here, we investigate the contributions of coevolutive, evolutive and stochastic information in determining protein-protein interactions (PPIs) based on primary sequences of two interacting protein families A and B. Specifically, under the assumption that coevolutive information is imprinted on the interacting amino acids of two proteins in contrast to other (evolutive and stochastic) sources spread over their sequences, we dissect those contributions in terms of compensatory mutations at physically-coupled and uncoupled amino acids of A and B. We find that physically-coupled amino-acids at short range distances store the largest per-contact mutual information content, with a significant fraction of that content resulting from coevolutive sources alone. The information stored in coupled amino acids is shown further to discriminate multi-sequence alignments (MSAs) with the largest expectation fraction of PPI matches – a conclusion that holds against various definitions of intermolecular contacts and binding modes. When compared to the informational content resulting from evolution at long-range interactions, the mutual information in physically-coupled amino-acids is the strongest signal to distinguish PPIs derived from cospeciation and likely, the unique indication in case of molecular coevolution in independent genomes as the evolutive information must vanish for uncorrelated proteins.
Collapse
|
23
|
Phylogenetic correlations can suffice to infer protein partners from sequences. PLoS Comput Biol 2019; 15:e1007179. [PMID: 31609984 PMCID: PMC6812855 DOI: 10.1371/journal.pcbi.1007179] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2019] [Revised: 10/24/2019] [Accepted: 09/25/2019] [Indexed: 12/30/2022] Open
Abstract
Determining which proteins interact together is crucial to a systems-level understanding of the cell. Recently, algorithms based on Direct Coupling Analysis (DCA) pairwise maximum-entropy models have allowed to identify interaction partners among paralogous proteins from sequence data. This success of DCA at predicting protein-protein interactions could be mainly based on its known ability to identify pairs of residues that are in contact in the three-dimensional structure of protein complexes and that coevolve to remain physicochemically complementary. However, interacting proteins possess similar evolutionary histories. What is the role of purely phylogenetic correlations in the performance of DCA-based methods to infer interaction partners? To address this question, we employ controlled synthetic data that only involve phylogeny and no interactions or contacts. We find that DCA accurately identifies the pairs of synthetic sequences that share evolutionary history. While phylogenetic correlations confound the identification of contacting residues by DCA, they are thus useful to predict interacting partners among paralogs. We find that DCA performs as well as phylogenetic methods to this end, and slightly better than them with large and accurate training sets. Employing DCA or phylogenetic methods within an Iterative Pairing Algorithm (IPA) allows to predict pairs of evolutionary partners without a training set. We further demonstrate the ability of these various methods to correctly predict pairings among real paralogous proteins with genome proximity but no known direct physical interaction, illustrating the importance of phylogenetic correlations in natural data. However, for physically interacting and strongly coevolving proteins, DCA and mutual information outperform phylogenetic methods. We finally discuss how to distinguish physically interacting proteins from proteins that only share a common evolutionary history. Many biologically important protein-protein interactions are conserved over evolutionary time scales. This leads to two different signals that can be used to computationally predict interactions between protein families and to identify specific interaction partners. First, the shared evolutionary history leads to highly similar phylogenetic relationships between interacting proteins of the two families. Second, the need to keep the interaction surfaces of partner proteins biophysically compatible causes a correlated amino-acid usage of interface residues. Employing simulated data, we show that the shared history alone can be used to detect partner proteins. Similar accuracies are achieved by algorithms comparing phylogenetic relationships and by methods based on Direct Coupling Analysis (DCA), which are primarily known for their ability to detect the second type of signal. Using natural sequence data, we show that in cases with shared evolutionary history but without known physical interactions, both methods work with similar accuracy, while for some physically interacting systems, DCA and mutual information outperform phylogenetic methods. We propose methods allowing both to predict interactions between protein families and to find interacting partners among paralogs.
Collapse
|
24
|
Wu Q, Lan Y, Cao X, Yao H, Qiao D, Xu H, Cao Y. Characterization and diverse evolution patterns of glycerol-3-phosphate dehydrogenase family genes in Dunaliella salina. Gene 2019; 710:161-169. [DOI: 10.1016/j.gene.2019.05.056] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2018] [Revised: 05/15/2019] [Accepted: 05/29/2019] [Indexed: 12/29/2022]
|
25
|
Wang X, Yu B, Ma A, Chen C, Liu B, Ma Q. Protein-protein interaction sites prediction by ensemble random forests with synthetic minority oversampling technique. Bioinformatics 2019; 35:2395-2402. [PMID: 30520961 PMCID: PMC6612859 DOI: 10.1093/bioinformatics/bty995] [Citation(s) in RCA: 86] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2018] [Revised: 11/19/2018] [Accepted: 12/03/2018] [Indexed: 11/14/2022] Open
Abstract
MOTIVATION The prediction of protein-protein interaction (PPI) sites is a key to mutation design, catalytic reaction and the reconstruction of PPI networks. It is a challenging task considering the significant abundant sequences and the imbalance issue in samples. RESULTS A new ensemble learning-based method, Ensemble Learning of synthetic minority oversampling technique (SMOTE) for Unbalancing samples and RF algorithm (EL-SMURF), was proposed for PPI sites prediction in this study. The sequence profile feature and the residue evolution rates were combined for feature extraction of neighboring residues using a sliding window, and the SMOTE was applied to oversample interface residues in the feature space for the imbalance problem. The Multi-dimensional Scaling feature selection method was implemented to reduce feature redundancy and subset selection. Finally, the Random Forest classifiers were applied to build the ensemble learning model, and the optimal feature vectors were inserted into EL-SMURF to predict PPI sites. The performance validation of EL-SMURF on two independent validation datasets showed 77.1% and 77.7% accuracy, which were 6.2-15.7% and 6.1-18.9% higher than the other existing tools, respectively. AVAILABILITY AND IMPLEMENTATION The source codes and data used in this study are publicly available at http://github.com/QUST-AIBBDRC/EL-SMURF/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Xiaoying Wang
- College of Mathematics and Physics, Qingdao University of Science and Technology, Qingdao, China
- School of Mathematics, Shandong University, Jinan, China
- Artificial Intelligence and Biomedical Big Data Research Center, Qingdao University of Science and Technology, Qingdao, China
| | - Bin Yu
- College of Mathematics and Physics, Qingdao University of Science and Technology, Qingdao, China
- Artificial Intelligence and Biomedical Big Data Research Center, Qingdao University of Science and Technology, Qingdao, China
- School of Life Sciences, University of Science and Technology of China, Hefei, China
| | - Anjun Ma
- Bioinformatics and Mathematical Biosciences Lab, Department of Agronomy, Horticulture and Plant Science, South Dakota State University, Brookings, SD, USA
- Department Biomedical Informatics, College of Medicine, The Ohio State University, Columbus, OH, USA
| | - Cheng Chen
- College of Mathematics and Physics, Qingdao University of Science and Technology, Qingdao, China
- Artificial Intelligence and Biomedical Big Data Research Center, Qingdao University of Science and Technology, Qingdao, China
| | - Bingqiang Liu
- School of Mathematics, Shandong University, Jinan, China
| | - Qin Ma
- Bioinformatics and Mathematical Biosciences Lab, Department of Agronomy, Horticulture and Plant Science, South Dakota State University, Brookings, SD, USA
- Department Biomedical Informatics, College of Medicine, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
26
|
Ding Z, Kihara D. Computational identification of protein-protein interactions in model plant proteomes. Sci Rep 2019; 9:8740. [PMID: 31217453 PMCID: PMC6584649 DOI: 10.1038/s41598-019-45072-8] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2019] [Accepted: 05/30/2019] [Indexed: 12/12/2022] Open
Abstract
Protein-protein interactions (PPIs) play essential roles in many biological processes. A PPI network provides crucial information on how biological pathways are structured and coordinated from individual protein functions. In the past two decades, large-scale PPI networks of a handful of organisms were determined by experimental techniques. However, these experimental methods are time-consuming, expensive, and are not easy to perform on new target organisms. Large-scale PPI data is particularly sparse in plant organisms. Here, we developed a computational approach for detecting PPIs trained and tested on known PPIs of Arabidopsis thaliana and applied to three plants, Arabidopsis thaliana, Glycine max (soybean), and Zea mays (maize) to discover new PPIs on a genome-scale. Our method considers a variety of features including protein sequences, gene co-expression, functional association, and phylogenetic profiles. This is the first work where a PPI prediction method was developed for is the first PPI prediction method applied on benchmark datasets of Arabidopsis. The method showed a high prediction accuracy of over 90% and very high precision of close to 1.0. We predicted 50,220 PPIs in Arabidopsis thaliana, 13,175,414 PPIs in corn, and 13,527,834 PPIs in soybean. Newly predicted PPIs were classified into three confidence levels according to the availability of existing supporting evidence and discussed. Predicted PPIs in the three plant genomes are made available for future reference.
Collapse
Affiliation(s)
- Ziyun Ding
- Department of Biological Sciences, Purdue University, West Lafayette, IN, 47907, USA.
| | - Daisuke Kihara
- Department of Biological Sciences, Purdue University, West Lafayette, IN, 47907, USA.
- Department of Computer Science, Purdue University, West Lafayette, IN, 47907, USA.
- Department of Pediatrics, University of Cincinnati, Cincinnati, OH, 45229, USA.
| |
Collapse
|
27
|
Pathogenicity of the H1N1 influenza virus enhanced by functional synergy between the NPV100I and NAD248N pair. PLoS One 2019; 14:e0217691. [PMID: 31150476 PMCID: PMC6544299 DOI: 10.1371/journal.pone.0217691] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2018] [Accepted: 05/16/2019] [Indexed: 11/20/2022] Open
Abstract
By comparing and measuring covariations of viral protein sequences from isolates of the 2009 pH1N1 influenza A virus (IAV), specific substitutions that co-occur in the NP-NA pair were identified. To investigate the effect of these co-occurring substitution pairs, the V100I substitution in NP and the D248N substitution in NA were introduced into laboratory-adapted WSN IAVs. The recombinant WSN with the covarying NPV100I-NAD248N pair exhibited enhanced pathogenicity, as characterized by increased viral production, increased death and inflammation of host cells, and high mortality in infected mice. Although direct interactions between the NPV100I and NAD248N proteins were not detected, the RNA-binding ability of NPV100I was increased, which was further strengthened by NAD248N, in expression-plasmid-transfected cells. Additionally, the NAD248N protein was frequently recruited within lipid rafts, indirectly affecting the RNA-binding ability of NP as well as viral release. Altogether, our data indicate that the covarying NPV100I-NAD248N pair obtained from 2009 pH1N1 IAV sequence information function together to synergistically augment viral assembly and release, which may explain the observed enhanced viral pathogenicity.
Collapse
|
28
|
Pandey S, Phale PS, Bhaumik P. Structural modulation of a periplasmic sugar-binding protein probes into its evolutionary ancestry. J Struct Biol 2018; 204:498-506. [PMID: 30244006 DOI: 10.1016/j.jsb.2018.09.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Revised: 09/10/2018] [Accepted: 09/19/2018] [Indexed: 10/28/2022]
Abstract
Substrate-binding proteins (SBPs) are periplasmic proteins consisting of two α/β domains joined by a hinge region with specificity towards cognate ligands. Based on three-dimensional fold, sugar-specific SBPs have been classified into cluster B and cluster D-I. The analysis of sequences and structures of sugar-binding pocket of cluster D-I SBPs revealed the presence of extra residues on two loops (L1, L2) and a helix (H1) in few members of this family, that binds specifically to monosaccharides. Presence of conserved histidine in L2 and tryptophan in H1 can be considered as the identity marks for the cluster D-I monosaccharide-binding SBPs. A glucose binding protein (ppGBP) from Pseudomonas putida CSV86 was found to contain a structural fold similar to oligosaccharide-binding cluster D-I SBPs, but functionally binds to only glucose due to constriction of its binding pocket mainly by L2 (375-382). ppGBP with partial deletion of L2 (ppGBPΔL2) was created, crystallized and biochemical characterization was performed. Compared to wild type ppGBP, the ppGBPΔL2 structure showed widening of the glucose-binding pocket with ∼80% lower glucose binding. Our results show that the substrate specificity of SBPs can be altered by modulating the size of the binding pocket. Based on this, we propose a sub classification of cluster D-I SBPs into (i) cluster D-I(a)-monosaccharide-binding SBPs and (ii) cluster D-I(b)-oligosaccharide-binding SBPs. This study also provides the direct structural and functional correlation indicating that divergence of proteins may occur through insertions or deletions of sequences in the already existing SBPs leading to evolution at the functional level.
Collapse
Affiliation(s)
- Suman Pandey
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra 400076, India
| | - Prashant S Phale
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra 400076, India
| | - Prasenjit Bhaumik
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra 400076, India.
| |
Collapse
|
29
|
Hersh TA, Dimond AL, Ruth BA, Lupica NV, Bruce JC, Kelley JM, King BL, Lutton BV. A role for the CXCR4-CXCL12 axis in the little skate, Leucoraja erinacea. Am J Physiol Regul Integr Comp Physiol 2018; 315:R218-R229. [PMID: 29641231 PMCID: PMC6139610 DOI: 10.1152/ajpregu.00322.2017] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The interaction between C-X-C chemokine receptor type 4 (CXCR4) and its cognate ligand C-X-C motif chemokine ligand 12 (CXCL12) plays a critical role in regulating hematopoietic stem cell activation and subsequent cellular mobilization. Extensive studies of these genes have been conducted in mammals, but much less is known about the expression and function of CXCR4 and CXCL12 in non-mammalian vertebrates. In the present study, we identify simultaneous expression of CXCR4 and CXCL12 orthologs in the epigonal organ (the primary hematopoietic tissue) of the little skate, Leucoraja erinacea. Genetic and phylogenetic analyses were functionally supported by significant mobilization of leukocytes following administration of Plerixafor, a CXCR4 antagonist and clinically important drug. Our results provide evidence that, as in humans, Plerixafor disrupts CXCR4/CXCL12 binding in the little skate, facilitating release of leukocytes into the bloodstream. Our study illustrates the value of the little skate as a model organism, particularly in studies of hematopoiesis and potentially for preclinical research on hematological and vascular disorders.
Collapse
Affiliation(s)
- Taylor A Hersh
- Mount Desert Island Biological Laboratory , Bar Harbor, Maine
| | - Alexandria L Dimond
- Mount Desert Island Biological Laboratory , Bar Harbor, Maine
- School of Arts and Sciences, Endicott College , Beverly, Massachusetts
| | - Brittany A Ruth
- Mount Desert Island Biological Laboratory , Bar Harbor, Maine
- School of Arts and Sciences, Endicott College , Beverly, Massachusetts
| | - Noah V Lupica
- Mount Desert Island Biological Laboratory , Bar Harbor, Maine
| | - Jacob C Bruce
- Mount Desert Island Biological Laboratory , Bar Harbor, Maine
| | - John M Kelley
- School of Arts and Sciences, Endicott College , Beverly, Massachusetts
- Beth Israel Deaconess Medical Center, Program in Placebo Studies, Harvard Medical School , Boston, Massachusetts
| | - Benjamin L King
- Department of Molecular and Biomedical Sciences, University of Maine , Orono, Maine
| | - Bram V Lutton
- Mount Desert Island Biological Laboratory , Bar Harbor, Maine
- School of Arts and Sciences, Endicott College , Beverly, Massachusetts
| |
Collapse
|
30
|
Ding Z, Kihara D. Computational Methods for Predicting Protein-Protein Interactions Using Various Protein Features. CURRENT PROTOCOLS IN PROTEIN SCIENCE 2018; 93:e62. [PMID: 29927082 PMCID: PMC6097941 DOI: 10.1002/cpps.62] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Understanding protein-protein interactions (PPIs) in a cell is essential for learning protein functions, pathways, and mechanism of diseases. PPIs are also important targets for developing drugs. Experimental methods, both small-scale and large-scale, have identified PPIs in several model organisms. However, results cover only a part of PPIs of organisms; moreover, there are many organisms whose PPIs have not yet been investigated. To complement experimental methods, many computational methods have been developed that predict PPIs from various characteristics of proteins. Here we provide an overview of literature reports to classify computational PPI prediction methods that consider different features of proteins, including protein sequence, genomes, protein structure, function, PPI network topology, and those which integrate multiple methods. © 2018 by John Wiley & Sons, Inc.
Collapse
Affiliation(s)
- Ziyun Ding
- Department of Biological Science, Purdue University, West Lafayette, IN, 47907 USA
| | - Daisuke Kihara
- Department of Biological Science, Purdue University, West Lafayette, IN, 47907 USA
- Department of Computer Science, Purdue University, West Lafayette, IN, 47907 USA
- Corresponding author: DK; , Phone: 1-765-496-2284 (DK)
| |
Collapse
|
31
|
Gribouval L, Sourdaine P, Lareyre JJ, Bellaiche J, Le Gac F, Mazan S, Guiardiere C, Auvray P, Gautier A. The nanos1 gene was duplicated in early Vertebrates and the two paralogs show different gonadal expression profiles in a shark. Sci Rep 2018; 8:6942. [PMID: 29720681 PMCID: PMC5932020 DOI: 10.1038/s41598-018-24643-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2017] [Accepted: 04/04/2018] [Indexed: 11/23/2022] Open
Abstract
Nanos are RNA-binding proteins playing crucial roles in germ cell development and maintenance. Based on phylogenetic and synteny analyses, this study reveals that nanos1 gene has undergone multiple duplications and gene copies losses in Vertebrates. Chondrichthyan species display two nanos1 genes (named nanos1A/1B), which were both retrieved in some Osteichthyes at basal positions in Sarcopterygii and Actinopterygii lineages. In contrast, Teleosts have lost nanos1A but duplicated nanos1B leading to the emergence of two ohnologs (nanos1Ba/1Bb), whereas Tetrapods have lost nanos1B gene. The two successive nanos gene duplications may result from the second and third whole genome duplication events at the basis of Vertebrates and Teleosts respectively. The expression profiles of nanos1A and nanos1B paralogs were characterized in the dogfish, Scyliorhinus canicula. Nanos1A was strongly expressed in brain and also localized in all germ cell types in the polarized testis. In contrast, nanos1B was detected in testis with the highest expression in the germinative zone. In addition, Nanos1B protein was predominantly located in the nuclei of male germinal cells. In the ovary, both paralogs were detected in germinal and somatic cells. Our study opens new perspectives concerning the complex evolution of nanos1 paralogs and their potential distinct roles in Vertebrates gonads.
Collapse
Affiliation(s)
- Laura Gribouval
- Normandie University, UNICAEN, Sorbonne Universités, MNHN, UPMC University Paris 06, UA, CNRS, IRD, Biologie des Organismes et Ecosystèmes Aquatiques (BOREA), CS14032, 14032 CAEN, Cedex 5, France
- KELIA, Parc Technopolitain Atalante Saint Malo, 35400, Saint Malo, France
| | - Pascal Sourdaine
- Normandie University, UNICAEN, Sorbonne Universités, MNHN, UPMC University Paris 06, UA, CNRS, IRD, Biologie des Organismes et Ecosystèmes Aquatiques (BOREA), CS14032, 14032 CAEN, Cedex 5, France
| | - Jean-Jacques Lareyre
- INRA UPR1037, Laboratory of Fish Physiology and Genomics, BIOSIT, Ouest-Genopole, Campus de Beaulieu, 35042, Rennes, France
| | - Johanna Bellaiche
- INRA UPR1037, Laboratory of Fish Physiology and Genomics, BIOSIT, Ouest-Genopole, Campus de Beaulieu, 35042, Rennes, France
| | - Florence Le Gac
- INRA UPR1037, Laboratory of Fish Physiology and Genomics, BIOSIT, Ouest-Genopole, Campus de Beaulieu, 35042, Rennes, France
| | - Sylvie Mazan
- CNRS-UPMC-Sorbonne Universités, UMR 7232, Observatoire océanologique, 66650, Banyuls sur mer, France
| | - Cécile Guiardiere
- KELIA, Parc Technopolitain Atalante Saint Malo, 35400, Saint Malo, France
| | - Pierrïck Auvray
- KELIA, Parc Technopolitain Atalante Saint Malo, 35400, Saint Malo, France
| | - Aude Gautier
- Normandie University, UNICAEN, Sorbonne Universités, MNHN, UPMC University Paris 06, UA, CNRS, IRD, Biologie des Organismes et Ecosystèmes Aquatiques (BOREA), CS14032, 14032 CAEN, Cedex 5, France.
| |
Collapse
|
32
|
Liu L, Wu Y, Liao Z, Xiong J, Wu F, Xu J, Lan H, Tang Q, Zhou S, Liu Y, Lu Y. Evolutionary conservation and functional divergence of the LFK gene family play important roles in the photoperiodic flowering pathway of land plants. Heredity (Edinb) 2017; 120:310-328. [PMID: 29225355 DOI: 10.1038/s41437-017-0006-5] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2017] [Revised: 08/15/2017] [Accepted: 08/22/2017] [Indexed: 12/22/2022] Open
Abstract
ZEITLUPE (ZTL), LOV KELCH PROTEIN 2 (LKP2), and FLAVIN-BINDING KELCH REPEAT F-BOX 1 (FKF1)-blue-light photoreceptors-play important roles in regulating the circadian clock and photoperiodic flowering pathway in plants. In this study, phylogenetic analysis revealed that the LOV (Light, Oxygen, or Voltage) and Kelch repeat-containing F-box (LFK) gene family can be classified into two clades, ZTL/LKP2 and FKF1, with clear differentiation between monocots and dicots within each clade. The LFK family genes underwent strong purifying selection; however, signatures of positive selection to adapt to local conditions still existed in 18 specific codons. In 87 diverse maize inbred lines, significant differences were identified (P ≤ 0.01) for days to female flowering between the haplotypes consisting of eight positive selection sites at ZmFKF1b corresponding to tropical and temperate maize groups of the phylogenetic tree, indicating a key role of ZmFKF1b in maize adaptive evolution. In addition, positive coevolution was detected in the domains of the LFK family for long-term cooperation to targets. The Type-I and Type-II functional divergence analysis revealed subfunctionalization or neofunctionalization of the LFKs, and the ZTL subfamily is most likely to maintain the ancestral function of LFKs. Over 50% of critical amino acid sites involved in the functional divergence were identified in the Kelch repeat domain, resulting in the distinction of substrates for ubiquitination and degradation. These results suggest that evolutionary conservation contributes to the maintenance of critical physiological functions, whereas functional divergence after duplication helps to generate diverse molecular regulation mechanisms.
Collapse
Affiliation(s)
- Ling Liu
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Yuanqi Wu
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Zhengqiao Liao
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Jing Xiong
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Fengkai Wu
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Jie Xu
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Hai Lan
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Qiling Tang
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Shufeng Zhou
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China.,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China
| | - Yaxi Liu
- Triticeae Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China
| | - Yanli Lu
- Maize Research Institute, Sichuan Agricultural University, 611130, Wenjiang, Sichuan, China. .,Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Wenjiang, China.
| |
Collapse
|
33
|
Gupta A, Patil S, Vijayakumar R, Kondabagil K. The Polyphyletic Origins of Primase-Helicase Bifunctional Proteins. J Mol Evol 2017; 85:188-204. [PMID: 29143083 DOI: 10.1007/s00239-017-9816-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2017] [Accepted: 10/28/2017] [Indexed: 11/28/2022]
Abstract
We studied the evolutionary relationships of different primase-helicase bifunctional proteins, found mostly in viruses, virophages, plasmids, and organellar genomes, by phylogeny and correlation analysis. Our study suggests independent origins of primase-helicase bifunctional proteins resulting from multiple fusion events between genes encoding primase and helicase domains of different families. The correlation analysis further indicated strong functional dependencies of domains in the bifunctional proteins that are part of smaller genomes and plasmids. Bifunctional proteins found in some bacterial genomes exhibited weak coevolution probably suggesting that these are the non-functional remnants of the proteins acquired via horizontal transfer. We have put forward possible scenarios for the origin of primase-helicase bifunctional proteins in large eukaryotic DNA viruses and virophages.
Collapse
Affiliation(s)
- Ankita Gupta
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| | - Supriya Patil
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| | - Ramya Vijayakumar
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| | - Kiran Kondabagil
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India.
| |
Collapse
|
34
|
Abstract
The study of evolutionary relationships among protein sequences was one of the first applications of bioinformatics. Since then, and accompanying the wealth of biological data produced by genome sequencing and other high-throughput techniques, the use of bioinformatics in general and phylogenetics in particular has been gaining ground in the study of protein and proteome evolution. Nowadays, the use of phylogenetics is instrumental not only to infer the evolutionary relationships among species and their genome sequences, but also to reconstruct ancestral states of proteins and proteomes and hence trace the paths followed by evolution. Here I survey recent progress in the elucidation of mechanisms of protein and proteome evolution in which phylogenetics has played a determinant role.
Collapse
Affiliation(s)
- Toni Gabaldón
- Bioinformatics Department, Centro de Investigación Principe Felipe
| |
Collapse
|
35
|
Meysman P, Titeca K, Eyckerman S, Tavernier J, Goethals B, Martens L, Valkenborg D, Laukens K. Protein complex analysis: From raw protein lists to protein interaction networks. MASS SPECTROMETRY REVIEWS 2017; 36:600-614. [PMID: 26709718 DOI: 10.1002/mas.21485] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/04/2015] [Accepted: 11/17/2015] [Indexed: 06/05/2023]
Abstract
The elucidation of molecular interaction networks is one of the pivotal challenges in the study of biology. Affinity purification-mass spectrometry and other co-complex methods have become widely employed experimental techniques to identify protein complexes. These techniques typically suffer from a high number of false negatives and false positive contaminants due to technical shortcomings and purification biases. To support a diverse range of experimental designs and approaches, a large number of computational methods have been proposed to filter, infer and validate protein interaction networks from experimental pull-down MS data. Nevertheless, this expansion of available methods complicates the selection of the most optimal ones to support systems biology-driven knowledge extraction. In this review, we give an overview of the most commonly used computational methods to process and interpret co-complex results, and we discuss the issues and unsolved problems that still exist within the field. © 2015 Wiley Periodicals, Inc. Mass Spec Rev 36:600-614, 2017.
Collapse
Affiliation(s)
- Pieter Meysman
- Advanced Database Research and Modelling (ADReM), Department of Mathematics and Computer Science, University of Antwerp, Antwerp, Belgium
- Biomedical Informatics Research Center Antwerp (biomina), University of Antwerp/Antwerp University Hospital, Edegem, Belgium
| | - Kevin Titeca
- Department of Medical Protein Research, VIB, B-9000 Ghent, Belgium
- Department of Biochemistry, Ghent University, B-9000 Ghent, Belgium
| | - Sven Eyckerman
- Department of Medical Protein Research, VIB, B-9000 Ghent, Belgium
- Department of Biochemistry, Ghent University, B-9000 Ghent, Belgium
| | - Jan Tavernier
- Department of Medical Protein Research, VIB, B-9000 Ghent, Belgium
- Department of Biochemistry, Ghent University, B-9000 Ghent, Belgium
| | - Bart Goethals
- Advanced Database Research and Modelling (ADReM), Department of Mathematics and Computer Science, University of Antwerp, Antwerp, Belgium
| | - Lennart Martens
- Department of Medical Protein Research, VIB, B-9000 Ghent, Belgium
- Department of Biochemistry, Ghent University, B-9000 Ghent, Belgium
| | - Dirk Valkenborg
- Flemish Institute for Technological Research (VITO), Mol, Belgium
- IBioStat, Hasselt University, Hasselt, Belgium
- CFP-CeProMa, University of Antwerp, Antwerp, Belgium
| | - Kris Laukens
- Advanced Database Research and Modelling (ADReM), Department of Mathematics and Computer Science, University of Antwerp, Antwerp, Belgium
- Biomedical Informatics Research Center Antwerp (biomina), University of Antwerp/Antwerp University Hospital, Edegem, Belgium
| |
Collapse
|
36
|
Peng J, Wu Y, Tian X, Pang J, Kuai L, Cao F, Qin X, Zhong J, Li X, Li Y, Sun X, Chen L, Jiang Y. High-Throughput Sequencing and Co-Expression Network Analysis of lncRNAs and mRNAs in Early Brain Injury Following Experimental Subarachnoid Haemorrhage. Sci Rep 2017; 7:46577. [PMID: 28417961 PMCID: PMC5394545 DOI: 10.1038/srep46577] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2016] [Accepted: 03/17/2017] [Indexed: 12/22/2022] Open
Abstract
Subarachnoid haemorrhage (SAH) is a fatal neurovascular disease following cerebral aneurysm rupture with high morbidity and mortality rates. Long non-coding RNAs (lncRNAs) are a type of mammalian genome transcript, are abundantly expressed in the brain and are involved in many nervous system diseases. However, little is currently known regarding the influence of lncRNAs in early brain injury (EBI) after SAH. This study analysed the expression profiles of lncRNAs and mRNAs in SAH brain tissues of mice using high-throughput sequencing. The results showed a remarkable difference in lncRNA and mRNA transcripts between SAH and control brains. Approximately 617 lncRNA transcripts and 441 mRNA transcripts were aberrantly expressed at 24 hours after SAH. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that the differentially expressed mRNAs were mostly involved in inflammation. Based on the lncRNA/mRNA co-expression network, knockdown of fantom3_F730004F19 reduced the mRNA and protein levels of CD14 and toll-like receptor 4 (TLR4) and attenuated inflammation in BV-2 microglia cells. These results indicate that lncRNA fantom3_F730004F19 may be associated with microglia induced inflammation via the TLR signaling pathway in EBI following SAH. LncRNA represent a potential therapeutic target for the prognosis, diagnosis, and treatment of SAH.
Collapse
Affiliation(s)
- Jianhua Peng
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, China
| | - Yue Wu
- Department of Neurosurgery, the First Affiliated Hospital of Chongqing Medical University, Chongqing, China
| | - Xiaocui Tian
- Chongqing Key Laboratory of Biochemistry and Molecular Pharmacology, College of Pharmacy, Chongqing, China
| | - Jinwei Pang
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, China
| | - Li Kuai
- Department of Ophthalmology, the Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan, China
| | - Fang Cao
- Department of Neurovascular Disease, the Affiliated Hospital of Zunyi Medical College, Zunyi, China
| | - Xinghu Qin
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, China
- Department of Neurosurgery, People’s Hospital of Deyang City, Deyang, China
| | - Jianjun Zhong
- Department of Neurosurgery, the First Affiliated Hospital of Chongqing Medical University, Chongqing, China
| | - Xinshen Li
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, China
| | - Yong Li
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, China
| | - Xiaochuan Sun
- Department of Neurosurgery, the First Affiliated Hospital of Chongqing Medical University, Chongqing, China
| | - Ligang Chen
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, China
| | - Yong Jiang
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, China
| |
Collapse
|
37
|
Mandloi S, Chakrabarti S. Protein sites with more coevolutionary connections tend to evolve slower, while more variable protein families acquire higher coevolutionary connections. F1000Res 2017; 6:453. [PMID: 28751967 PMCID: PMC5506539 DOI: 10.12688/f1000research.11251.2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 07/05/2017] [Indexed: 11/20/2022] Open
Abstract
Background: Amino acid exchanges within proteins sometimes compensate for one another and could therefore be co-evolved. It is essential to investigate the intricate relationship between the extent of coevolution and the evolutionary variability exerted at individual protein sites, as well as the whole protein. Methods: In this study, we have used a reliable set of coevolutionary connections (sites within 10Å spatial distance) and investigated their correlation with the evolutionary diversity within the respective protein sites. Results: Based on our observations, we propose an interesting hypothesis that higher numbers of coevolutionary connections are associated with lesser evolutionary variable protein sites, while higher numbers of the coevolutionary connections can be observed for a protein family that has higher evolutionary variability. Our findings also indicate that highly coevolved sites located in a solvent accessible state tend to be less evolutionary variable. This relationship reverts at the whole protein level where cytoplasmic and extracellular proteins show moderately higher anti-correlation between the number of coevolutionary connections and the average evolutionary conservation of the whole protein. Conclusions: Observations and hypothesis presented in this study provide intriguing insights towards understanding the critical relationship between coevolutionary and evolutionary changes observed within proteins. Our observations encourage further investigation to find out the reasons behind subtle variations in the relationship between coevolutionary connectivity and evolutionary diversity for proteins located at various cellular localizations and/or involved in different molecular-biological functions.
Collapse
Affiliation(s)
- Sapan Mandloi
- Department of Structural Biology and Bioinformatics Division, Council of Scientific and Industrial Research, Indian Institute of Chemical Biology, Kolkata, West Bengal, 700032, India
| | - Saikat Chakrabarti
- Department of Structural Biology and Bioinformatics Division, Council of Scientific and Industrial Research, Indian Institute of Chemical Biology, Kolkata, West Bengal, 700032, India
| |
Collapse
|
38
|
Holzer G, Markov GV, Laudet V. Evolution of Nuclear Receptors and Ligand Signaling. Curr Top Dev Biol 2017; 125:1-38. [DOI: 10.1016/bs.ctdb.2017.02.003] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
39
|
Hu L, Chan KCC. Extracting Coevolutionary Features from Protein Sequences for Predicting Protein-Protein Interactions. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2017; 14:155-166. [PMID: 26812730 DOI: 10.1109/tcbb.2016.2520923] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Knowing the ways proteins interact with each other are crucial to our understanding of the functional mechanisms of proteins. It is for this reason that different approaches have been developed in attempts to predict protein-protein interactions (PPIs) computationally. Among them, the sequence-based approaches are preferred to the others as they do not require any information about protein properties to perform their tasks. Instead, most sequence-based approaches make use of feature extraction methods to extract features directly from protein sequences so that for each protein sequence, we can construct a feature vector. The feature vectors of every pair of proteins are then concatenated to form two classes of interacting and non-interacting proteins. The prediction of whether or not two proteins interact with each other is then formulated as a classification problem. How accurate PPI predictions can be made therefore depends on how good the features are that can be extracted from the protein sequences to allow interacting or non-interacting to be best distinguished. To do so, instead of extracting such features from individual protein sequences independently of the other protein in the same pair, we propose to jointly consider features from both sequences in a protein pair during the feature extraction process through using a novel coevolutionary feature extraction approach called CoFex. Coevolutionary features extracted by CoFex refer to the covariations found at coevolving positions. Based on the presence and absence of these coevolutionary features in the sequences of two proteins, feature vectors can be composed for pairs of proteins rather than individual proteins. The experiment results show that CoFex is a promising feature extraction approach and can improve the performance of PPI prediction.
Collapse
|
40
|
Chowdhury A, Rakshit P, Konar A. Protein-protein interaction network prediction using stochastic learning automata induced differential evolution. Appl Soft Comput 2016. [DOI: 10.1016/j.asoc.2016.08.053] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
41
|
Yang Z, Liu L, Fang H, Li P, Xu S, Cao W, Xu C, Huang J, Zhou Y. Origin of the plant Tm-1-like gene via two independent horizontal transfer events and one gene fusion event. Sci Rep 2016; 6:33691. [PMID: 27647002 PMCID: PMC5028733 DOI: 10.1038/srep33691] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2016] [Accepted: 08/31/2016] [Indexed: 01/07/2023] Open
Abstract
The Tomato mosaic virus (ToMV) resistance gene Tm-1 encodes a direct inhibitor of ToMV RNA replication to protect tomato from infection. The plant Tm-1-like (Tm-1L) protein is predicted to contain an uncharacterized N-terminal UPF0261 domain and a C-terminal TIM-barrel signal transduction (TBST) domain. Homologous searches revealed that proteins containing both of these two domains are mainly present in charophyte green algae and land plants but absent from glaucophytes, red algae and chlorophyte green algae. Although Tm-1 homologs are widely present in bacteria, archaea and fungi, UPF0261- and TBST-domain-containing proteins are generally encoded by different genes in these linages. A co-evolution analysis also suggested a putative interaction between UPF0261- and TBST-domain-containing proteins. Phylogenetic analyses based on homologs of these two domains revealed that plants have acquired UPF0261- and TBST-domain-encoding genes through two independent horizontal gene transfer (HGT) events before the origin of land plants from charophytes. Subsequently, gene fusion occurred between these two horizontally acquired genes and resulted in the origin of the Tm-1L gene in streptophytes. Our results demonstrate a novel evolutionary mechanism through which the recipient organism may acquire genes with functional interaction through two different HGT events and further fuse them into one functional gene.
Collapse
Affiliation(s)
- Zefeng Yang
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Key Laboratory of Plant Functional Genomics of the Ministry of Education, Yangzhou University, Yangzhou, 225009, China
| | - Li Liu
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Key Laboratory of Plant Functional Genomics of the Ministry of Education, Yangzhou University, Yangzhou, 225009, China
| | - Huimin Fang
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Key Laboratory of Plant Functional Genomics of the Ministry of Education, Yangzhou University, Yangzhou, 225009, China
| | - Pengcheng Li
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Key Laboratory of Plant Functional Genomics of the Ministry of Education, Yangzhou University, Yangzhou, 225009, China
| | - Shuhui Xu
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Key Laboratory of Plant Functional Genomics of the Ministry of Education, Yangzhou University, Yangzhou, 225009, China
| | - Wei Cao
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Key Laboratory of Plant Functional Genomics of the Ministry of Education, Yangzhou University, Yangzhou, 225009, China
| | - Chenwu Xu
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Key Laboratory of Plant Functional Genomics of the Ministry of Education, Yangzhou University, Yangzhou, 225009, China
| | - Jinling Huang
- Department of Biology, East Carolina University, Greenville, NC, 27858, USA
| | - Yong Zhou
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Key Laboratory of Plant Functional Genomics of the Ministry of Education, Yangzhou University, Yangzhou, 225009, China
| |
Collapse
|
42
|
Zhong J, Jiang L, Cheng C, Huang Z, Zhang H, Liu H, He J, Cao F, Peng J, Jiang Y, Sun X. Altered expression of long non-coding RNA and mRNA in mouse cortex after traumatic brain injury. Brain Res 2016; 1646:589-600. [PMID: 27380725 DOI: 10.1016/j.brainres.2016.07.002] [Citation(s) in RCA: 51] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2016] [Revised: 06/27/2016] [Accepted: 07/01/2016] [Indexed: 01/01/2023]
Abstract
BACKGROUND AND OBJECTIVE The present study aims to detect the altered lncRNA expression in the mouse cortex after traumatic brain injury (TBI). We also simultaneously detected the altered mRNA profile to further analyze the possible function of lncRNA. METHOD C57BL/6 mice (n=18) were used to construct a controlled cortical impact model. At 24h post-TBI, the cortex around injury site was collected and the total RNA was extracted to construct the cDNA library. RNA sequencing (RNA-seq) was carried out followed by RT-PCR for confirmation. Bioinformatic analysis (including GO analysis, KEGG pathway and co-expression analysis) also were performed. RESULTS A total of 64,530 transcripts were detected in the current sequencing study, in which 27,457 transcripts were identified as mRNA and 37,073 transcripts as lncRNA. A total of 1580 mRNAs (1430 up-regulated and 150 down-regulated) and 823 lncRNAs (667 up-regulated and 156 down-regulated) were significantly changed according to the criteria ( (|)log2((fold change))|>1 and P<0.05). These altered mRNAs were mainly related to inflammatory and immunological activity, metabolism, neuronal and vascular network. The expression of single lncRNA may be related with several mRNAs, and so was the mRNA. Also, a total of 360 new mRNAs and 8041 new lncRNAs were identified. The good reproducibility and reliability of RNA-seq were confirmed by RT-PCR. CONCLUSION Numerous lncRNAs and mRNAs were significantly altered in mouse cortex around the injury site 24h after TBI. Our present data may provide a promising approach for further study about TBI.
Collapse
Affiliation(s)
- Jianjun Zhong
- Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing 400016, China
| | - Li Jiang
- Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing 400016, China
| | - Chongjie Cheng
- Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing 400016, China
| | - Zhijian Huang
- Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing 400016, China
| | - Hongrong Zhang
- Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing 400016, China
| | - Han Liu
- Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing 400016, China
| | - Junchi He
- Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing 400016, China
| | - Fang Cao
- Department of Cerebrovascular, the First Affiliated Hospital of Zunyi Medical College, Zunyi, Guizhou 653000, China
| | - Jianhua Peng
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, Sichaun 646000, China
| | - Yong Jiang
- Department of Neurosurgery, the Affiliated Hospital of Southwest Medical University, Luzhou, Sichaun 646000, China
| | - Xiaochuan Sun
- Department of Neurosurgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing 400016, China
| |
Collapse
|
43
|
Chowdhury A, Rakshit P, Konar A. Prediction of protein-protein interaction network using a multi-objective optimization approach. J Bioinform Comput Biol 2016; 14:1650008. [PMID: 26846814 DOI: 10.1142/s0219720016500086] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Protein-Protein Interactions (PPIs) are very important as they coordinate almost all cellular processes. This paper attempts to formulate PPI prediction problem in a multi-objective optimization framework. The scoring functions for the trial solution deal with simultaneous maximization of functional similarity, strength of the domain interaction profiles, and the number of common neighbors of the proteins predicted to be interacting. The above optimization problem is solved using the proposed Firefly Algorithm with Nondominated Sorting. Experiments undertaken reveal that the proposed PPI prediction technique outperforms existing methods, including gene ontology-based Relative Specific Similarity, multi-domain-based Domain Cohesion Coupling method, domain-based Random Decision Forest method, Bagging with REP Tree, and evolutionary/swarm algorithm-based approaches, with respect to sensitivity, specificity, and F1 score.
Collapse
Affiliation(s)
- Archana Chowdhury
- 1 Artificial Intelligence Laboratory, Department of Electronics and Telecommunication Engineering, Jadavpur University, Kolkata, India
| | - Pratyusha Rakshit
- 1 Artificial Intelligence Laboratory, Department of Electronics and Telecommunication Engineering, Jadavpur University, Kolkata, India
| | - Amit Konar
- 1 Artificial Intelligence Laboratory, Department of Electronics and Telecommunication Engineering, Jadavpur University, Kolkata, India
| |
Collapse
|
44
|
Abstract
It is long known that peptide neurotoxins derived from a diversity of venomous animals evolve by positive selection following gene duplication, yet a force that drives their adaptive evolution remains a mystery. By using maximum-likelihood models of codon substitution, we analyzed molecular adaptation in scorpion sodium channel toxins from a specific species and found ten positively selected sites, six of which are located at the core-domain of scorpion α-toxins, a region known to interact with two adjacent loops in the voltage-sensor domain (DIV) of sodium channels, as validated by our newly constructed computational model of toxin-channel complex. Despite the lack of positive selection signals in these two loops, they accumulated extensive sequence variations by relaxed purifying selection in prey and predators of scorpions. The evolutionary variability in the toxin-bound regions of sodium channels indicates that accelerated substitutions in the multigene family of scorpion toxins is a consequence of dealing with the target diversity. This work presents an example of atypical co-evolution between animal toxins and their molecular targets, in which toxins suffered from more prominent selective pressure from the channels of their competitors. Our discovery helps explain the evolutionary rationality of gene duplication of toxins in a specific venomous species.
Collapse
|
45
|
Gao H, Yu X, Dou Y, Wang J. New Measurement for Correlation of Co-evolution Relationship of Subsequences in Protein. Interdiscip Sci 2015; 7:364-72. [PMID: 26396121 DOI: 10.1007/s12539-015-0024-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2014] [Revised: 04/08/2014] [Accepted: 04/16/2014] [Indexed: 11/26/2022]
Abstract
Many computational tools have been developed to measure the protein residues co-evolution. Most of them only focus on co-evolution for pairwise residues in a protein sequence. However, number of residues participate in co-evolution might be multiple. And some co-evolved residues are clustered in several distinct regions in primary structure. Therefore, the co-evolution among the adjacent residues and the correlation between the distinct regions offer insights into function and evolution of the protein and residues. Subsequence is used to represent the adjacent multiple residues in one distinct region. In the paper, co-evolution relationship in each subsequence is represented by mutual information matrix (MIM). Then, Pearson's correlation coefficient: R value is developed to measure the similarity correlation of two MIMs. MSAs from Catalytic Data Base (Catalytic Site Atlas, CSA) are used for testing. R value characterizes a specific class of residues. In contrast to individual pairwise co-evolved residues, adjacent residues without high individual MI values are found since the co-evolved relationship among them is similar to that among another set of adjacent residues. These subsequences possess some flexibility in the composition of side chains, such as the catalyzed environment.
Collapse
Affiliation(s)
- Hongyun Gao
- School of Mathematical Sciences, Dalian University of Technology, Dalian, 116024, China
- Information and Engineering College, Dalian University, Dalian, 116622, China
| | - Xiaoqing Yu
- College of Sciences, Shanghai Institute of Technology, Shanghai, 201418, China
| | - Yongchao Dou
- Center for Plant Science and Innovation, School of Biological Sciences, University of Nebraska, Lincoln, NE, 68588, USA
| | - Jun Wang
- Department of Mathematics, Shanghai Normal University, Shanghai, 200234, China.
| |
Collapse
|
46
|
Heavner ME, Qiu WG, Cheng HP. Phylogenetic Co-Occurrence of ExoR, ExoS, and ChvI, Components of the RSI Bacterial Invasion Switch, Suggests a Key Adaptive Mechanism Regulating the Transition between Free-Living and Host-Invading Phases in Rhizobiales. PLoS One 2015; 10:e0135655. [PMID: 26309130 PMCID: PMC4550343 DOI: 10.1371/journal.pone.0135655] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2015] [Accepted: 07/23/2015] [Indexed: 11/18/2022] Open
Abstract
Both bacterial symbionts and pathogens rely on their host-sensing mechanisms to activate the biosynthetic pathways necessary for their invasion into host cells. The Gram-negative bacterium Sinorhizobium meliloti relies on its RSI (ExoR-ExoS-ChvI) Invasion Switch to turn on the production of succinoglycan, an exopolysaccharide required for its host invasion. Recent whole-genome sequencing efforts have uncovered putative components of RSI-like invasion switches in many other symbiotic and pathogenic bacteria. To explore the possibility of the existence of a common invasion switch, we have conducted a phylogenomic survey of orthologous ExoR, ExoS, and ChvI tripartite sets in more than ninety proteobacterial genomes. Our analyses suggest that functional orthologs of the RSI invasion switch co-exist in Rhizobiales, an order characterized by numerous invasive species, but not in the order’s close relatives. Phylogenomic analyses and reconstruction of orthologous sets of the three proteins in Alphaproteobacteria confirm Rhizobiales-specific gene synteny and congruent RSI evolutionary histories. Evolutionary analyses further revealed site-specific substitutions correlated specifically to either animal-bacteria or plant-bacteria associations. Lineage restricted conservation of any one specialized gene is in itself an indication of species adaptation. However, the orthologous phylogenetic co-occurrence of all interacting partners within this single signaling pathway strongly suggests that the development of the RSI switch was a key adaptive mechanism. The RSI invasion switch, originally found in S. meliloti, is a characteristic of the Rhizobiales, and potentially a conserved crucial activation step that may be targeted to control host invasion by pathogenic bacterial species.
Collapse
Affiliation(s)
- Mary Ellen Heavner
- Biochemistry Program, The Graduate Center, City University of New York, New York, New York, United States of America
| | - Wei-Gang Qiu
- Biological Sciences Department, Hunter College, City University of New York, New York, New York, United States of America
| | - Hai-Ping Cheng
- Biochemistry Program, The Graduate Center, City University of New York, New York, New York, United States of America
- Biological Sciences Department, Lehman College, City University of New York, Bronx, New York, United States of America
- * E-mail:
| |
Collapse
|
47
|
Identification of Protein–Protein Interactions by Detecting Correlated Mutation at the Interface. J Chem Inf Model 2015; 55:2042-9. [DOI: 10.1021/acs.jcim.5b00320] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
|
48
|
Wang YF, Yan JJ, Tseng YC, Chen RD, Hwang PP. Molecular Physiology of an Extra-renal Cl(-) Uptake Mechanism for Body Fluid Cl(-) Homeostasis. Int J Biol Sci 2015; 11:1190-203. [PMID: 26327813 PMCID: PMC4551755 DOI: 10.7150/ijbs.11737] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2015] [Accepted: 07/25/2015] [Indexed: 01/07/2023] Open
Abstract
The development of an ion regulatory mechanism for body fluid homeostasis was an important trait for vertebrates during the evolution from aquatic to terrestrial life. The homeostatic mechanism of Cl- in aquatic fish appears to be similar to that of terrestrial vertebrates; however, the mechanism in non-mammalian vertebrates is poorly understood. Unlike in mammals, in which the kidney plays a central role, in most fish species, the gill is responsible for the maintenance of Cl- homeostasis via Cl- transport uptake mechanisms. Previous studies in zebrafish identified Na+-Cl- cotransporter (NCC) 2b-expressing cells in the gills and skin as the major ionocytes responsible for Cl- uptake, similar to distal convoluted tubular cells in mammalian kidney. However, the mechanism by which basolateral ions exit from NCC cells is still unclear. Of the in situ hybridization signals of twelve members of the clc Cl- channel family, only that of clc-2c exhibited an ionocyte pattern in the gill and embryonic skin. Double in situ hybridization/immunocytochemistry confirmed colocalization of apical NCC2b with basolateral CLC-2c. Acclimation to a low Cl- environment increased mRNA expression of both clc-2c and ncc2b, and also the protein expression of CLC-2c in embryos and adult gills. Loss-of-function of clc-2c resulted in a significant decrease in whole body Cl- content in zebrafish embryos, a phenotype similar to that of ncc2b mutants; this finding suggests a role for CLC-2c in Cl- uptake. Translational knockdown of clc-2c stimulated ncc2b mRNA expression and vice versa, revealing cooperation between these two transporters in the context of zebrafish Cl- homeostasis. Further comparative genomic and phylogenetic analyses revealed that zebrafish CLC-2c is a fish-specific isoform that diverged from a kidney-predominant homologue, in the same manner as NCC2b and its counterparts (NCCs). Several lines of molecular and cellular physiological evidences demonstrated the cofunctional role of apical NCC2b and basolateral CLC-2c in the gill/skin Cl- uptake pathway. Taking the phylogenetic evidence into consideration, fish-specific NCC2b and CLC-2c may have coevolved to perform extra-renal Cl- uptake during the evolution of vertebrates in an aquatic environment.
Collapse
Affiliation(s)
- Yi-Fang Wang
- 1. Institute of Fishery Science, National Taiwan University, Taipei city, Taiwan ; 2. Institute of Cellular and Organismic Biology, Academia Sinica, Taipei city, Taiwan
| | - Jia-Jiun Yan
- 2. Institute of Cellular and Organismic Biology, Academia Sinica, Taipei city, Taiwan
| | - Yung-Che Tseng
- 3. Department of Life Science, National Taiwan Normal University, Taipei city, Taiwan
| | - Ruo-Dong Chen
- 2. Institute of Cellular and Organismic Biology, Academia Sinica, Taipei city, Taiwan
| | - Pung-Pung Hwang
- 1. Institute of Fishery Science, National Taiwan University, Taipei city, Taiwan ; 2. Institute of Cellular and Organismic Biology, Academia Sinica, Taipei city, Taiwan
| |
Collapse
|
49
|
Ramilowski JA, Goldberg T, Harshbarger J, Kloppman E, Lizio M, Satagopam VP, Itoh M, Kawaji H, Carninci P, Rost B, Forrest ARR. A draft network of ligand-receptor-mediated multicellular signalling in human. Nat Commun 2015; 6:7866. [PMID: 26198319 PMCID: PMC4525178 DOI: 10.1038/ncomms8866] [Citation(s) in RCA: 522] [Impact Index Per Article: 58.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2014] [Accepted: 06/19/2015] [Indexed: 02/07/2023] Open
Abstract
Cell-to-cell communication across multiple cell types and tissues strictly governs proper functioning of metazoans and extensively relies on interactions between secreted ligands and cell-surface receptors. Herein, we present the first large-scale map of cell-to-cell communication between 144 human primary cell types. We reveal that most cells express tens to hundreds of ligands and receptors to create a highly connected signalling network through multiple ligand-receptor paths. We also observe extensive autocrine signalling with approximately two-thirds of partners possibly interacting on the same cell type. We find that plasma membrane and secreted proteins have the highest cell-type specificity, they are evolutionarily younger than intracellular proteins, and that most receptors had evolved before their ligands. We provide an online tool to interactively query and visualize our networks and demonstrate how this tool can reveal novel cell-to-cell interactions with the prediction that mast cells signal to monoblastic lineages via the CSF1-CSF1R interacting pair.
Collapse
Affiliation(s)
- Jordan A. Ramilowski
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045 Japan
| | - Tatyana Goldberg
- Department for Bioinformatics and Computational Biology-I12, Technische Universität München (TUM), Boltzmannstrasse 3, 85748 Garching, Germany
- TUM Graduate School, Center of Doctoral Studies in Informatics and its Applications (CeDoSIA), Boltzmannstrasse 11, 85748 Garching, Germany
| | - Jayson Harshbarger
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045 Japan
| | - Edda Kloppman
- Department for Bioinformatics and Computational Biology-I12, Technische Universität München (TUM), Boltzmannstrasse 3, 85748 Garching, Germany
| | - Marina Lizio
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045 Japan
| | - Venkata P. Satagopam
- Luxembourg Centre for Systems Biomedicine, Campus Belval, 7 Avenue des Hauts Fourneaux, L-4362 Belval, Luxembourg
| | - Masayoshi Itoh
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045 Japan
- RIKEN Preventive Medicine and Diagnosis Innovation Program, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Hideya Kawaji
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045 Japan
- RIKEN Preventive Medicine and Diagnosis Innovation Program, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Piero Carninci
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045 Japan
| | - Burkhard Rost
- Department for Bioinformatics and Computational Biology-I12, Technische Universität München (TUM), Boltzmannstrasse 3, 85748 Garching, Germany
- TUM Graduate School, Center of Doctoral Studies in Informatics and its Applications (CeDoSIA), Boltzmannstrasse 11, 85748 Garching, Germany
| | - Alistair R. R. Forrest
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045 Japan
- Harry Perkins Institute of Medical Research, QEII Medical Centre and Centre for Medical Research, the University of Western Australia, PO Box 7214, 6 Verdun Street, Nedlands, Perth, Western Australia 6008, Australia
| |
Collapse
|
50
|
Cao C, Zhao J, Doughty EK, Migliorini M, Strickland DK, Kann MG, Zhang L. Mac-1 Regulates IL-13 Activity in Macrophages by Directly Interacting with IL-13Rα1. J Biol Chem 2015; 290:21642-51. [PMID: 26160172 DOI: 10.1074/jbc.m115.645796] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2015] [Indexed: 11/06/2022] Open
Abstract
Mac-1 exhibits a unique inhibitory activity toward IL-13-induced JAK/STAT activation and thereby regulates macrophage to foam cell transformation. However, the underlying molecular mechanism is unknown. In this study, we report the identification of IL-13Rα1, a component of the IL-13 receptor (IL-13R), as a novel ligand of integrin Mac-1, using a co-evolution-based algorithm. Biochemical analyses demonstrated that recombinant IL-13Rα1 binds Mac-1 in a purified system and supports Mac-1-mediated cell adhesion. Co-immunoprecipitation experiments revealed that endogenous Mac-1 forms a complex with IL-13Rα1 in solution, and confocal fluorescence microscopy demonstrated that these two receptors co-localize with each other on the surface of macrophages. Moreover, we found that genetic inactivation of Mac-1 promotes IL-13-induced JAK/STAT activation in macrophages, resulting in enhanced polarization along the alternative activation pathway. Importantly, we observed that Mac-1(-/-) macrophages exhibit increased expression of foam cell differentiation markers including 15-lipoxygenase and lectin-type oxidized LDL receptor-1 both in vitro and in vivo. Indeed, we found that Mac-1(-/-)LDLR(-/-) mice develop significantly more foam cells than control LDLR(-/-) mice, using an in vivo model of foam cell formation. Together, our data establish for the first time a molecular mechanism by which Mac-1 regulates the signaling activity of IL-13 in macrophages. This newly identified IL-13Rα1/Mac-1-dependent pathway may offer novel targets for therapeutic intervention in the future.
Collapse
Affiliation(s)
| | | | - Emily K Doughty
- the Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, Maryland 21250
| | - Mary Migliorini
- Surgery, Center for Vascular and Inflammatory Diseases, the University of Maryland, School of Medicine, Baltimore, Maryland 21201 and
| | - Dudley K Strickland
- Surgery, Center for Vascular and Inflammatory Diseases, the University of Maryland, School of Medicine, Baltimore, Maryland 21201 and
| | - Maricel G Kann
- the Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, Maryland 21250
| | - Li Zhang
- From the Departments of Physiology and
| |
Collapse
|