1
|
Baker M, Zhang X, Maciel-Guerra A, Babaarslan K, Dong Y, Wang W, Hu Y, Renney D, Liu L, Li H, Hossain M, Heeb S, Tong Z, Pearcy N, Zhang M, Geng Y, Zhao L, Hao Z, Senin N, Chen J, Peng Z, Li F, Dottorini T. Convergence of resistance and evolutionary responses in Escherichia coli and Salmonella enterica co-inhabiting chicken farms in China. Nat Commun 2024; 15:206. [PMID: 38182559 PMCID: PMC10770378 DOI: 10.1038/s41467-023-44272-1] [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: 04/04/2023] [Accepted: 12/06/2023] [Indexed: 01/07/2024] Open
Abstract
Sharing of genetic elements among different pathogens and commensals inhabiting same hosts and environments has significant implications for antimicrobial resistance (AMR), especially in settings with high antimicrobial exposure. We analysed 661 Escherichia coli and Salmonella enterica isolates collected within and across hosts and environments, in 10 Chinese chicken farms over 2.5 years using data-mining methods. Most isolates within same hosts possessed the same clinically relevant AMR-carrying mobile genetic elements (plasmids: 70.6%, transposons: 78%), which also showed recent common evolution. Supervised machine learning classifiers revealed known and novel AMR-associated mutations and genes underlying resistance to 28 antimicrobials, primarily associated with resistance in E. coli and susceptibility in S. enterica. Many were essential and affected same metabolic processes in both species, albeit with varying degrees of phylogenetic penetration. Multi-modal strategies are crucial to investigate the interplay of mobilome, resistance and metabolism in cohabiting bacteria, especially in ecological settings where community-driven resistance selection occurs.
Collapse
Affiliation(s)
- Michelle Baker
- School of Veterinary Medicine and Science, University of Nottingham, College Road, Sutton Bonington, Loughborough, Leicestershire, LE12 5RD, UK
| | - Xibin Zhang
- Shandong New Hope Liuhe Group Co. Ltd. and Qingdao Key Laboratory of Animal Feed Safety, Qingdao, Shandong, 266000, P.R. China
| | - Alexandre Maciel-Guerra
- School of Veterinary Medicine and Science, University of Nottingham, College Road, Sutton Bonington, Loughborough, Leicestershire, LE12 5RD, UK
| | - Kubra Babaarslan
- School of Veterinary Medicine and Science, University of Nottingham, College Road, Sutton Bonington, Loughborough, Leicestershire, LE12 5RD, UK
| | - Yinping Dong
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, 100021, P. R. China
| | - Wei Wang
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, 100021, P. R. China
| | - Yujie Hu
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, 100021, P. R. China
| | - David Renney
- Nimrod Veterinary Products Limited, 2, Wychwood Court, Cotswold Business Village, Moreton-in-Marsh, GL56 0JQ, London, UK
| | - Longhai Liu
- Shandong Kaijia Food Co. Ltd, Weifang, P. R. China
| | - Hui Li
- Luoyang Center for Disease Control and Prevention, No. 9, Zhenghe Road, Luolong District, Luoyang City, Henan Province, Luolong, 471000, P. R. China
| | - Maqsud Hossain
- School of Veterinary Medicine and Science, University of Nottingham, College Road, Sutton Bonington, Loughborough, Leicestershire, LE12 5RD, UK
| | - Stephan Heeb
- School of Life Sciences, University of Nottingham, East Drive, Nottingham, Nottinghamshire, NG7 2RD, UK
| | - Zhiqin Tong
- Luoyang Center for Disease Control and Prevention, No. 9, Zhenghe Road, Luolong District, Luoyang City, Henan Province, Luolong, 471000, P. R. China
| | - Nicole Pearcy
- School of Veterinary Medicine and Science, University of Nottingham, College Road, Sutton Bonington, Loughborough, Leicestershire, LE12 5RD, UK
- School of Life Sciences, University of Nottingham, East Drive, Nottingham, Nottinghamshire, NG7 2RD, UK
| | - Meimei Zhang
- Liaoning Provincial Center for Disease Control and Prevention, No. 168, Jinfeng Street, Hunnan District, Shenyang City, Liaoning Province, 110072, P. R. China
| | - Yingzhi Geng
- Liaoning Provincial Center for Disease Control and Prevention, No. 168, Jinfeng Street, Hunnan District, Shenyang City, Liaoning Province, 110072, P. R. China
| | - Li Zhao
- Agricultural Biopharmaceutical Laboratory, College of Chemistry and Pharmaceutical Sciences, Qingdao Agricultural University, No. 700 Changcheng Road, Chengyang District, Qingdao City, Shandong Province, 266109, P. R. China
| | - Zhihui Hao
- Chinese Veterinary Medicine Innovation Center, College of Veterinary Medicine, China Agricultural University, Haidian District, Beijing City, 100193, P. R. China
| | - Nicola Senin
- Department of Engineering, University of Perugia, Perugia, I06125, Italy
| | - Junshi Chen
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, 100021, P. R. China
| | - Zixin Peng
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, 100021, P. R. China.
| | - Fengqin Li
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, 100021, P. R. China.
| | - Tania Dottorini
- School of Veterinary Medicine and Science, University of Nottingham, College Road, Sutton Bonington, Loughborough, Leicestershire, LE12 5RD, UK.
- Centre for Smart Food Research, Nottingham Ningbo China Beacons of Excellence Research and Innovation Institute, University of Nottingham Ningbo China, Ningbo, 315100, P. R. China.
| |
Collapse
|
2
|
Tang QY, Ren W, Wang J, Kaneko K. The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database. Mol Biol Evol 2022; 39:6701686. [PMID: 36108094 PMCID: PMC9550990 DOI: 10.1093/molbev/msac197] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.
Collapse
Affiliation(s)
| | - Weitong Ren
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Jun Wang
- School of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, People’s Republic of China
| | | |
Collapse
|
3
|
Jakubec D, Skoda P, Krivak R, Novotny M, Hoksza D. PrankWeb 3: accelerated ligand-binding site predictions for experimental and modelled protein structures. Nucleic Acids Res 2022; 50:W593-W597. [PMID: 35609995 DOI: 10.1093/nar/gkac389] [Citation(s) in RCA: 48] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Revised: 04/15/2022] [Accepted: 05/06/2022] [Indexed: 11/13/2022] Open
Abstract
Knowledge of protein-ligand binding sites (LBSs) enables research ranging from protein function annotation to structure-based drug design. To this end, we have previously developed a stand-alone tool, P2Rank, and the web server PrankWeb (https://prankweb.cz/) for fast and accurate LBS prediction. Here, we present significant enhancements to PrankWeb. First, a new, more accurate evolutionary conservation estimation pipeline based on the UniRef50 sequence database and the HMMER3 package is introduced. Second, PrankWeb now allows users to enter UniProt ID to carry out LBS predictions in situations where no experimental structure is available by utilizing the AlphaFold model database. Additionally, a range of minor improvements has been implemented. These include the ability to deploy PrankWeb and P2Rank as Docker containers, support for the mmCIF file format, improved public REST API access, or the ability to batch download the LBS predictions for the whole PDB archive and parts of the AlphaFold database.
Collapse
Affiliation(s)
- David Jakubec
- Department of Software Engineering, Faculty of Mathematics and Physics, Charles University, Czech Republic
| | - Petr Skoda
- Department of Software Engineering, Faculty of Mathematics and Physics, Charles University, Czech Republic
| | - Radoslav Krivak
- Department of Software Engineering, Faculty of Mathematics and Physics, Charles University, Czech Republic
| | - Marian Novotny
- Department of Cell Biology, Faculty of Science, Charles University, Czech Republic
| | - David Hoksza
- Department of Software Engineering, Faculty of Mathematics and Physics, Charles University, Czech Republic
| |
Collapse
|
4
|
Kimura M, Imai K, Morinaka Y, Hosono-Sakuma Y, Horton P, Imamoto N. Distinct mutations in importin-β family nucleocytoplasmic transport receptors transportin-SR and importin-13 affect specific cargo binding. Sci Rep 2021; 11:15649. [PMID: 34341383 PMCID: PMC8329185 DOI: 10.1038/s41598-021-94948-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2020] [Accepted: 07/20/2021] [Indexed: 01/25/2023] Open
Abstract
Importin-(Imp)β family nucleocytoplasmic transport receptors (NTRs) are supposed to bind to their cargoes through interaction between a confined interface on an NTR and a nuclear localization or export signal (NLS/NES) on a cargo. Although consensus NLS/NES sequence motifs have been defined for cargoes of some NTRs, many experimentally identified cargoes of those NTRs lack those motifs, and consensus NLSs/NESs have been reported for only a few NTRs. Crystal structures of NTR-cargo complexes have exemplified 3D structure-dependent binding of cargoes lacking a consensus NLS/NES to different sites on an NTR. Since only a limited number of NTR-cargo interactions have been studied, whether most cargoes lacking a consensus NLS/NES bind to the same confined interface or to various sites on an NTR is still unclear. Addressing this issue, we generated four mutants of transportin-(Trn)SR, of which many cargoes lack a consensus NLS, and eight mutants of Imp13, where no consensus NLS has been defined, and we analyzed their binding to as many as 40 cargo candidates that we previously identified by a nuclear import reaction-based method. The cargoes bind differently to the NTR mutants, suggesting that positions on an NTR contribute differently to the binding of respective cargoes.
Collapse
Affiliation(s)
- Makoto Kimura
- Cellular Dynamics Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan.
| | - Kenichiro Imai
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, Japan.
- Molecular Profiling Research Center for Drug Discovery, National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, Japan.
| | - Yuriko Morinaka
- Cellular Dynamics Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
| | - Yoshiko Hosono-Sakuma
- Cellular Dynamics Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
| | - Paul Horton
- Department of Computer Science and Information Engineering, National Cheng Kung University, Tainan City, Taiwan
| | - Naoko Imamoto
- Cellular Dynamics Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan.
| |
Collapse
|
5
|
Arsiccio A, Shea JE. Pressure Unfolding of Proteins: New Insights into the Role of Bound Water. J Phys Chem B 2021; 125:8431-8442. [PMID: 34310136 DOI: 10.1021/acs.jpcb.1c04398] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
High pressures can be detrimental for protein stability, resulting in unfolding and loss of function. This phenomenon occurs because the unfolding transition is accompanied by a decrease in volume, which is typically attributed to the elimination of cavities that are present within the native state as a result of packing defects. We present a novel computational approach that enables the study of pressure unfolding in atomistically detailed protein models in implicit solvent. We include the effect of pressure using a transfer free energy term that allows us to decouple the effect of protein residues and bound water molecules on the volume change upon unfolding. We discuss molecular dynamics simulations results using this protocol for two model proteins, Trp-cage and staphylococcal nuclease (SNase). We find that the volume reduction of bound water is the key energetic term that drives protein denaturation under the effect of pressure, for both Trp-cage and SNase. However, we note differences in unfolding mechanisms between the smaller Trp-cage and the larger SNase protein. Indeed, the unfolding of SNase, but not Trp-cage, is seen to be further accompanied by a reduction in the volume of internal cavities. Our results indicate that, for small peptides, like Trp-cage, pressure denaturation is driven by the increase in solvent accessibility upon unfolding, and the subsequent increase in the number of bound water molecules. For larger proteins, like SNase, the cavities within the native fold act as weak spots, determining the overall resistance to pressure denaturation. Our simulations display a striking agreement with the pressure-unfolding profile experimentally obtained for SNase and represent a promising approach for a computationally efficient and accurate exploration of pressure-induced denaturation of proteins.
Collapse
Affiliation(s)
- Andrea Arsiccio
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States.,Department of Physics, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
6
|
Vymětal J, Jakubec D, Galgonek J, Vondrášek J. Amino Acid Interactions (INTAA) web server v2.0: a single service for computation of energetics and conservation in biomolecular 3D structures. Nucleic Acids Res 2021; 49:W15-W20. [PMID: 34019656 PMCID: PMC8262704 DOI: 10.1093/nar/gkab377] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2021] [Revised: 04/10/2021] [Accepted: 04/27/2021] [Indexed: 12/29/2022] Open
Abstract
Interactions among amino acid residues are the principal contributor to the stability of the three-dimensional structure of a protein. The Amino Acid Interactions (INTAA) web server (https://bioinfo.uochb.cas.cz/INTAA/) has established itself as a unique computational resource, which enables users to calculate the contribution of individual residues in a biomolecular structure to its total energy using a molecular mechanical scoring function. In this update, we describe major additions to the web server which help solidify its position as a robust, comprehensive resource for biomolecular structure analysis. Importantly, a new continuum solvation model was introduced, allowing more accurate representation of electrostatic interactions in aqueous media. In addition, a low-overhead pipeline for the estimation of evolutionary conservation in protein chains has been added. New visualization options were introduced as well, allowing users to easily switch between and interrelate the energetic and evolutionary views of the investigated structures.
Collapse
Affiliation(s)
- Jiří Vymětal
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Praha 6, 160 00, Czech Republic
| | - David Jakubec
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Praha 6, 160 00, Czech Republic.,Department of Software Engineering, Faculty of Mathematics and Physics, Charles University, Praha 1, 118 00, Czech Republic
| | - Jakub Galgonek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Praha 6, 160 00, Czech Republic
| | - Jiří Vondrášek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Praha 6, 160 00, Czech Republic
| |
Collapse
|
7
|
Wang W, Baker M, Hu Y, Xu J, Yang D, Maciel-Guerra A, Xue N, Li H, Yan S, Li M, Bai Y, Dong Y, Peng Z, Ma J, Li F, Dottorini T. Whole-Genome Sequencing and Machine Learning Analysis of Staphylococcus aureus from Multiple Heterogeneous Sources in China Reveals Common Genetic Traits of Antimicrobial Resistance. mSystems 2021; 6:e0118520. [PMID: 34100643 PMCID: PMC8579812 DOI: 10.1128/msystems.01185-20] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Accepted: 05/10/2021] [Indexed: 12/27/2022] Open
Abstract
Staphylococcus aureus is a worldwide leading cause of numerous diseases ranging from food-poisoning to lethal infections. Methicillin-resistant S. aureus (MRSA) has been found capable of acquiring resistance to most antimicrobials. MRSA is ubiquitous and diverse even in terms of antimicrobial resistance (AMR) profiles, posing a challenge for treatment. Here, we present a comprehensive study of S. aureus in China, addressing epidemiology, phylogenetic reconstruction, genomic characterization, and identification of AMR profiles. The study analyzes 673 S. aureus isolates from food as well as from hospitalized and healthy individuals. The isolates have been collected over a 9-year period, between 2010 and 2018, from 27 provinces across China. By whole-genome sequencing, Bayesian divergence analysis, and supervised machine learning, we reconstructed the phylogeny of the isolates and compared them to references from other countries. We identified 72 sequence types (STs), of which, 29 were novel. We found 81 MRSA lineages by multilocus sequence type (MLST), spa, staphylococcal cassette chromosome mec element (SCCmec), and Panton-Valentine leukocidin (PVL) typing. In addition, novel variants of SCCmec type IV hosting extra metal and antimicrobial resistance genes, as well as a new SCCmec type, were found. New Bayesian dating of the split times of major clades showed that ST9, ST59, and ST239 in China and European countries fell in different branches, whereas this pattern was not observed for the ST398 clone. On the contrary, the clonal transmission of ST398 was more intermixed in regard to geographic origin. Finally, we identified genetic determinants of resistance to 10 antimicrobials, discriminating drug-resistant bacteria from susceptible strains in the cohort. Our results reveal the emergence of Chinese MRSA lineages enriched of AMR determinants that share similar genetic traits of antimicrobial resistance across human and food, hinting at a complex scenario of evolving transmission routes. IMPORTANCE Little information is available on the epidemiology and characterization of Staphylococcus aureus in China. The role of food is a cause of major concern: staphylococcal foodborne diseases affect thousands every year, and the presence of resistant Staphylococcus strains on raw retail meat products is well documented. We studied a large heterogeneous data set of S. aureus isolates from many provinces of China, isolated from food as well as from individuals. Our large whole-genome collection represents a unique catalogue that can be easily meta-analyzed and integrated with further studies and adds to the library of S. aureus sequences in the public domain in a currently underrepresented geographical region. The new Bayesian dating of the split times of major drug-resistant enriched clones is relevant in showing that Chinese and European methicillin-resistant S. aureus (MRSA) have evolved differently. Our machine learning approach, across a large number of antibiotics, shows novel determinants underlying resistance and reveals frequent resistant traits in specific clonal complexes, highlighting the importance of particular clonal complexes in China. Our findings substantially expand what is known of the evolution and genetic determinants of resistance in food-associated S. aureus in China and add crucial information for whole-genome sequencing (WGS)-based surveillance of S. aureus.
Collapse
Affiliation(s)
- Wei Wang
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Michelle Baker
- School of Veterinary Medicine and Science, University of Nottingham, Sutton Bonington, Leicestershire, United Kingdom
| | - Yue Hu
- School of Veterinary Medicine and Science, University of Nottingham, Sutton Bonington, Leicestershire, United Kingdom
| | - Jin Xu
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Dajin Yang
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | | | - Ning Xue
- School of Veterinary Medicine and Science, University of Nottingham, Sutton Bonington, Leicestershire, United Kingdom
| | - Hui Li
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Shaofei Yan
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Menghan Li
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Yao Bai
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Yinping Dong
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Zixin Peng
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Jinjing Ma
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
- School of Chemistry and Chemical Engineering, Anqing Normal University, Anqing, Anhui, China
| | - Fengqin Li
- NHC Key Laboratory of Food Safety Risk Assessment, China National Center for Food Safety Risk Assessment, Beijing, China
| | - Tania Dottorini
- School of Veterinary Medicine and Science, University of Nottingham, Sutton Bonington, Leicestershire, United Kingdom
| |
Collapse
|
8
|
De Maio N, Alekseyenko AV, Coleman-Smith WJ, Pardi F, Suchard MA, Tamuri AU, Truszkowski J, Goldman N. A phylogenetic approach for weighting genetic sequences. BMC Bioinformatics 2021; 22:285. [PMID: 34049487 PMCID: PMC8164272 DOI: 10.1186/s12859-021-04183-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Accepted: 05/04/2021] [Indexed: 11/22/2022] Open
Abstract
BACKGROUND Many important applications in bioinformatics, including sequence alignment and protein family profiling, employ sequence weighting schemes to mitigate the effects of non-independence of homologous sequences and under- or over-representation of certain taxa in a dataset. These schemes aim to assign high weights to sequences that are 'novel' compared to the others in the same dataset, and low weights to sequences that are over-represented. RESULTS We formalise this principle by rigorously defining the evolutionary 'novelty' of a sequence within an alignment. This results in new sequence weights that we call 'phylogenetic novelty scores'. These scores have various desirable properties, and we showcase their use by considering, as an example application, the inference of character frequencies at an alignment column-important, for example, in protein family profiling. We give computationally efficient algorithms for calculating our scores and, using simulations, show that they are versatile and can improve the accuracy of character frequency estimation compared to existing sequence weighting schemes. CONCLUSIONS Our phylogenetic novelty scores can be useful when an evolutionarily meaningful system for adjusting for uneven taxon sampling is desired. They have numerous possible applications, including estimation of evolutionary conservation scores and sequence logos, identification of targets in conservation biology, and improving and measuring sequence alignment accuracy.
Collapse
Affiliation(s)
- Nicola De Maio
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, UK
| | - Alexander V. Alekseyenko
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, UK
- Present Address: Department of Public Health Sciences, Medical University of South Carolina, Charleston, SC USA
| | - William J. Coleman-Smith
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, UK
| | - Fabio Pardi
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, UK
- Present Address: LIRMM, University of Montpellier, CNRS, Montpellier, France
| | - Marc A. Suchard
- Departments of Biostatistics, Biomathematics and Human Genetics, University of California, Los Angeles, CA USA
| | - Asif U. Tamuri
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, UK
- Present Address: Research IT Services, University College London, London, UK
| | - Jakub Truszkowski
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, UK
- Present Address: RBC Borealis AI, Waterloo, ON Canada
| | - Nick Goldman
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, UK
| |
Collapse
|
9
|
Maiangwa J, Hamdan SH, Mohamad Ali MS, Salleh AB, Zaliha Raja Abd Rahman RN, Shariff FM, Leow TC. Enhancing the stability of Geobacillus zalihae T1 lipase in organic solvents and insights into the structural stability of its variants. J Mol Graph Model 2021; 105:107897. [PMID: 33770705 DOI: 10.1016/j.jmgm.2021.107897] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Revised: 03/04/2021] [Accepted: 03/04/2021] [Indexed: 11/28/2022]
Abstract
Critical to the applications of proteins in non-aqueous enzymatic processes is their structural dynamics in relation to solvent polarity. A pool of mutants derived from Geobacillus zalihae T1 lipase was screened in organic solvents (methanol, ethanol, propanol, butanol and pentanol) resulting in the selection of six mutants at initial screening (A83D/K251E, R21C, G35D/S195 N, K84R/R103C/M121I/T272 M and R106H/G327S). Site-directed mutagenesis further yielded quadruple mutants A83D/M121I/K251E/G327S and A83D/M121I/S195 N/T272 M, both of which had improved activity after incubation in methanol. The km and kcat values of these mutants vary marginally with the wild-type enzyme in the methanol/substrate mixture. Thermally induced unfolding of mutants was accompanied with some loss of secondary structure content. The root mean square deviations (RMSD) and B-factors revealed that changes in the structural organization are intertwined with an interplay of the protein backbone with organic solvents. Spatially exposed charged residues showed correlations between the solvation dynamics of the methanol solvent and the hydrophobicity of the residues. The short distances of the radial distribution function provided the required distances for hydrogen bond formation and hydrophobic interactions. These dynamic changes demonstrate newly formed structural interactions could be targeted and incorporated experimentally on the basis of solvent mobility and mutant residues.
Collapse
Affiliation(s)
- Jonathan Maiangwa
- Department of Cell and Molecular Biology, Enzyme Microbial Technology Research Center, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia Serdang, 43400, UPM Serdang, Selangor, Malaysia; Department of Microbiology Kaduna State University, Nigeria; Enzyme Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia Serdang, 43400, UPM Serdang, Selangor, Malaysia
| | - Siti Hajar Hamdan
- Department of Biochemistry, Enzyme Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia Serdang, 43400, UPM Serdang, Selangor, Malaysia
| | - Mohd Shukuri Mohamad Ali
- Department of Biochemistry, Enzyme Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia Serdang, 43400, UPM Serdang, Selangor, Malaysia
| | - Abu Bakar Salleh
- Enzyme Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia Serdang, 43400, UPM Serdang, Selangor, Malaysia
| | - Raja Noor Zaliha Raja Abd Rahman
- Department of Microbiology, Enzyme Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia, 43400, UPM Serdang, Selangor, Malaysia
| | - Fairolniza Mohd Shariff
- Institute of Bioscience, 43400, UPM Serdang, Universiti Putra Malaysia Serdang, Selangor, Malaysia
| | - Thean Chor Leow
- Department of Cell and Molecular Biology, Enzyme Microbial Technology Research Center, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia Serdang, 43400, UPM Serdang, Selangor, Malaysia; Enzyme Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Science, Universiti Putra Malaysia Serdang, 43400, UPM Serdang, Selangor, Malaysia; Institute of Bioscience, 43400, UPM Serdang, Universiti Putra Malaysia Serdang, Selangor, Malaysia.
| |
Collapse
|
10
|
The Roles of Protein Structure, Taxon Sampling, and Model Complexity in Phylogenomics: A Case Study Focused on Early Animal Divergences. BIOPHYSICA 2021. [DOI: 10.3390/biophysica1020008] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Despite the long history of using protein sequences to infer the tree of life, the potential for different parts of protein structures to retain historical signal remains unclear. We propose that it might be possible to improve analyses of phylogenomic datasets by incorporating information about protein structure. We test this idea using the position of the root of Metazoa (animals) as a model system. We examined the distribution of “strongly decisive” sites (alignment positions that support a specific tree topology) in a dataset comprising >1500 proteins and almost 100 taxa. The proportion of each class of strongly decisive sites in different structural environments was very sensitive to the model used to analyze the data when a limited number of taxa were used but they were stable when taxa were added. As long as enough taxa were analyzed, sites in all structural environments supported the same topology regardless of whether standard tree searches or decisive sites were used to select the optimal tree. However, the use of decisive sites revealed a difference between the support for minority topologies for sites in different structural environments: buried sites and sites in sheet and coil environments exhibited equal support for the minority topologies, whereas solvent-exposed and helix sites had unequal numbers of sites, supporting the minority topologies. This suggests that the relatively slowly evolving buried, sheet, and coil sites are giving an accurate picture of the true species tree and the amount of conflict among gene trees. Taken as a whole, this study indicates that phylogenetic analyses using sites in different structural environments can yield different topologies for the deepest branches in the animal tree of life and that analyzing larger numbers of taxa eliminates this conflict. More broadly, our results highlight the desirability of incorporating information about protein structure into phylogenomic analyses.
Collapse
|
11
|
Pandey A, Braun EL. Phylogenetic Analyses of Sites in Different Protein Structural Environments Result in Distinct Placements of the Metazoan Root. BIOLOGY 2020; 9:E64. [PMID: 32231097 PMCID: PMC7235752 DOI: 10.3390/biology9040064] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/24/2019] [Revised: 03/09/2020] [Accepted: 03/20/2020] [Indexed: 12/23/2022]
Abstract
Phylogenomics, the use of large datasets to examine phylogeny, has revolutionized the study of evolutionary relationships. However, genome-scale data have not been able to resolve all relationships in the tree of life; this could reflect, at least in part, the poor-fit of the models used to analyze heterogeneous datasets. Some of the heterogeneity may reflect the different patterns of selection on proteins based on their structures. To test that hypothesis, we developed a pipeline to divide phylogenomic protein datasets into subsets based on secondary structure and relative solvent accessibility. We then tested whether amino acids in different structural environments had distinct signals for the topology of the deepest branches in the metazoan tree. We focused on a dataset that appeared to have a mixture of signals and we found that the most striking difference in phylogenetic signal reflected relative solvent accessibility. Analyses of exposed sites (residues located on the surface of proteins) yielded a tree that placed ctenophores sister to all other animals whereas sites buried inside proteins yielded a tree with a sponge+ctenophore clade. These differences in phylogenetic signal were not ameliorated when we conducted analyses using a set of maximum-likelihood profile mixture models. These models are very similar to the Bayesian CAT model, which has been used in many analyses of deep metazoan phylogeny. In contrast, analyses conducted after recoding amino acids to limit the impact of deviations from compositional stationarity increased the congruence in the estimates of phylogeny for exposed and buried sites; after recoding amino acid trees estimated using the exposed and buried site both supported placement of ctenophores sister to all other animals. Although the central conclusion of our analyses is that sites in different structural environments yield distinct trees when analyzed using models of protein evolution, our amino acid recoding analyses also have implications for metazoan evolution. Specifically, our results add to the evidence that ctenophores are the sister group of all other animals and they further suggest that the placozoa+cnidaria clade found in some other studies deserves more attention. Taken as a whole, these results provide striking evidence that it is necessary to achieve a better understanding of the constraints due to protein structure to improve phylogenetic estimation.
Collapse
Affiliation(s)
- Akanksha Pandey
- Department of Biology, University of Florida, Gainesville, FL 32611, USA;
| | - Edward L. Braun
- Department of Biology, University of Florida, Gainesville, FL 32611, USA;
- Genetics Institute, University of Florida, Gainesville, FL 32611, USA
| |
Collapse
|
12
|
Pollack JD, Gerard D, Makhatadze GI, Pearl DK. Evolutionary conservation and structural localizations suggest a physical trace of metabolism’s progressive geochronological emergence. J Biomol Struct Dyn 2019; 38:3700-3719. [DOI: 10.1080/07391102.2019.1679666] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- J. Dennis Pollack
- Department of Molecular Virology, Immunology and Medical Genetics, College of Medicine, The Ohio State University, Columbus, Ohio, USA
| | - David Gerard
- Department of Mathematics and Statistics, American University, Washington, DC, USA
| | - George I. Makhatadze
- Department of Biological Sciences, Rensselaer Polytechnic Institute, Troy, New York, USA
| | - Dennis K. Pearl
- Department of Statistics, Penn State University, University Park, Pennsylvania, USA
| |
Collapse
|
13
|
Abstract
Homologous sequence alignments contain important information about the constraints that shape protein family evolution. Correlated changes between different residues, for instance, can be highly predictive of physical contacts within three-dimensional structures. Detecting such co-evolutionary signals via direct coupling analysis is particularly challenging given the shared phylogenetic history and uneven sampling of different lineages from which protein sequences are derived. Current best practices for mitigating such effects include sequence-identity-based weighting of input sequences and post-hoc re-scaling of evolutionary coupling scores. However, numerous weighting schemes have been previously developed for other applications, and it is unknown whether any of these schemes may better account for phylogenetic artifacts in evolutionary coupling analyses. Here, we show across a dataset of 150 diverse protein families that the current best practices out-perform several alternative sequence- and tree-based weighting methods. Nevertheless, we find that sequence weighting in general provides only a minor benefit relative to post-hoc transformations that re-scale the derived evolutionary couplings. While our findings do not rule out the possibility that an as-yet-untested weighting method may show improved results, the similar predictive accuracies that we observe across distinct weighting methods suggests that there may be little room for further improvement on top of existing strategies.
Collapse
|
14
|
Manguy J, Shields DC. Implications of kappa-casein evolutionary diversity for the self-assembly and aggregation of casein micelles. ROYAL SOCIETY OPEN SCIENCE 2019; 6:190939. [PMID: 31824707 PMCID: PMC6837221 DOI: 10.1098/rsos.190939] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/23/2019] [Accepted: 09/24/2019] [Indexed: 06/10/2023]
Abstract
Milk alpha-, beta- and kappa-casein proteins assemble into casein micelles in breast epithelial cells. The glycomacropeptide (GMP) tails of kappa-casein that extend from the surface of the micelle are key to assembly and aggregation. Aggregation is triggered by stomach pepsin cleavage of GMP from para-kappa-casein (PKC). While one casein micelle model emphasizes the importance of hydrophobic interactions, another focuses on polar residues. We performed an evolutionary analysis of kappa-casein primary sequence and predicted features that potentially impact on protein interactions. We noted more rapid change in the earlier period (166 to 60 Ma). Pepsin and plasmin cleavage sites were avoided in the GMP, which may partly explain its amino acid composition. Short tandem repeats have led to modest expansions of PKC, and to large GMP expansions, suggesting the GMP is less length constrained. Amino acid compositional constraints were assessed across species. Polarity and hydrophobicity properties were insufficient to explain differences between PKC and GMP. Among polar residues, threonine dominates the GMP, compared to serine, probably reflecting its preference for O-glycosylation over phosphorylation. Glutamine, enriched in the bovine PQ-rich region, is not positionally conserved in other species. Among hydrophobic residues, isoleucine is clearly preferred over leucine in the GMP, and patches of hydrophobicity are not markedly positionally conserved. PKC tyrosine and charged residues showed stronger conservation of position, suggesting a role for pi-interactions, seen in other structurally dynamic protein membraneless assemblies. Independent acquisitions of cysteines are consistent with a trend of increasing stabilization of multimers by covalent disulphide bonds, over evolutionary time. In conclusion, kappa-casein compositional and positional constraints appear to be influenced by modification preferences, protease evasion and protein-protein interactions.
Collapse
Affiliation(s)
- Jean Manguy
- UCD Conway Institute, University College Dublin, Belfield, Dublin 4, Ireland
- School of Medicine, University College Dublin, Belfield, Dublin 4, Ireland
- Food for Health Ireland, University College Dublin, Belfield, Dublin 4, Ireland
| | - Denis C. Shields
- UCD Conway Institute, University College Dublin, Belfield, Dublin 4, Ireland
- School of Medicine, University College Dublin, Belfield, Dublin 4, Ireland
- Food for Health Ireland, University College Dublin, Belfield, Dublin 4, Ireland
| |
Collapse
|
15
|
Tran VDT, Moretti S, Coste AT, Amorim-Vaz S, Sanglard D, Pagni M. Condition-specific series of metabolic sub-networks and its application for gene set enrichment analysis. Bioinformatics 2019; 35:2258-2266. [PMID: 30445518 PMCID: PMC6596900 DOI: 10.1093/bioinformatics/bty929] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2018] [Revised: 10/16/2018] [Accepted: 11/09/2018] [Indexed: 11/12/2022] Open
Abstract
MOTIVATION Genome-scale metabolic networks and transcriptomic data represent complementary sources of knowledge about an organism's metabolism, yet their integration to achieve biological insight remains challenging. RESULTS We investigate here condition-specific series of metabolic sub-networks constructed by successively removing genes from a comprehensive network. The optimal order of gene removal is deduced from transcriptomic data. The sub-networks are evaluated via a fitness function, which estimates their degree of alteration. We then consider how a gene set, i.e. a group of genes contributing to a common biological function, is depleted in different series of sub-networks to detect the difference between experimental conditions. The method, named metaboGSE, is validated on public data for Yarrowia lipolytica and mouse. It is shown to produce GO terms of higher specificity compared to popular gene set enrichment methods like GSEA or topGO. AVAILABILITY AND IMPLEMENTATION The metaboGSE R package is available at https://CRAN.R-project.org/package=metaboGSE. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Van Du T Tran
- Vital-IT Group, SIB Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Sébastien Moretti
- Vital-IT Group, SIB Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
- Evolutionary Bioinformatics Group, SIB Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Alix T Coste
- Institute of Microbiology, University Hospital and University of Lausanne, CH-1015 Lausanne, Switzerland
| | - Sara Amorim-Vaz
- Institute of Microbiology, University Hospital and University of Lausanne, CH-1015 Lausanne, Switzerland
| | - Dominique Sanglard
- Institute of Microbiology, University Hospital and University of Lausanne, CH-1015 Lausanne, Switzerland
| | - Marco Pagni
- Vital-IT Group, SIB Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| |
Collapse
|
16
|
Guthrie L, Wolfson S, Kelly L. The human gut chemical landscape predicts microbe-mediated biotransformation of foods and drugs. eLife 2019; 8:42866. [PMID: 31184303 PMCID: PMC6559788 DOI: 10.7554/elife.42866] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2018] [Accepted: 05/10/2019] [Indexed: 12/14/2022] Open
Abstract
Microbes are nature's chemists, capable of producing and metabolizing a diverse array of compounds. In the human gut, microbial biochemistry can be beneficial, for example vitamin production and complex carbohydrate breakdown; or detrimental, such as the reactivation of an inactive drug metabolite leading to patient toxicity. Identifying clinically relevant microbiome metabolism requires linking microbial biochemistry and ecology with patient outcomes. Here we present MicrobeFDT, a resource which clusters chemically similar drug and food compounds and links these compounds to microbial enzymes and known toxicities. We demonstrate that compound structural similarity can serve as a proxy for toxicity, enzyme sharing, and coarse-grained functional similarity. MicrobeFDT allows users to flexibly interrogate microbial metabolism, compounds of interest, and toxicity profiles to generate novel hypotheses of microbe-diet-drug-phenotype interactions that influence patient outcomes. We validate one such hypothesis experimentally, using MicrobeFDT to reveal unrecognized gut microbiome metabolism of the ovarian cancer drug altretamine.
Collapse
Affiliation(s)
- Leah Guthrie
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, New York, United States
| | - Sarah Wolfson
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, New York, United States
| | - Libusha Kelly
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, New York, United States.,Department of Microbiology and Immunology, Albert Einstein College of Medicine, New York, United States
| |
Collapse
|
17
|
Aun E, Brauer A, Kisand V, Tenson T, Remm M. A k-mer-based method for the identification of phenotype-associated genomic biomarkers and predicting phenotypes of sequenced bacteria. PLoS Comput Biol 2018; 14:e1006434. [PMID: 30346947 PMCID: PMC6211763 DOI: 10.1371/journal.pcbi.1006434] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2018] [Revised: 11/01/2018] [Accepted: 08/15/2018] [Indexed: 11/18/2022] Open
Abstract
We have developed an easy-to-use and memory-efficient method called PhenotypeSeeker that (a) identifies phenotype-specific k-mers, (b) generates a k-mer-based statistical model for predicting a given phenotype and (c) predicts the phenotype from the sequencing data of a given bacterial isolate. The method was validated on 167 Klebsiella pneumoniae isolates (virulence), 200 Pseudomonas aeruginosa isolates (ciprofloxacin resistance) and 459 Clostridium difficile isolates (azithromycin resistance). The phenotype prediction models trained from these datasets obtained the F1-measure of 0.88 on the K. pneumoniae test set, 0.88 on the P. aeruginosa test set and 0.97 on the C. difficile test set. The F1-measures were the same for assembled sequences and raw sequencing data; however, building the model from assembled genomes is significantly faster. On these datasets, the model building on a mid-range Linux server takes approximately 3 to 5 hours per phenotype if assembled genomes are used and 10 hours per phenotype if raw sequencing data are used. The phenotype prediction from assembled genomes takes less than one second per isolate. Thus, PhenotypeSeeker should be well-suited for predicting phenotypes from large sequencing datasets. PhenotypeSeeker is implemented in Python programming language, is open-source software and is available at GitHub (https://github.com/bioinfo-ut/PhenotypeSeeker/).
Collapse
Affiliation(s)
- Erki Aun
- Department of Bioinformatics, Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia
- * E-mail:
| | - Age Brauer
- Department of Bioinformatics, Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia
| | - Veljo Kisand
- Institute of Technology, University of Tartu, Tartu, Estonia
| | - Tanel Tenson
- Institute of Technology, University of Tartu, Tartu, Estonia
| | - Maido Remm
- Department of Bioinformatics, Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia
| |
Collapse
|
18
|
Rivas E, Clements J, Eddy SR. A statistical test for conserved RNA structure shows lack of evidence for structure in lncRNAs. Nat Methods 2016; 14:45-48. [PMID: 27819659 PMCID: PMC5554622 DOI: 10.1038/nmeth.4066] [Citation(s) in RCA: 231] [Impact Index Per Article: 28.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2016] [Accepted: 09/14/2016] [Indexed: 12/14/2022]
Abstract
Many functional RNAs have an evolutionarily conserved secondary structure. Conservation of RNA base pairing induces pairwise covariations in sequence alignments. We developed a computational method, R-scape (RNA Structural Covariation Above Phylogenetic Expectation), that quantitatively tests whether covariation analysis supports the presence of a conserved RNA secondary structure. R-scape analysis finds no statistically significant support for proposed secondary structures of the long noncoding RNAs HOTAIR, SRA, and Xist.
Collapse
Affiliation(s)
- Elena Rivas
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, Massachusetts, USA
| | - Jody Clements
- Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, Virginia, USA
| | - Sean R Eddy
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, Massachusetts, USA.,Howard Hughes Medical Institute, Harvard University, Cambridge, Massachusetts, USA.,FAS Center for Systems Biology, Harvard University, Cambridge, Massachusetts, USA.,John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts, USA
| |
Collapse
|
19
|
Korostelev YD, Zharov IA, Mironov AA, Rakhmaininova AB, Gelfand MS. Identification of Position-Specific Correlations between DNA-Binding Domains and Their Binding Sites. Application to the MerR Family of Transcription Factors. PLoS One 2016; 11:e0162681. [PMID: 27690309 PMCID: PMC5045206 DOI: 10.1371/journal.pone.0162681] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2015] [Accepted: 08/26/2016] [Indexed: 11/25/2022] Open
Abstract
The large and increasing volume of genomic data analyzed by comparative methods provides information about transcription factors and their binding sites that, in turn, enables statistical analysis of correlations between factors and sites, uncovering mechanisms and evolution of specific protein-DNA recognition. Here we present an online tool, Prot-DNA-Korr, designed to identify and analyze crucial protein-DNA pairs of positions in a family of transcription factors. Correlations are identified by analysis of mutual information between columns of protein and DNA alignments. The algorithm reduces the effects of common phylogenetic history and of abundance of closely related proteins and binding sites. We apply it to five closely related subfamilies of the MerR family of bacterial transcription factors that regulate heavy metal resistance systems. We validate the approach using known 3D structures of MerR-family proteins in complexes with their cognate DNA binding sites and demonstrate that a significant fraction of correlated positions indeed form specific side-chain-to-base contacts. The joint distribution of amino acids and nucleotides hence may be used to predict changes of specificity for point mutations in transcription factors.
Collapse
Affiliation(s)
- Yuriy D. Korostelev
- A.A. Kharkevich Institute for Information Transmission Problems, Russian Academy of Sciences, 19-1 Bolshoy Karetny pereulok, Moscow, Russia, 127994
- Department of Bioengineering and Bioinformatics, Moscow State University, 1-73 Vorobievy Gory, Moscow, Russia, 119991
| | - Ilya A. Zharov
- A.A. Kharkevich Institute for Information Transmission Problems, Russian Academy of Sciences, 19-1 Bolshoy Karetny pereulok, Moscow, Russia, 127994
| | - Andrey A. Mironov
- A.A. Kharkevich Institute for Information Transmission Problems, Russian Academy of Sciences, 19-1 Bolshoy Karetny pereulok, Moscow, Russia, 127994
- Department of Bioengineering and Bioinformatics, Moscow State University, 1-73 Vorobievy Gory, Moscow, Russia, 119991
| | - Alexandra B. Rakhmaininova
- A.A. Kharkevich Institute for Information Transmission Problems, Russian Academy of Sciences, 19-1 Bolshoy Karetny pereulok, Moscow, Russia, 127994
| | - Mikhail S. Gelfand
- A.A. Kharkevich Institute for Information Transmission Problems, Russian Academy of Sciences, 19-1 Bolshoy Karetny pereulok, Moscow, Russia, 127994
- Department of Bioengineering and Bioinformatics, Moscow State University, 1-73 Vorobievy Gory, Moscow, Russia, 119991
- * E-mail:
| |
Collapse
|
20
|
Shahmoradi A, Wilke CO. Dissecting the roles of local packing density and longer-range effects in protein sequence evolution. Proteins 2016; 84:841-54. [PMID: 26990194 PMCID: PMC5292938 DOI: 10.1002/prot.25034] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2015] [Revised: 02/01/2016] [Accepted: 02/24/2016] [Indexed: 11/07/2022]
Abstract
What are the structural determinants of protein sequence evolution? A number of site-specific structural characteristics have been proposed, most of which are broadly related to either the density of contacts or the solvent accessibility of individual residues. Most importantly, there has been disagreement in the literature over the relative importance of solvent accessibility and local packing density for explaining site-specific sequence variability in proteins. We show that this discussion has been confounded by the definition of local packing density. The most commonly used measures of local packing, such as contact number and the weighted contact number, represent the combined effects of local packing density and longer-range effects. As an alternative, we propose a truly local measure of packing density around a single residue, based on the Voronoi cell volume. We show that the Voronoi cell volume, when calculated relative to the geometric center of amino-acid side chains, behaves nearly identically to the relative solvent accessibility, and each individually can explain, on average, approximately 34% of the site-specific variation in evolutionary rate in a data set of 209 enzymes. An additional 10% of variation can be explained by nonlocal effects that are captured in the weighted contact number. Consequently, evolutionary variation at a site is determined by the combined effects of the immediate amino-acid neighbors of that site and effects mediated by more distant amino acids. We conclude that instead of contrasting solvent accessibility and local packing density, future research should emphasize on the relative importance of immediate contacts and longer-range effects on evolutionary variation. Proteins 2016; 84:841-854. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Amir Shahmoradi
- Department of Physics, The University of Texas at Austin
- Center for Computational Biology and Bioinformatics, The University
of Texas at Austin
- Institute for Cellular and Molecular Biology, The University of
Texas at Austin
| | - Claus O. Wilke
- Center for Computational Biology and Bioinformatics, The University
of Texas at Austin
- Institute for Cellular and Molecular Biology, The University of
Texas at Austin
- Department of Integrative Biology, The University of Texas at
Austin
| |
Collapse
|
21
|
Abstract
Although protein folding and stability have been well explored under simplified conditions in vitro, it is yet unclear how these basic self-organization events are modulated by the crowded interior of live cells. To find out, we use here in-cell NMR to follow at atomic resolution the thermal unfolding of a β-barrel protein inside mammalian and bacterial cells. Challenging the view from in vitro crowding effects, we find that the cells destabilize the protein at 37 °C but with a conspicuous twist: While the melting temperature goes down the cold unfolding moves into the physiological regime, coupled to an augmented heat-capacity change. The effect seems induced by transient, sequence-specific, interactions with the cellular components, acting preferentially on the unfolded ensemble. This points to a model where the in vivo influence on protein behavior is case specific, determined by the individual protein's interplay with the functionally optimized "interaction landscape" of the cellular interior.
Collapse
|
22
|
Janssen S, Giegerich R. Ambivalent covariance models. BMC Bioinformatics 2015; 16:178. [PMID: 26017195 PMCID: PMC4504443 DOI: 10.1186/s12859-015-0569-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2014] [Accepted: 04/10/2015] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Evolutionary variations let us define a set of similar nucleic acid sequences as a family if these different molecules execute a common function. Capturing their sequence variation by using e. g. position specific scoring matrices significantly improves sensitivity of detection tools. Members of a functional (non-coding) RNA family are affected by these variations not only on the sequence, but also on the structural level. For example, some transfer-RNAs exhibit a fifth helix in addition to the typical cloverleaf structure. Current covariance models - the unrivaled homology search approach for structured RNA - do not benefit from structural variation within a family, but rather penalize it. This leads to artificial subdivision of families and loss of information in the RFAM database. RESULTS We propose an extension to the fundamental architecture of covariance models to allow for several, compatible consensus structures. The resulting models are called ambivalent covariance models. Evaluation on several RFAM families shows that coalescence of structural variation within a family by using ambivalent consensus models is superior to subdividing the family into multiple classical covariance models. CONCLUSION A prototype and source code is available at http://bibiserv.cebitec.uni-bielefeld.de/acms.
Collapse
Affiliation(s)
- Stefan Janssen
- Practical Computer Science, Faculty of Technology, Bielefeld University, Universitätsstraße 25, Bielefeld, 33615, Germany.
| | - Robert Giegerich
- Practical Computer Science, Faculty of Technology, Bielefeld University, Universitätsstraße 25, Bielefeld, 33615, Germany.
| |
Collapse
|
23
|
Abstract
The molecular mechanisms underlying pressure-induced protein denaturation can be analyzed based on the pressure-dependent differences in the apparent volume occupied by amino acids inside the protein and when exposed to water in an unfolded conformation. This chapter presents a volumetric analysis of the peptide group and the 20 naturally occurring amino acid side chains in the interior of the native state, the micelle-like interior of the pressure-induced denatured state, and in the unfolded conformation modeled by low-molecular analogs of proteins. The transfer of a peptide group from the protein interior to water becomes increasingly favorable as pressure increases. This observation classifies solvation of peptide groups as a major driving force in pressure-induced protein denaturation. Polar side chains do not appear to exhibit significant pressure-dependent changes in their preference for the protein interior or solvent. The transfer of nonpolar side chains from the protein interior to water becomes more unfavorable as pressure increases. An inference can be drawn that a sizeable population of nonpolar side chains remains buried inside a solvent-inaccessible core of the pressure-induced denatured state. At elevated pressures this core, owing to the absence of structural constraints, may become packed almost as tightly as the interior of the native state. The presence and partial disappearance of large intraglobular voids is another driving force facilitating pressure-induced protein denaturation. Volumetric data presented here have implications for the kinetics of protein folding and shed light on the nature of the folding transition state ensembles.
Collapse
|
24
|
Valach M, Burger G, Gray MW, Lang BF. Widespread occurrence of organelle genome-encoded 5S rRNAs including permuted molecules. Nucleic Acids Res 2014; 42:13764-77. [PMID: 25429974 PMCID: PMC4267664 DOI: 10.1093/nar/gku1266] [Citation(s) in RCA: 95] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2014] [Revised: 11/14/2014] [Accepted: 11/17/2014] [Indexed: 11/13/2022] Open
Abstract
5S Ribosomal RNA (5S rRNA) is a universal component of ribosomes, and the corresponding gene is easily identified in archaeal, bacterial and nuclear genome sequences. However, organelle gene homologs (rrn5) appear to be absent from most mitochondrial and several chloroplast genomes. Here, we re-examine the distribution of organelle rrn5 by building mitochondrion- and plastid-specific covariance models (CMs) with which we screened organelle genome sequences. We not only recover all organelle rrn5 genes annotated in GenBank records, but also identify more than 50 previously unrecognized homologs in mitochondrial genomes of various stramenopiles, red algae, cryptomonads, malawimonads and apusozoans, and surprisingly, in the apicoplast (highly derived plastid) genomes of the coccidian pathogens Toxoplasma gondii and Eimeria tenella. Comparative modeling of RNA secondary structure reveals that mitochondrial 5S rRNAs from brown algae adopt a permuted triskelion shape that has not been seen elsewhere. Expression of the newly predicted rrn5 genes is confirmed experimentally in 10 instances, based on our own and published RNA-Seq data. This study establishes that particularly mitochondrial 5S rRNA has a much broader taxonomic distribution and a much larger structural variability than previously thought. The newly developed CMs will be made available via the Rfam database and the MFannot organelle genome annotator.
Collapse
MESH Headings
- Coccidia/genetics
- Databases, Nucleic Acid
- Genes, Mitochondrial
- Genes, rRNA
- Genome, Mitochondrial
- Genome, Plastid
- Nucleic Acid Conformation
- Phaeophyceae/genetics
- RNA/chemistry
- RNA/genetics
- RNA, Mitochondrial
- RNA, Ribosomal, 5S/chemistry
- RNA, Ribosomal, 5S/classification
- RNA, Ribosomal, 5S/genetics
- Sequence Analysis, RNA
- Stramenopiles/genetics
Collapse
Affiliation(s)
- Matus Valach
- Department of Biochemistry and Robert-Cedergren Centre of Bioinformatics and Genomics, Université de Montréal, Montréal, QC, H3C 3J7, Canada
| | - Gertraud Burger
- Department of Biochemistry and Robert-Cedergren Centre of Bioinformatics and Genomics, Université de Montréal, Montréal, QC, H3C 3J7, Canada
| | - Michael W Gray
- Department of Biochemistry and Molecular Biology, Dalhousie University, Halifax, NS, B3H 4B2, Canada
| | - B Franz Lang
- Department of Biochemistry and Robert-Cedergren Centre of Bioinformatics and Genomics, Université de Montréal, Montréal, QC, H3C 3J7, Canada
| |
Collapse
|
25
|
Zaucha J, Stahlhacke J, Oates ME, Thurlby N, Rackham OJL, Fang H, Smithers B, Gough J. A Proteome Quality Index. Environ Microbiol 2014; 17:4-9. [DOI: 10.1111/1462-2920.12622] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Affiliation(s)
- Jan Zaucha
- Department of Computer Science; University of Bristol; Bristol UK
- Bristol Centre for Complexity Sciences; University of Bristol; Bristol UK
| | | | - Matt E. Oates
- Department of Computer Science; University of Bristol; Bristol UK
| | - Natalie Thurlby
- Department of Computer Science; University of Bristol; Bristol UK
- Bristol Centre for Complexity Sciences; University of Bristol; Bristol UK
| | - Owen J. L. Rackham
- Medical Research Council Clinical Sciences Centre; Faculty of Medicine; Hammersmith Hospital; Imperial College London; London UK
| | - Hai Fang
- Department of Computer Science; University of Bristol; Bristol UK
| | - Ben Smithers
- Department of Computer Science; University of Bristol; Bristol UK
| | - Julian Gough
- Department of Computer Science; University of Bristol; Bristol UK
| |
Collapse
|
26
|
Pervouchine DD. IRBIS: a systematic search for conserved complementarity. RNA (NEW YORK, N.Y.) 2014; 20:1519-31. [PMID: 25142064 PMCID: PMC4174434 DOI: 10.1261/rna.045088.114] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Accepted: 06/26/2014] [Indexed: 05/28/2023]
Abstract
IRBIS is a computational pipeline for detecting conserved complementary regions in unaligned orthologous sequences. Unlike other methods, it follows the "first-fold-then-align" principle in which all possible combinations of complementary k-mers are searched for simultaneous conservation. The novel trimming procedure reduces the size of the search space and improves the performance to the point where large-scale analyses of intra- and intermolecular RNA-RNA interactions become possible. In this article, I provide a rigorous description of the method, benchmarking on simulated and real data, and a set of stringent predictions of intramolecular RNA structure in placental mammals, drosophilids, and nematodes. I discuss two particular cases of long-range RNA structures that are likely to have a causal effect on single- and multiple-exon skipping, one in the mammalian gene Dystonin and the other in the insect gene Ca-α1D. In Dystonin, one of the two complementary boxes contains a binding site of Rbfox protein similar to one recently described in Enah gene. I also report that snoRNAs and long noncoding RNAs (lncRNAs) have a high capacity of base-pairing to introns of protein-coding genes, suggesting possible involvement of these transcripts in splicing regulation. I also find that conserved sequences that occur equally likely on both strands of DNA (e.g., transcription factor binding sites) contribute strongly to the false-discovery rate and, therefore, would confound every such analysis. IRBIS is an open-source software that is available at http://genome.crg.es/~dmitri/irbis/.
Collapse
Affiliation(s)
- Dmitri D Pervouchine
- Centre for Genomic Regulation and UPF, Barcelona 08003, Spain Faculty of Bioengineering and Bioinformatics, Moscow State University, 119992 Moscow, Russia
| |
Collapse
|
27
|
Nelson JW, Sudarsan N, Furukawa K, Weinberg Z, Wang JX, Breaker RR. Riboswitches in eubacteria sense the second messenger c-di-AMP. Nat Chem Biol 2013; 9:834-9. [PMID: 24141192 PMCID: PMC3830699 DOI: 10.1038/nchembio.1363] [Citation(s) in RCA: 210] [Impact Index Per Article: 19.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2013] [Accepted: 09/10/2013] [Indexed: 11/28/2022]
Abstract
Cyclic di-adenosine monophosphate (c-di-AMP) is a recently discovered bacterial second messenger implicated in the control of cell wall metabolism, osmotic stress responses and sporulation. However, the mechanisms by which c-di-AMP triggers these physiological responses have remained largely unknown. Notably, a candidate riboswitch class called ydaO associates with numerous genes involved in these same processes. Although a representative ydaO motif RNA recently was reported to weakly bind ATP, we report that numerous members of this noncoding RNA class selectively respond to c-di-AMP with subnanomolar affinity. Our findings resolve the mystery regarding the primary ligand for this extremely common riboswitch class and expose a major portion of the super-regulon of genes that are controlled by the widespread bacterial second messenger c-di-AMP.
Collapse
Affiliation(s)
- James W. Nelson
- Department of Chemistry, Yale University, Box 208107, New Haven, CT 06520, USA
| | - Narasimhan Sudarsan
- Howard Hughes Medical Institute, Yale University, Box 208103, New Haven, CT 06520, USA
- Department of Molecular, Cellular and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA
| | - Kazuhiro Furukawa
- Department of Molecular, Cellular and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA
| | - Zasha Weinberg
- Howard Hughes Medical Institute, Yale University, Box 208103, New Haven, CT 06520, USA
- Department of Molecular, Cellular and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA
| | - Joy X. Wang
- Department of Molecular, Cellular and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA
| | - Ronald R. Breaker
- Howard Hughes Medical Institute, Yale University, Box 208103, New Haven, CT 06520, USA
- Department of Molecular, Cellular and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA
- Department of Molecular Biophysics and Biochemistry, Yale University, Box 208103, New Haven, CT 06520, USA
| |
Collapse
|
28
|
Kim S, Jin D, Lee H. Predicting drug-target interactions using drug-drug interactions. PLoS One 2013; 8:e80129. [PMID: 24278248 PMCID: PMC3836969 DOI: 10.1371/journal.pone.0080129] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2013] [Accepted: 09/30/2013] [Indexed: 11/18/2022] Open
Abstract
Computational methods for predicting drug-target interactions have become important in drug research because they can help to reduce the time, cost, and failure rates for developing new drugs. Recently, with the accumulation of drug-related data sets related to drug side effects and pharmacological data, it has became possible to predict potential drug-target interactions. In this study, we focus on drug-drug interactions (DDI), their adverse effects ([Formula: see text]) and pharmacological information ([Formula: see text]), and investigate the relationship among chemical structures, side effects, and DDIs from several data sources. In this study, [Formula: see text] data from the STITCH database, [Formula: see text] from drugs.com, and drug-target pairs from ChEMBL and SIDER were first collected. Then, by applying two machine learning approaches, a support vector machine (SVM) and a kernel-based L1-norm regularized logistic regression (KL1LR), we showed that DDI is a promising feature in predicting drug-target interactions. Next, the accuracies of predicting drug-target interactions using DDI were compared to those obtained using the chemical structure and side effects based on the SVM and KL1LR approaches, showing that DDI was the data source contributing the most for predicting drug-target interactions.
Collapse
Affiliation(s)
- Shinhyuk Kim
- School of Information and Communications, Gwangju Institute of Science and Technology, Gwangju, South Korea
| | - Daeyong Jin
- School of Information and Communications, Gwangju Institute of Science and Technology, Gwangju, South Korea
| | - Hyunju Lee
- School of Information and Communications, Gwangju Institute of Science and Technology, Gwangju, South Korea
- * E-mail:
| |
Collapse
|
29
|
Abstract
BACKGROUND Many genome-wide association studies focus on associating single loci with target phenotypes. However, in the setting of rare variation, accumulating sufficient samples to assess these associations can be difficult. Moreover, multiple variations in a gene or a set of genes within a pathway may all contribute to the phenotype, suggesting that the aggregation of variations found over the gene or pathway may be useful for improving the power to detect associations. RESULTS Here, we present a method for aggregating single nucleotide polymorphisms (SNPs) along biologically relevant pathways in order to seek genetic associations with phenotypes. Our method uses all available genetic variants and does not remove those in linkage disequilibrium (LD). Instead, it uses a novel SNP weighting scheme to down-weight the contributions of correlated SNPs. We apply our method to three cohorts of patients taking warfarin: two European descent cohorts and an African American cohort. Although the clinical covariates and key pharmacogenetic loci for warfarin have been characterized, our association metric identifies a significant association with mutations distributed throughout the pathway of warfarin metabolism. We improve dose prediction after using all known clinical covariates and pharmacogenetic variants in VKORC1 and CYP2C9. In particular, we find that at least 1% of the missing heritability in warfarin dose may be due to the aggregated effects of variations in the warfarin metabolic pathway, even though the SNPs do not individually show a significant association. CONCLUSIONS Our method allows researchers to study aggregative SNP effects in an unbiased manner by not preselecting SNPs. It retains all the available information by accounting for LD-structure through weighting, which eliminates the need for LD pruning.
Collapse
|
30
|
Osipovitch DC, Parker AS, Makokha CD, Desrosiers J, Kett WC, Moise L, Bailey-Kellogg C, Griswold KE. Design and analysis of immune-evading enzymes for ADEPT therapy. Protein Eng Des Sel 2012; 25:613-23. [PMID: 22898588 DOI: 10.1093/protein/gzs044] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
The unparalleled specificity and activity of therapeutic proteins has reshaped many aspects of modern clinical practice, and aggressive development of new protein drugs promises a continued revolution in disease therapy. As a result of their biological origins, however, therapeutic proteins present unique design challenges for the biomolecular engineer. For example, protein drugs are subject to immune surveillance within the patient's body; this anti-drug immune response can compromise therapeutic efficacy and even threaten patient safety. Thus, there is a growing demand for broadly applicable protein deimmunization strategies. We have recently developed optimization algorithms that integrate computational prediction of T-cell epitopes and bioinformatics-based assessment of the structural and functional consequences of epitope-deleting mutations. Here, we describe the first experimental validation of our deimmunization algorithms using Enterobacter cloacae P99 β-lactamase, a component of antibody-directed enzyme prodrug cancer therapies. Compared with wild-type or a previously deimmunized variant, our computationally optimized sequences exhibited significantly less in vitro binding to human type II major histocompatibility complex immune molecules. At the same time, our globally optimal design exhibited wild-type catalytic proficiency. We conclude that our deimmunization algorithms guide the protein engineer towards promising immunoevasive candidates and thereby have the potential to streamline biotherapeutic development.
Collapse
Affiliation(s)
- Daniel C Osipovitch
- Program in Experimental and Molecular Medicine, Geisel School of Medicine, Dartmouth, Hanover, NH 03784, USA
| | | | | | | | | | | | | | | |
Collapse
|
31
|
Kalinina OV, Oberwinkler H, Glass B, Kräusslich HG, Russell RB, Briggs JAG. Computational identification of novel amino-acid interactions in HIV Gag via correlated evolution. PLoS One 2012; 7:e42468. [PMID: 22879995 PMCID: PMC3411748 DOI: 10.1371/journal.pone.0042468] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2012] [Accepted: 07/09/2012] [Indexed: 12/31/2022] Open
Abstract
Pairs of amino acid positions that evolve in a correlated manner are proposed to play important roles in protein structure or function. Methods to detect them might fare better with families for which sequences of thousands of closely related homologs are available than families with only a few distant relatives. We applied co-evolution analysis to thousands of sequences of HIV Gag, finding that the most significantly co-evolving positions are proximal in the quaternary structures of the viral capsid. A reduction in infectivity caused by mutating one member of a significant pair could be rescued by a compensatory mutation of the other.
Collapse
Affiliation(s)
- Olga V. Kalinina
- CellNetworks, Bioquant, University of Heidelberg, Heidelberg, Germany
| | - Heike Oberwinkler
- Department of Infectious Diseases, Virology, Universitätsklinikum Heidelberg, Heidelberg, Germany
| | - Bärbel Glass
- Department of Infectious Diseases, Virology, Universitätsklinikum Heidelberg, Heidelberg, Germany
| | - Hans-Georg Kräusslich
- CellNetworks, Bioquant, University of Heidelberg, Heidelberg, Germany
- Department of Infectious Diseases, Virology, Universitätsklinikum Heidelberg, Heidelberg, Germany
| | - Robert B. Russell
- CellNetworks, Bioquant, University of Heidelberg, Heidelberg, Germany
| | - John A. G. Briggs
- Structural and Computational Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany
| |
Collapse
|
32
|
Nikulova AA, Favorov AV, Sutormin RA, Makeev VJ, Mironov AA. CORECLUST: identification of the conserved CRM grammar together with prediction of gene regulation. Nucleic Acids Res 2012; 40:e93. [PMID: 22422836 PMCID: PMC3384346 DOI: 10.1093/nar/gks235] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Identification of transcriptional regulatory regions and tracing their internal organization are important for understanding the eukaryotic cell machinery. Cis-regulatory modules (CRMs) of higher eukaryotes are believed to possess a regulatory ‘grammar’, or preferred arrangement of binding sites, that is crucial for proper regulation and thus tends to be evolutionarily conserved. Here, we present a method CORECLUST (COnservative REgulatory CLUster STructure) that predicts CRMs based on a set of positional weight matrices. Given regulatory regions of orthologous and/or co-regulated genes, CORECLUST constructs a CRM model by revealing the conserved rules that describe the relative location of binding sites. The constructed model may be consequently used for the genome-wide prediction of similar CRMs, and thus detection of co-regulated genes, and for the investigation of the regulatory grammar of the system. Compared with related methods, CORECLUST shows better performance at identification of CRMs conferring muscle-specific gene expression in vertebrates and early-developmental CRMs in Drosophila.
Collapse
Affiliation(s)
- Anna A Nikulova
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, 1-73 Leninskie Gory, Moscow 119991, Russia.
| | | | | | | | | |
Collapse
|
33
|
Abstract
An overwhelming array of structural variants has evolved from a comparatively small number of protein structural domains; which has in turn facilitated an expanse of functional derivatives. Herein, I review the primary mechanisms which have contributed to the vastness of our existing, and expanding, protein repertoires. Protein function prediction strategies, both sequence and structure based, are also discussed and their associated strengths and weaknesses assessed.
Collapse
Affiliation(s)
- Roy D Sleator
- Department of Biological Sciences, Cork Institute of Technology, Cork, Ireland.
| |
Collapse
|
34
|
Abstract
The recent explosion in the number and diversity of novel proteins identified by the large-scale "omics" technologies poses new and important questions to the blossoming field of systems biology--what are all these proteins, how did they come about, and most importantly, what do they do? From a comparatively small number of protein structural domains a staggering array of structural variants has evolved, which has in turn facilitated an expanse of functional derivatives. This review considers the primary mechanisms that have contributed to the vastness of our existing, and expanding, protein repertoires, while also outlining the protocols available for elucidating their true biological function. The various function prediction programs available, both sequence and structure based, are discussed and their associated strengths and weaknesses outlined.
Collapse
Affiliation(s)
- Roy D Sleator
- Department of Biological Sciences, Cork Institute of Technology, Bishopstown, Cork, Ireland.
| |
Collapse
|
35
|
Baker JL, Sudarsan N, Weinberg Z, Roth A, Stockbridge RB, Breaker RR. Widespread genetic switches and toxicity resistance proteins for fluoride. Science 2011; 335:233-235. [PMID: 22194412 DOI: 10.1126/science.1215063] [Citation(s) in RCA: 316] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
Most riboswitches are metabolite-binding RNA structures located in bacterial messenger RNAs where they control gene expression. We have discovered a riboswitch class in many bacterial and archaeal species whose members are selectively triggered by fluoride but reject other small anions, including chloride. These fluoride riboswitches activate expression of genes that encode putative fluoride transporters, enzymes that are known to be inhibited by fluoride, and additional proteins of unknown function. Our findings indicate that most organisms are naturally exposed to toxic levels of fluoride and that many species use fluoride-sensing RNAs to control the expression of proteins that alleviate the deleterious effects of this anion.
Collapse
Affiliation(s)
- Jenny L Baker
- Department of Chemistry, Yale University, Box 208103, New Haven, CT 06520, USA
| | - Narasimhan Sudarsan
- Howard Hughes Medical Institute, Yale University, Box 208103, New Haven, CT 06520, USA.,Department of Molecular, Cellular, and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA
| | - Zasha Weinberg
- Howard Hughes Medical Institute, Yale University, Box 208103, New Haven, CT 06520, USA.,Department of Molecular, Cellular, and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA
| | - Adam Roth
- Howard Hughes Medical Institute, Yale University, Box 208103, New Haven, CT 06520, USA.,Department of Molecular, Cellular, and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA
| | - Randy B Stockbridge
- Department of Biochemistry and Howard Hughes Medical Institute, Brandeis University, Waltham, MA 02454, USA
| | - Ronald R Breaker
- Howard Hughes Medical Institute, Yale University, Box 208103, New Haven, CT 06520, USA.,Department of Molecular, Cellular, and Developmental Biology, Yale University, Box 208103, New Haven, CT 06520, USA.,Department of Molecular Biophysics and Biochemistry, Yale University, Box 208103, New Haven, CT 06520, USA
| |
Collapse
|
36
|
Ye X, Wang G, Altschul SF. An assessment of substitution scores for protein profile-profile comparison. Bioinformatics 2011; 27:3356-63. [PMID: 21998158 PMCID: PMC3232366 DOI: 10.1093/bioinformatics/btr565] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2011] [Revised: 09/22/2011] [Accepted: 10/06/2011] [Indexed: 11/14/2022] Open
Abstract
MOTIVATION Pairwise protein sequence alignments are generally evaluated using scores defined as the sum of substitution scores for aligning amino acids to one another, and gap scores for aligning runs of amino acids in one sequence to null characters inserted into the other. Protein profiles may be abstracted from multiple alignments of protein sequences, and substitution and gap scores have been generalized to the alignment of such profiles either to single sequences or to other profiles. Although there is widespread agreement on the general form substitution scores should take for profile-sequence alignment, little consensus has been reached on how best to construct profile-profile substitution scores, and a large number of these scoring systems have been proposed. Here, we assess a variety of such substitution scores. For this evaluation, given a gold standard set of multiple alignments, we calculate the probability that a profile column yields a higher substitution score when aligned to a related than to an unrelated column. We also generalize this measure to sets of two or three adjacent columns. This simple approach has the advantages that it does not depend primarily upon the gold-standard alignment columns with the weakest empirical support, and that it does not need to fit gap and offset costs for use with each substitution score studied. RESULTS A simple symmetrization of mean profile-sequence scores usually performed the best. These were followed closely by several specific scoring systems constructed using a variety of rationales. CONTACT altschul@ncbi.nlm.nih.gov SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Xugang Ye
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894, USA
| | | | | |
Collapse
|
37
|
Yip KY, Utz L, Sitwell S, Hu X, Sidhu SS, Turk BE, Gerstein M, Kim PM. Identification of specificity determining residues in peptide recognition domains using an information theoretic approach applied to large-scale binding maps. BMC Biol 2011; 9:53. [PMID: 21835011 PMCID: PMC3224579 DOI: 10.1186/1741-7007-9-53] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2011] [Accepted: 08/11/2011] [Indexed: 01/06/2023] Open
Abstract
Background Peptide Recognition Domains (PRDs) are commonly found in signaling proteins. They mediate protein-protein interactions by recognizing and binding short motifs in their ligands. Although a great deal is known about PRDs and their interactions, prediction of PRD specificities remains largely an unsolved problem. Results We present a novel approach to identifying these Specificity Determining Residues (SDRs). Our algorithm generalizes earlier information theoretic approaches to coevolution analysis, to become applicable to this problem. It leverages the growing wealth of binding data between PRDs and large numbers of random peptides, and searches for PRD residues that exhibit strong evolutionary covariation with some positions of the statistical profiles of bound peptides. The calculations involve only information from sequences, and thus can be applied to PRDs without crystal structures. We applied the approach to PDZ, SH3 and kinase domains, and evaluated the results using both residue proximity in co-crystal structures and verified binding specificity maps from mutagenesis studies. Discussion Our predictions were found to be strongly correlated with the physical proximity of residues, demonstrating the ability of our approach to detect physical interactions of the binding partners. Some high-scoring pairs were further confirmed to affect binding specificity using previous experimental results. Combining the covariation results also allowed us to predict binding profiles with higher reliability than two other methods that do not explicitly take residue covariation into account. Conclusions The general applicability of our approach to the three different domain families demonstrated in this paper suggests its potential in predicting binding targets and assisting the exploration of binding mechanisms.
Collapse
Affiliation(s)
- Kevin Y Yip
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | | | | | | | | | | | | | | |
Collapse
|
38
|
Fedonin GG, Rakhmaninova AB, Korostelev YD, Laikova ON, Gelfand MS. Machine learning study of DNA binding by transcription factors from the LacI family. Mol Biol 2011. [DOI: 10.1134/s0026893311040054] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
39
|
Bureeva S, Nikolsky Y. Quantitative knowledge-based analysis in compound safety assessment. Expert Opin Drug Metab Toxicol 2011; 7:287-98. [DOI: 10.1517/17425255.2011.553191] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
40
|
Weinberg Z, Breaker RR. R2R--software to speed the depiction of aesthetic consensus RNA secondary structures. BMC Bioinformatics 2011; 12:3. [PMID: 21205310 PMCID: PMC3023696 DOI: 10.1186/1471-2105-12-3] [Citation(s) in RCA: 193] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2010] [Accepted: 01/04/2011] [Indexed: 11/30/2022] Open
Abstract
Background With continuing identification of novel structured noncoding RNAs, there is an increasing need to create schematic diagrams showing the consensus features of these molecules. RNA structural diagrams are typically made either with general-purpose drawing programs like Adobe Illustrator, or with automated or interactive programs specific to RNA. Unfortunately, the use of applications like Illustrator is extremely time consuming, while existing RNA-specific programs produce figures that are useful, but usually not of the same aesthetic quality as those produced at great cost in Illustrator. Additionally, most existing RNA-specific applications are designed for drawing single RNA molecules, not consensus diagrams. Results We created R2R, a computer program that facilitates the generation of aesthetic and readable drawings of RNA consensus diagrams in a fraction of the time required with general-purpose drawing programs. Since the inference of a consensus RNA structure typically requires a multiple-sequence alignment, the R2R user annotates the alignment with commands directing the layout and annotation of the RNA. R2R creates SVG or PDF output that can be imported into Adobe Illustrator, Inkscape or CorelDRAW. R2R can be used to create consensus sequence and secondary structure models for novel RNA structures or to revise models when new representatives for known RNA classes become available. Although R2R does not currently have a graphical user interface, it has proven useful in our efforts to create 100 schematic models of distinct noncoding RNA classes. Conclusions R2R makes it possible to obtain high-quality drawings of the consensus sequence and structural models of many diverse RNA structures with a more practical amount of effort. R2R software is available at http://breaker.research.yale.edu/R2R and as an Additional file.
Collapse
Affiliation(s)
- Zasha Weinberg
- Howard Hughes Medical Institute, Yale University, Box 208103, New Haven, CT 06520-8103, USA.
| | | |
Collapse
|
41
|
Vinogradov SN, Fernández I, Hoogewijs D, Arredondo-Peter R. Phylogenetic relationships of 3/3 and 2/2 hemoglobins in Archaeplastida genomes to bacterial and other eukaryote hemoglobins. MOLECULAR PLANT 2011; 4:42-58. [PMID: 20952597 DOI: 10.1093/mp/ssq040] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
Land plants and algae form a supergroup, the Archaeplastida, believed to be monophyletic. We report the results of an analysis of the phylogeny of putative globins in the currently available genomes to bacterial and other eukaryote hemoglobins (Hbs). Archaeplastida genomes have 3/3 and 2/2 Hbs, with the land plant genomes having group 2 2/2 Hbs, except for the unexpected occurrence of two group 1 2/2 Hbs in Ricinus communis. Bayesian analysis shows that plant 3/3 Hbs are related to vertebrate neuroglobins and bacterial flavohemoglobins (FHbs). We sought to define the bacterial groups, whose ancestors shared the precursors of Archaeplastida Hbs, via Bayesian and neighbor-joining analyses based on COBALT alignment of representative sets of bacterial 3/3 FHb-like globins and group 1 and 2 2/2 Hbs with the corresponding Archaeplastida Hbs. The results suggest that the Archaeplastida 3/3 and group 1 2/2 Hbs could have originated from the horizontal gene transfers (HGTs) that accompanied the two generally accepted endosymbioses of a proteobacterium and a cyanobacterium with a eukaryote ancestor. In contrast, the origin of the group 2 2/2 Hbs unexpectedly appears to involve HGT from a bacterium ancestral to Chloroflexi, Deinococcales, Bacilli, and Actinomycetes. Furthermore, although intron positions and phases are mostly conserved among the land plant 3/3 and 2/2 globin genes, introns are absent in the algal 3/3 genes and intron positions and phases are highly variable in their 2/2 genes. Thus, introns are irrelevant to globin evolution in Archaeplastida.
Collapse
Affiliation(s)
- Serge N Vinogradov
- Department of Biochemistry and Molecular Biology, Wayne State University, School of Medicine, Detroit, MI 48201, USA.
| | | | | | | |
Collapse
|
42
|
Kumar A, Cowen L. Recognition of beta-structural motifs using hidden Markov models trained with simulated evolution. Bioinformatics 2010; 26:i287-93. [PMID: 20529918 PMCID: PMC2881384 DOI: 10.1093/bioinformatics/btq199] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022] Open
Abstract
Motivation: One of the most successful methods to date for recognizing protein sequences that are evolutionarily related, has been profile hidden Markov models. However, these models do not capture pairwise statistical preferences of residues that are hydrogen bonded in β-sheets. We thus explore methods for incorporating pairwise dependencies into these models. Results: We consider the remote homology detection problem for β-structural motifs. In particular, we ask if a statistical model trained on members of only one family in a SCOP β-structural superfamily, can recognize members of other families in that superfamily. We show that HMMs trained with our pairwise model of simulated evolution achieve nearly a median 5% improvement in AUC for β-structural motif recognition as compared to ordinary HMMs. Availability: All datasets and HMMs are available at: http://bcb.cs.tufts.edu/pairwise/ Contact:anoop.kumar@tufts.edu; lenore.cowen@tufts.edu
Collapse
Affiliation(s)
- Anoop Kumar
- Department of Computer Science, Tufts University, Medford, MA, USA.
| | | |
Collapse
|
43
|
Altschul SF, Wootton JC, Zaslavsky E, Yu YK. The construction and use of log-odds substitution scores for multiple sequence alignment. PLoS Comput Biol 2010; 6:e1000852. [PMID: 20657661 PMCID: PMC2904766 DOI: 10.1371/journal.pcbi.1000852] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2009] [Accepted: 06/03/2010] [Indexed: 01/18/2023] Open
Abstract
Most pairwise and multiple sequence alignment programs seek alignments with optimal scores. Central to defining such scores is selecting a set of substitution scores for aligned amino acids or nucleotides. For local pairwise alignment, substitution scores are implicitly of log-odds form. We now extend the log-odds formalism to multiple alignments, using Bayesian methods to construct "BILD" ("Bayesian Integral Log-odds") substitution scores from prior distributions describing columns of related letters. This approach has been used previously only to define scores for aligning individual sequences to sequence profiles, but it has much broader applicability. We describe how to calculate BILD scores efficiently, and illustrate their uses in Gibbs sampling optimization procedures, gapped alignment, and the construction of hidden Markov model profiles. BILD scores enable automated selection of optimal motif and domain model widths, and can inform the decision of whether to include a sequence in a multiple alignment, and the selection of insertion and deletion locations. Other applications include the classification of related sequences into subfamilies, and the definition of profile-profile alignment scores. Although a fully realized multiple alignment program must rely upon more than substitution scores, many existing multiple alignment programs can be modified to employ BILD scores. We illustrate how simple BILD score based strategies can enhance the recognition of DNA binding domains, including the Api-AP2 domain in Toxoplasma gondii and Plasmodium falciparum.
Collapse
Affiliation(s)
- Stephen F Altschul
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland, United States of America.
| | | | | | | |
Collapse
|
44
|
Liu Y, Schmidt B, Maskell DL. MSAProbs: multiple sequence alignment based on pair hidden Markov models and partition function posterior probabilities. Bioinformatics 2010; 26:1958-64. [PMID: 20576627 DOI: 10.1093/bioinformatics/btq338] [Citation(s) in RCA: 188] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
MOTIVATION Multiple sequence alignment is of central importance to bioinformatics and computational biology. Although a large number of algorithms for computing a multiple sequence alignment have been designed, the efficient computation of highly accurate multiple alignments is still a challenge. RESULTS We present MSAProbs, a new and practical multiple alignment algorithm for protein sequences. The design of MSAProbs is based on a combination of pair hidden Markov models and partition functions to calculate posterior probabilities. Furthermore, two critical bioinformatics techniques, namely weighted probabilistic consistency transformation and weighted profile-profile alignment, are incorporated to improve alignment accuracy. Assessed using the popular benchmarks: BAliBASE, PREFAB, SABmark and OXBENCH, MSAProbs achieves statistically significant accuracy improvements over the existing top performing aligners, including ClustalW, MAFFT, MUSCLE, ProbCons and Probalign. Furthermore, MSAProbs is optimized for multi-core CPUs by employing a multi-threaded design, leading to a competitive execution time compared to other aligners. AVAILABILITY The source code of MSAProbs, written in C++, is freely and publicly available from http://msaprobs.sourceforge.net.
Collapse
Affiliation(s)
- Yongchao Liu
- School of Computer Engineering, Nanyang Technological University, Singapore.
| | | | | |
Collapse
|
45
|
Parker AS, Zheng W, Griswold KE, Bailey-Kellogg C. Optimization algorithms for functional deimmunization of therapeutic proteins. BMC Bioinformatics 2010; 11:180. [PMID: 20380721 PMCID: PMC2873530 DOI: 10.1186/1471-2105-11-180] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2009] [Accepted: 04/09/2010] [Indexed: 11/10/2022] Open
Abstract
Background To develop protein therapeutics from exogenous sources, it is necessary to mitigate the risks of eliciting an anti-biotherapeutic immune response. A key aspect of the response is the recognition and surface display by antigen-presenting cells of epitopes, short peptide fragments derived from the foreign protein. Thus, developing minimal-epitope variants represents a powerful approach to deimmunizing protein therapeutics. Critically, mutations selected to reduce immunogenicity must not interfere with the protein's therapeutic activity. Results This paper develops methods to improve the likelihood of simultaneously reducing the anti-biotherapeutic immune response while maintaining therapeutic activity. A dynamic programming approach identifies optimal and near-optimal sets of conservative point mutations to minimize the occurrence of predicted T-cell epitopes in a target protein. In contrast with existing methods, those described here integrate analysis of immunogenicity and stability/activity, are broadly applicable to any protein class, guarantee global optimality, and provide sufficient flexibility for users to limit the total number of mutations and target MHC alleles of interest. The input is simply the primary amino acid sequence of the therapeutic candidate, although crystal structures and protein family sequence alignments may also be input when available. The output is a scored list of sets of point mutations predicted to reduce the protein's immunogenicity while maintaining structure and function. We demonstrate the effectiveness of our approach in a number of case study applications, showing that, in general, our best variants are predicted to be better than those produced by previous deimmunization efforts in terms of either immunogenicity or stability, or both factors. Conclusions By developing global optimization algorithms leveraging well-established immunogenicity and stability prediction techniques, we provide the protein engineer with a mechanism for exploring the favorable sequence space near a targeted protein therapeutic. Our mechanism not only helps identify designs more likely to be effective, but also provides insights into the interrelated implications of design choices.
Collapse
Affiliation(s)
- Andrew S Parker
- Department of Computer Science, Dartmouth College, Hanover, NH 03755, USA
| | | | | | | |
Collapse
|
46
|
Abstract
From a comparatively small number of protein structural domains a staggering array of structural variants has evolved which has, in turn, facilitated an expanse of functional derivatives. Herein I review the primary mechanisms which have contributed to the vastness of our existing, and expanding, protein repertoires.
Collapse
Affiliation(s)
- Roy D Sleator
- Department of Biological Sciences, Cork Institute of Technology.
| |
Collapse
|
47
|
Chalikian TV, Macgregor RB. Origins of Pressure-Induced Protein Transitions. J Mol Biol 2009; 394:834-42. [DOI: 10.1016/j.jmb.2009.10.020] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2009] [Revised: 09/29/2009] [Accepted: 10/12/2009] [Indexed: 10/20/2022]
|
48
|
Gardner PP. The use of covariance models to annotate RNAs in whole genomes. BRIEFINGS IN FUNCTIONAL GENOMICS AND PROTEOMICS 2009; 8:444-50. [PMID: 19833700 DOI: 10.1093/bfgp/elp042] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In this review we discuss bioinformatic issues in non-coding RNA analysis, in particular the annotation of genome sequences using covariance models. Some recent innovations for improving the speed and accuracy of covariance models is discussed.
Collapse
Affiliation(s)
- Paul P Gardner
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton CB10 1SA, UK.
| |
Collapse
|
49
|
Meyer MM, Ames TD, Smith DP, Weinberg Z, Schwalbach MS, Giovannoni SJ, Breaker RR. Identification of candidate structured RNAs in the marine organism 'Candidatus Pelagibacter ubique'. BMC Genomics 2009; 10:268. [PMID: 19531245 PMCID: PMC2704228 DOI: 10.1186/1471-2164-10-268] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2009] [Accepted: 06/16/2009] [Indexed: 02/04/2023] Open
Abstract
Background Metagenomic sequence data are proving to be a vast resource for the discovery of biological components. Yet analysis of this data to identify functional RNAs lags behind efforts to characterize protein diversity. The genome of 'Candidatus Pelagibacter ubique' HTCC 1062 is the closest match for approximately 20% of marine metagenomic sequence reads. It is also small, contains little non-coding DNA, and has strikingly low GC content. Results To aid the discovery of RNA motifs within the marine metagenome we exploited the genomic properties of 'Cand. P. ubique' by targeting our search to long intergenic regions (IGRs) with relatively high GC content. Analysis of known RNAs (rRNA, tRNA, riboswitches etc.) shows that structured RNAs are significantly enriched in such IGRs. To identify additional candidate structured RNAs, we examined other IGRs with similar characteristics from 'Cand. P. ubique' using comparative genomics approaches in conjunction with marine metagenomic data. Employing this strategy, we discovered four candidate structured RNAs including a new riboswitch class as well as three additional likely cis-regulatory elements that precede genes encoding ribosomal proteins S2 and S12, and the cytoplasmic protein component of the signal recognition particle. We also describe four additional potential RNA motifs with few or no examples occurring outside the metagenomic data. Conclusion This work begins the process of identifying functional RNA motifs present in the metagenomic data and illustrates how existing completed genomes may be used to aid in this task.
Collapse
Affiliation(s)
- Michelle M Meyer
- Department of Molecular Cellular and Developmental Biology, Yale University, New Haven, CT 06520, USA.
| | | | | | | | | | | | | |
Collapse
|
50
|
Prediction of deleterious non-synonymous single nucleotide polymorphisms of genes related to ethanol-induced toxicity. Toxicol Lett 2009; 187:99-114. [DOI: 10.1016/j.toxlet.2009.02.007] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2008] [Revised: 02/05/2009] [Accepted: 02/09/2009] [Indexed: 12/30/2022]
|