1
|
Ose NJ, Campitelli P, Modi T, Kazan IC, Kumar S, Ozkan SB. Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics. eLife 2024; 12:RP92063. [PMID: 38713502 PMCID: PMC11076047 DOI: 10.7554/elife.92063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2024] Open
Abstract
We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.
Collapse
Affiliation(s)
- Nicholas James Ose
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Paul Campitelli
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Tushar Modi
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - I Can Kazan
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Sudhir Kumar
- Institute for Genomics and Evolutionary Medicine, Temple UniversityPhiladelphiaUnited States
- Department of Biology, Temple UniversityPhiladelphiaUnited States
- Center for Genomic Medicine Research, King Abdulaziz UniversityJeddahSaudi Arabia
| | - Sefika Banu Ozkan
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| |
Collapse
|
2
|
Khrustalev VV, Khrustaleva TA, Popinako AV. Germline mutations directions are different between introns of the same gene: case study of the gene coding for amyloid-beta precursor protein. Genetica 2023; 151:61-73. [PMID: 36129589 DOI: 10.1007/s10709-022-00166-6] [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: 07/23/2022] [Accepted: 09/08/2022] [Indexed: 02/01/2023]
Abstract
Amyloid-beta precursor protein (APP) is highly conserved in mammals. This feature allowed us to compare nucleotide usage biases in fourfold degenerated sites along the length of its coding region for 146 species of mammals and birds in search of fragments with significant deviations. Even though cytosine usage has the highest value in fourfold degenerated sites in APP coding region from all tested placental mammals, in contrast to marsupial mammals with the bias toward thymine usage, the most frequent germline and somatic mutations in human APP coding region are C to T and G to A transitions. The same mutational AT-pressure is characteristic for germline mutations in introns of human APP gene. However, surprisingly, there are several exceptional introns with deviations in germline mutations rates. The most of those introns surround exons with exceptional biases in nucleotide usage in fourfold degenerated sites. Existence of such fragments in exons 4 and 5, as well as in exon 14, can be connected with the presence of lncRNA genes in complementary strand of DNA. Exceptional nucleotide usage bias in exons 16 and 17 that contain a sequence encoding amyloid-beta peptides can be explained either by the presence of yet unmapped lncRNA(s), or by the autonomous expression of a short mRNA that encodes just C-terminal part of the APP providing an alternative source of amyloid-beta peptides. This hypothesis is supported by the increased rate of T to C transitions in introns 16-17 and 17-18 of Human APP gene relatively to other introns.
Collapse
Affiliation(s)
| | | | - Anna Vladimirovna Popinako
- Bach Institute of Biochemistry, Research Center of Biotechnology of the Russian Academy of Sciences, Moscow, Russian Federation
| |
Collapse
|
3
|
Langille BL, Tierney SM, Bertozzi T, Beasley-Hall PG, Bradford TM, Fagan-Jeffries EP, Hyde J, Leijs R, Richardson M, Saint KM, Stringer DN, Villastrigo A, Humphreys WF, Austin AD, Cooper SJB. Parallel decay of vision genes in subterranean water beetles. Mol Phylogenet Evol 2022; 173:107522. [PMID: 35595008 DOI: 10.1016/j.ympev.2022.107522] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2022] [Revised: 04/19/2022] [Accepted: 04/20/2022] [Indexed: 11/19/2022]
Abstract
In the framework of neutral theory of molecular evolution, genes specific to the development and function of eyes in subterranean animals living in permanent darkness are expected to evolve by relaxed selection, ultimately becoming pseudogenes. However, definitive empirical evidence for the role of neutral processes in the loss of vision over evolutionary time remains controversial. In previous studies, we characterized an assemblage of independently-evolved water beetle (Dytiscidae) species from a subterranean archipelago in Western Australia, where parallel vision and eye loss have occurred. Using a combination of transcriptomics and exon capture, we present evidence of parallel coding sequence decay, resulting from the accumulation of frameshift mutations and premature stop codons, in eight phototransduction genes (arrestins, opsins, ninaC and transient receptor potential channel genes) in 32 subterranean species in contrast to surface species, where these genes have open reading frames. Our results provide strong evidence to support neutral evolutionary processes as a major contributing factor to the loss of phototransduction genes in subterranean animals, with the ultimate fate being the irreversible loss of a light detection system.
Collapse
Affiliation(s)
- Barbara L Langille
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia.
| | - Simon M Tierney
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Hawkesbury Institute for the Environment, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia
| | - Terry Bertozzi
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia
| | - Perry G Beasley-Hall
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia
| | - Tessa M Bradford
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia
| | - Erinn P Fagan-Jeffries
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia
| | - Josephine Hyde
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Western Australia Department of Biodiversity Conservation and Attractions, Kensington, WA 6151, Australia
| | - Remko Leijs
- Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia
| | - Matthew Richardson
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia
| | - Kathleen M Saint
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia
| | - Danielle N Stringer
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia
| | - Adrián Villastrigo
- Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia; Institute of Evolutionary Biology, Passeig Marítim de la Barceloneta, 37-49, 08003, Spain
| | - William F Humphreys
- Western Australian Museum, Locked Bag 40, Welshpool DC, WA 6986, Australia; School of Animal Biology, University of Western Australia, Nedlands, Western Australia, Australia
| | - Andrew D Austin
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia
| | - Steven J B Cooper
- Australian Centre for Evolutionary Biology and Biodiversity, Department of Ecology and Evolution, School of Biological Sciences, University of Adelaide, South Australia 5005, Australia; Evolutionary Biology Unit, South Australian Museum, North Terrace, Adelaide, South Australia 5000, Australia
| |
Collapse
|
4
|
Wang Q, Paskevicius T, Filbert A, Qin W, Kim HJ, Chen XZ, Tang J, Dacks JB, Agellon LB, Michalak M. Phylogenetic and biochemical analysis of calsequestrin structure and association of its variants with cardiac disorders. Sci Rep 2020; 10:18115. [PMID: 33093545 PMCID: PMC7582152 DOI: 10.1038/s41598-020-75097-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Accepted: 10/08/2020] [Indexed: 12/16/2022] Open
Abstract
Calsequestrin is among the most abundant proteins in muscle sarcoplasmic reticulum and displays a high capacity but a low affinity for Ca2+ binding. In mammals, calsequestrin is encoded by two genes, CASQ1 and CASQ2, which are expressed almost exclusively in skeletal and cardiac muscles, respectively. Phylogenetic analysis indicates that calsequestrin is an ancient gene in metazoans, and that the duplication of the ancestral calsequestrin gene took place after the emergence of the lancelet. CASQ2 gene variants associated with catecholaminergic polymorphic ventricular tachycardia (CPVT) in humans are positively correlated with a high degree of evolutionary conservation across all calsequestrin homologues. The mutations are distributed in diverse locations of the calsequestrin protein and impart functional diversity but remarkably manifest in a similar phenotype in humans.
Collapse
Affiliation(s)
- Qian Wang
- Department of Biochemistry, University of Alberta, Edmonton, Alberta, T6G 2H7, Canada
| | - Tautvydas Paskevicius
- Department of Biochemistry, University of Alberta, Edmonton, Alberta, T6G 2H7, Canada
| | - Alexander Filbert
- Division of Infectious Disease, Department of Medicine, University of Alberta, Edmonton, AB, T6G 2G3, Canada
| | - Wenying Qin
- Institute of Biomedical and Pharmaceutical Sciences, Key Laboratory of Fermentation Engineering, Hubei Provincial Cooperative Innovation Center of Industrial Fermentation, Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, Hubei, China
| | - Hyeong Jin Kim
- Department of Biochemistry, University of Alberta, Edmonton, Alberta, T6G 2H7, Canada
| | - Xing-Zhen Chen
- Institute of Biomedical and Pharmaceutical Sciences, Key Laboratory of Fermentation Engineering, Hubei Provincial Cooperative Innovation Center of Industrial Fermentation, Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, Hubei, China.,Department of Physiology, University of Alberta, Edmonton, Alberta, T6G 2H7, Canada
| | - Jingfeng Tang
- Institute of Biomedical and Pharmaceutical Sciences, Key Laboratory of Fermentation Engineering, Hubei Provincial Cooperative Innovation Center of Industrial Fermentation, Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, Hubei, China
| | - Joel B Dacks
- Division of Infectious Disease, Department of Medicine, University of Alberta, Edmonton, AB, T6G 2G3, Canada.
| | - Luis B Agellon
- School of Dietetics and Human Nutrition, McGill University, Ste. Anne de Bellevue, Quebec, H9X 3V9, Canada.
| | - Marek Michalak
- Department of Biochemistry, University of Alberta, Edmonton, Alberta, T6G 2H7, Canada. .,Institute of Biomedical and Pharmaceutical Sciences, Key Laboratory of Fermentation Engineering, Hubei Provincial Cooperative Innovation Center of Industrial Fermentation, Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, Hubei, China.
| |
Collapse
|
5
|
Liu C, Zhao J, Lu W, Dai Y, Hockings J, Zhou Y, Nussinov R, Eng C, Cheng F. Individualized genetic network analysis reveals new therapeutic vulnerabilities in 6,700 cancer genomes. PLoS Comput Biol 2020; 16:e1007701. [PMID: 32101536 PMCID: PMC7062285 DOI: 10.1371/journal.pcbi.1007701] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Revised: 03/09/2020] [Accepted: 01/30/2020] [Indexed: 02/06/2023] Open
Abstract
Tumor-specific genomic alterations allow systematic identification of genetic interactions that promote tumorigenesis and tumor vulnerabilities, offering novel strategies for development of targeted therapies for individual patients. We develop an Individualized Network-based Co-Mutation (INCM) methodology by inspecting over 2.5 million nonsynonymous somatic mutations derived from 6,789 tumor exomes across 14 cancer types from The Cancer Genome Atlas. Our INCM analysis reveals a higher genetic interaction burden on the significantly mutated genes, experimentally validated cancer genes, chromosome regulatory factors, and DNA damage repair genes, as compared to human pan-cancer essential genes identified by CRISPR-Cas9 screenings on 324 cancer cell lines. We find that genes involved in the cancer type-specific genetic subnetworks identified by INCM are significantly enriched in established cancer pathways, and the INCM-inferred putative genetic interactions are correlated with patient survival. By analyzing drug pharmacogenomics profiles from the Genomics of Drug Sensitivity in Cancer database, we show that the network-predicted putative genetic interactions (e.g., BRCA2-TP53) are significantly correlated with sensitivity/resistance of multiple therapeutic agents. We experimentally validated that afatinib has the strongest cytotoxic activity on BT474 (IC50 = 55.5 nM, BRCA2 and TP53 co-mutant) compared to MCF7 (IC50 = 7.7 μM, both BRCA2 and TP53 wild type) and MDA-MB-231 (IC50 = 7.9 μM, BRCA2 wild type but TP53 mutant). Finally, drug-target network analysis reveals several potential druggable genetic interactions by targeting tumor vulnerabilities. This study offers a powerful network-based methodology for identification of candidate therapeutic pathways that target tumor vulnerabilities and prioritization of potential pharmacogenomics biomarkers for development of personalized cancer medicine. Recent efforts to map genetic interactions in tumor cells have suggested that tumor vulnerabilities can be exploited for development of novel targeted therapies. Tumor-specific genomic alterations derived from multi-center cancer genome projects allow identification of genetic interactions that promote tumor vulnerabilities, offering novel strategies for development of targeted cancer therapies. This study develops a novel Individualized Network-based Co-Mutation (termed INCM) methodology for quantifying the putative genetic interactions in cancer. Trained on over 2.5 million nonsynonymous somatic mutations derived from 6,789 tumor exomes across 14 cancer type, we found that genes identified in the cancer type-specific genetic subnetworks were significantly enriched in established cancer pathways. The network-predicted putative genetic interactions are correlated with patient survival. By analyzing drug pharmacogenomics profiles, we showed that the network-predicted putative genetic interactions (e.g., BRCA2-TP53) were significantly correlated with sensitivity/resistance of anticancer drugs (e.g., afatinib) and we experimentally validated it in breast cancer cell lines. Finally, drug-target network analysis reveals several potential druggable genetic interactions (e.g., PIK3CA-PTEN) by targeting tumor vulnerabilities. This study offers a generalizable network-based approach for comprehensive identification of candidate therapeutic pathways that target tumor vulnerabilities and prioritization of potential prognostic and pharmacogenomics biomarkers for development of personalized cancer medicine.
Collapse
Affiliation(s)
- Chuang Liu
- Alibaba Research Center for Complexity Sciences, Hangzhou Normal University, Hangzhou, Zhejiang, China
| | - Junfei Zhao
- Department of Systems Biology, Columbia University, New York, New York, United States of America
- Department of Biomedical Informatics, Columbia University, New York, New York, United States of America
| | - Weiqiang Lu
- Shanghai Key Laboratory of Regulatory Biology, Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, China
| | - Yao Dai
- Alibaba Research Center for Complexity Sciences, Hangzhou Normal University, Hangzhou, Zhejiang, China
| | - Jennifer Hockings
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America
| | - Yadi Zhou
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America
| | - Ruth Nussinov
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, Maryland, United States of America
- Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Charis Eng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America
- Department of Molecular Medicine, Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, Ohio, United States of America
- Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, Ohio, United States of America
- Taussig Cancer Institute, Cleveland Clinic, Cleveland, Ohio, United States of America
- Department of Genetics and Genome Sciences, Case Western Reserve University School of Medicine, Cleveland, Ohio, United States of America
| | - Feixiong Cheng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America
- Department of Molecular Medicine, Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, Ohio, United States of America
- Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, Ohio, United States of America
- * E-mail:
| |
Collapse
|
6
|
Patel R, Kumar S. On estimating evolutionary probabilities of population variants. BMC Evol Biol 2019; 19:133. [PMID: 31238981 PMCID: PMC6593550 DOI: 10.1186/s12862-019-1455-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2018] [Accepted: 06/06/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The evolutionary probability (EP) of an allele in a DNA or protein sequence predicts evolutionarily permissible (ePerm; EP ≥ 0.05) and forbidden (eForb; EP < 0.05) variants. EP of an allele represents an independent evolutionary expectation of observing an allele in a population based solely on the long-term substitution patterns captured in a multiple sequence alignment. In the neutral theory, EP and population frequencies can be compared to identify neutral and non-neutral alleles. This approach has been used to discover candidate adaptive polymorphisms in humans, which are eForbs segregating with high frequencies. The original method to compute EP requires the evolutionary relationships and divergence times of species in the sequence alignment (a timetree), which are not known with certainty for most datasets. This requirement impedes a general use of the original EP formulation. Here, we present an approach in which the phylogeny and times are inferred from the sequence alignment itself prior to the EP calculation. We evaluate if the modified EP approach produces results that are similar to those from the original method. RESULTS We compared EP estimates from the original and the modified approaches by using more than 18,000 protein sequence alignments containing orthologous sequences from 46 vertebrate species. For the original EP calculations, we used species relationships from UCSC and divergence times from TimeTree web resource, and the resulting EP estimates were considered to be the ground truth. We found that the modified approaches produced reasonable EP estimates for HGMD disease missense variant and 1000 Genomes Project missense variant datasets. Our results showed that reliable estimates of EP can be obtained without a priori knowledge of the sequence phylogeny and divergence times. We also found that, in order to obtain robust EP estimates, it is important to assemble a dataset with many sequences, sampling from a diversity of species groups. CONCLUSION We conclude that the modified EP approach will be generally applicable for alignments and enable the detection of potentially neutral, deleterious, and adaptive alleles in populations.
Collapse
Affiliation(s)
- Ravi Patel
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, PA, 19122, USA.,Department of Biology, Temple University, Philadelphia, PA, 19122, USA
| | - Sudhir Kumar
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, PA, 19122, USA. .,Department of Biology, Temple University, Philadelphia, PA, 19122, USA. .,Center for Excellence in Genome Medicine and Research, King Abdulaziz University, Jeddah, Saudi Arabia.
| |
Collapse
|