1
|
Kumar A K, Rathore RS. Categorization of hotspots into three types - weak, moderate and strong to distinguish protein-protein versus protein-peptide interactions. J Biomol Struct Dyn 2024; 42:9348-9360. [PMID: 37649387 DOI: 10.1080/07391102.2023.2252077] [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/15/2023] [Accepted: 08/18/2023] [Indexed: 09/01/2023]
Abstract
Protein-protein and protein-peptide interactions (PPI and PPepI) belong to a similar category of interactions, yet seemingly subtle differences exist among them. To characterize differences between protein-protein (PP) and protein-peptide (PPep) interactions, we have focussed on two important classes of residues-hotspot and anchor residues. Using implicit solvation-based free energy calculations, a very large-scale alanine scanning has been performed on benchmark datasets, consisting of over 5700 interface residues. The differences in the two categories are more pronounced, if the data were divided into three distinct types, namely - weak hotspots (having binding free energy loss upon Ala mutation, ΔΔG, ∼2-10 kcal/mol), moderate hotspots (ΔΔG, ∼10-20 kcal/mol) and strong hotspots (ΔΔG ≥ ∼20 kcal/mol). The analysis suggests that for PPI, weak hotspots are predominantly populated by polar and hydrophobic residues. The distribution shifts towards charged and polar residues for moderate hotspot and charged residues (principally Arg) are overwhelmingly present in the strong hotspot. On the other hand, in the PPepI dataset, the distribution shifts from predominantly hydrophobic and polar (in the weak type) to almost similar preference for polar, hydrophobic and charged residues (in moderate type) and finally the charged residue (Arg) and Trp are mostly occupied in the strong type. The preferred anchor residues in both categories are Arg, Tyr and Leu, possessing bulky side chain and which also strike a delicate balance between side chain flexibility and rigidity. The present knowledge should aid in effective design of biologics, by augmentation or disruption of PPIs with peptides or peptidomimetics.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Kiran Kumar A
- Department of Bioinformatics, School of Earth, Biological and Environmental Sciences, Central University of South Bihar, Gaya, India
| | - R S Rathore
- Department of Bioinformatics, School of Earth, Biological and Environmental Sciences, Central University of South Bihar, Gaya, India
| |
Collapse
|
2
|
Seiler J, Beye M. Honeybees' novel complementary sex-determining system: function and origin. Trends Genet 2024:S0168-9525(24)00178-1. [PMID: 39232877 DOI: 10.1016/j.tig.2024.08.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2024] [Revised: 08/06/2024] [Accepted: 08/07/2024] [Indexed: 09/06/2024]
Abstract
Complementary sex determination regulates female and male development in honeybees (Apis mellifera) via heterozygous versus homo-/hemizygous genotypes of the csd (complementary sex determiner) gene involving numerous naturally occurring alleles. This lineage-specific function offers a rare opportunity to understand an undescribed regulatory mechanism and the molecular evolutionary path leading to this mechanism. We reviewed recent advances in understanding how Csd recognizes different versus identical protein variants, how these variants regulate downstream pathways and sexual differentiation, and how this mechanism has evolved and been shaped by evolutionary forces. Finally, we highlighted the shared regulatory principles of sex determination despite the diversity of primary signals and demonstrated that lineage-specific mutations are very informative for characterizing newly evolved functions.
Collapse
Affiliation(s)
- Jana Seiler
- Institute of Evolutionary Genetics, Heinrich-Heine University, Düsseldorf, Germany
| | - Martin Beye
- Institute of Evolutionary Genetics, Heinrich-Heine University, Düsseldorf, Germany.
| |
Collapse
|
3
|
Venezian J, Bar-Yosef H, Ben-Arie Zilberman H, Cohen N, Kleifeld O, Fernandez-Recio J, Glaser F, Shiber A. Diverging co-translational protein complex assembly pathways are governed by interface energy distribution. Nat Commun 2024; 15:2638. [PMID: 38528060 DOI: 10.1038/s41467-024-46881-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 03/12/2024] [Indexed: 03/27/2024] Open
Abstract
Protein-protein interactions are at the heart of all cellular processes, with the ribosome emerging as a platform, orchestrating the nascent-chain interplay dynamics. Here, to study the characteristics governing co-translational protein folding and complex assembly, we combine selective ribosome profiling, imaging, and N-terminomics with all-atoms molecular dynamics. Focusing on conserved N-terminal acetyltransferases (NATs), we uncover diverging co-translational assembly pathways, where highly homologous subunits serve opposite functions. We find that only a few residues serve as "hotspots," initiating co-translational assembly interactions upon exposure at the ribosome exit tunnel. These hotspots are characterized by high binding energy, anchoring the entire interface assembly. Alpha-helices harboring hotspots are highly thermolabile, folding and unfolding during simulations, depending on their partner subunit to avoid misfolding. In vivo hotspot mutations disrupted co-translational complexation, leading to aggregation. Accordingly, conservation analysis reveals that missense NATs variants, causing neurodevelopmental and neurodegenerative diseases, disrupt putative hotspot clusters. Expanding our study to include phosphofructokinase, anthranilate synthase, and nucleoporin subcomplex, we employ AlphaFold-Multimer to model the complexes' complete structures. Computing MD-derived interface energy profiles, we find similar trends. Here, we propose a model based on the distribution of interface energy as a strong predictor of co-translational assembly.
Collapse
Affiliation(s)
- Johannes Venezian
- Faculty of Biology, Technion Israel institute of Technology, Haifa, Israel
| | - Hagit Bar-Yosef
- Faculty of Biology, Technion Israel institute of Technology, Haifa, Israel
| | | | - Noam Cohen
- Faculty of Biology, Technion Israel institute of Technology, Haifa, Israel
| | - Oded Kleifeld
- Faculty of Biology, Technion Israel institute of Technology, Haifa, Israel
| | - Juan Fernandez-Recio
- Instituto de Ciencias de la Vid y del Vino (ICVV), CSIC-Universidad de La Rioja-Gobierno de La Rioja, Logroño, Spain
| | - Fabian Glaser
- Lorry I. Lokey Interdisciplinary Center for Life Sciences & Engineering, Haifa, Israel
| | - Ayala Shiber
- Faculty of Biology, Technion Israel institute of Technology, Haifa, Israel.
| |
Collapse
|
4
|
Pushan SS, Samantaray M, Rajagopalan M, Ramaswamy A. Structural dynamics of influenza A (H1N1) hemagglutinin protein: a comparative study of Indian (2018) isolate with its evolutionary neighbor, Californian (2009) strain. J Biomol Struct Dyn 2024:1-14. [PMID: 38379377 DOI: 10.1080/07391102.2024.2317985] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2023] [Accepted: 02/07/2024] [Indexed: 02/22/2024]
Abstract
This work highlights the structure and dynamics of two trimeric HA proteins of the H1N1 virus from different origins, the pandemic Californian (HACal) and its closest Indian neighbor (HAInd), reported in 2009 and 2018, respectively. Because of mutation, HAInd acquires new N-glycosylation and epitope binding sites along with mutations at RBD, which might trigger an altered viral-host interaction mechanism. Molecular dynamics simulations performed on HA trimers for a period of 250 ns reveal the highly dynamic nature of HACal trimers inherited by the flexibility of HA monomers. In the trimer, the dynamics of one monomer are more pronounced compared to others, and the enhanced dynamics of RBD especially gain attention as they plays a key role during fusion. Conversely, the mutant HAInd trimer effectively establishes more H-bond interactions, and accordingly, the trimer undergoes more stabilized dynamics with a relatively lower amplitude of RBD dynamics, as endorsed by the reduced RMSD, Rg, and SASA variations. The cooperative and anti-cooperative motions dissected for the subdomains of both strains also reveal a prominent anticorrelative motion of RBD against other subdomains. In agreement, the free energy landscape of stable HAInd is also characterized by a single lowest wide energy basin instead of the two minimum energy basins of the HACal trimer. In essence, the mutant HAInd acquires a highly stable conformation with novel functional features, which calls for (i) further investigation on the emerging mutation-mediated variation in viral-host binding mechanism and (ii) the need for further design of site-specific potential inhibitors to face future challenges.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Shilpa Sri Pushan
- Department of Bioinformatics, Pondicherry University (A Central University), Kalapet, Puducherry, India
| | - Mahesh Samantaray
- Department of Bioinformatics, Pondicherry University (A Central University), Kalapet, Puducherry, India
| | - Muthukumaran Rajagopalan
- Department of Biological Sciences and Bioengineering, Indian Institute of Technology Kanpur, Kanpur
| | - Amutha Ramaswamy
- Department of Bioinformatics, Pondicherry University (A Central University), Kalapet, Puducherry, India
| |
Collapse
|
5
|
Guo XY, Yi L, Yang J, An HW, Yang ZX, Wang H. Self-assembly of peptide nanomaterials at biointerfaces: molecular design and biomedical applications. Chem Commun (Camb) 2024; 60:2009-2021. [PMID: 38275083 DOI: 10.1039/d3cc05811e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2024]
Abstract
Self-assembly is an important strategy for constructing ordered structures and complex functions in nature. Based on this, people can imitate nature and artificially construct functional materials with novel structures through the supermolecular self-assembly pathway of biological interfaces. Among the many assembly units, peptide molecular self-assembly has received widespread attention in recent years. In this review, we introduce the interactions (hydrophobic interaction, hydrogen bond, and electrostatic interaction) between peptide nanomaterials and biological interfaces, summarizing the latest advancements in multifunctional self-assembling peptide materials. We systematically demonstrate the assembly mechanisms of peptides at biological interfaces, such as proteins and cell membranes, while highlighting their application potential and challenges in fields like drug delivery, antibacterial strategies, and cancer therapy.
Collapse
Affiliation(s)
- Xin-Yuan Guo
- College of Chemistry, Huazhong Agricultural University, Shizishan 1, Hongshan District, Wuhan, 430070, China
- CAS Key Laboratory for Biomedical Effects of Nanomaterials and Nanosafety, CAS Center for Excellence in Nanoscience, National Center for Nanoscience and Technology (NCNST), Beijing, 100190, China.
| | - Li Yi
- CAS Key Laboratory for Biomedical Effects of Nanomaterials and Nanosafety, CAS Center for Excellence in Nanoscience, National Center for Nanoscience and Technology (NCNST), Beijing, 100190, China.
| | - Jia Yang
- CAS Key Laboratory for Biomedical Effects of Nanomaterials and Nanosafety, CAS Center for Excellence in Nanoscience, National Center for Nanoscience and Technology (NCNST), Beijing, 100190, China.
| | - Hong-Wei An
- CAS Key Laboratory for Biomedical Effects of Nanomaterials and Nanosafety, CAS Center for Excellence in Nanoscience, National Center for Nanoscience and Technology (NCNST), Beijing, 100190, China.
| | - Zi-Xin Yang
- College of Chemistry, Huazhong Agricultural University, Shizishan 1, Hongshan District, Wuhan, 430070, China
| | - Hao Wang
- CAS Key Laboratory for Biomedical Effects of Nanomaterials and Nanosafety, CAS Center for Excellence in Nanoscience, National Center for Nanoscience and Technology (NCNST), Beijing, 100190, China.
| |
Collapse
|
6
|
Parisi G, Piacentini R, Incocciati A, Bonamore A, Macone A, Rupert J, Zacco E, Miotto M, Milanetti E, Tartaglia GG, Ruocco G, Boffi A, Di Rienzo L. Design of protein-binding peptides with controlled binding affinity: the case of SARS-CoV-2 receptor binding domain and angiotensin-converting enzyme 2 derived peptides. Front Mol Biosci 2024; 10:1332359. [PMID: 38250735 PMCID: PMC10797010 DOI: 10.3389/fmolb.2023.1332359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2023] [Accepted: 12/14/2023] [Indexed: 01/23/2024] Open
Abstract
The development of methods able to modulate the binding affinity between proteins and peptides is of paramount biotechnological interest in view of a vast range of applications that imply designed polypeptides capable to impair or favour Protein-Protein Interactions. Here, we applied a peptide design algorithm based on shape complementarity optimization and electrostatic compatibility and provided the first experimental in vitro proof of the efficacy of the design algorithm. Focusing on the interaction between the SARS-CoV-2 Spike Receptor-Binding Domain (RBD) and the human angiotensin-converting enzyme 2 (ACE2) receptor, we extracted a 23-residues long peptide that structurally mimics the major interacting portion of the ACE2 receptor and designed in silico five mutants of such a peptide with a modulated affinity. Remarkably, experimental KD measurements, conducted using biolayer interferometry, matched the in silico predictions. Moreover, we investigated the molecular determinants that govern the variation in binding affinity through molecular dynamics simulation, by identifying the mechanisms driving the different values of binding affinity at a single residue level. Finally, the peptide sequence with the highest affinity, in comparison with the wild type peptide, was expressed as a fusion protein with human H ferritin (HFt) 24-mer. Solution measurements performed on the latter constructs confirmed that peptides still exhibited the expected trend, thereby enhancing their efficacy in RBD binding. Altogether, these results indicate the high potentiality of this general method in developing potent high-affinity vectors for hindering/enhancing protein-protein associations.
Collapse
Affiliation(s)
- Giacomo Parisi
- Department of Basic and Applied Sciences for Engineering (SBAI), Università“Sapienza”, Roma, Italy
| | - Roberta Piacentini
- Department of Biochemical Sciences “Alessandro Rossi Fanelli”, Università“Sapienza”, Roma, Italy
| | - Alessio Incocciati
- Department of Biochemical Sciences “Alessandro Rossi Fanelli”, Università“Sapienza”, Roma, Italy
| | - Alessandra Bonamore
- Department of Biochemical Sciences “Alessandro Rossi Fanelli”, Università“Sapienza”, Roma, Italy
| | - Alberto Macone
- Department of Biochemical Sciences “Alessandro Rossi Fanelli”, Università“Sapienza”, Roma, Italy
| | - Jakob Rupert
- Department of Biology and Biotechnologies “Charles Darwin”, Università“Sapienza”, Roma, Italy
- Centre for Human Technologies (CHT), Istituto Italiano di Tecnologia (IIT), Genova, Italy
| | - Elsa Zacco
- Centre for Human Technologies (CHT), Istituto Italiano di Tecnologia (IIT), Genova, Italy
| | - Mattia Miotto
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Roma, Italy
| | - Edoardo Milanetti
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Roma, Italy
- Department of Physics, Università“Sapienza”, Roma, Italy
| | - Gian Gaetano Tartaglia
- Department of Biology and Biotechnologies “Charles Darwin”, Università“Sapienza”, Roma, Italy
- Centre for Human Technologies (CHT), Istituto Italiano di Tecnologia (IIT), Genova, Italy
| | - Giancarlo Ruocco
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Roma, Italy
- Department of Physics, Università“Sapienza”, Roma, Italy
| | - Alberto Boffi
- Department of Biochemical Sciences “Alessandro Rossi Fanelli”, Università“Sapienza”, Roma, Italy
| | - Lorenzo Di Rienzo
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Roma, Italy
| |
Collapse
|
7
|
Rana MM, Nguyen DD. Geometric Graph Learning to Predict Changes in Binding Free Energy and Protein Thermodynamic Stability upon Mutation. J Phys Chem Lett 2023; 14:10870-10879. [PMID: 38032742 DOI: 10.1021/acs.jpclett.3c02679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2023]
Abstract
Accurate prediction of binding free energy changes upon mutations is vital for optimizing drugs, designing proteins, understanding genetic diseases, and cost-effective virtual screening. While machine learning methods show promise in this domain, achieving accuracy and generalization across diverse data sets remains a challenge. This study introduces Geometric Graph Learning for Protein-Protein Interactions (GGL-PPI), a novel approach integrating geometric graph representation and machine learning to forecast mutation-induced binding free energy changes. GGL-PPI leverages atom-level graph coloring and multiscale weighted colored geometric subgraphs to capture structural features of biomolecules, demonstrating superior performance on three standard data sets, namely, AB-Bind, SKEMPI 1.0, and SKEMPI 2.0 data sets. The model's efficacy extends to predicting protein thermodynamic stability in a blind test set, providing unbiased predictions for both direct and reverse mutations and showcasing notable generalization. GGL-PPI's precision in predicting changes in binding free energy and stability due to mutations enhances our comprehension of protein complexes, offering valuable insights for drug design endeavors.
Collapse
Affiliation(s)
- Md Masud Rana
- Department of Mathematics, University of Kentucky, Lexington, Kentucky 40506, United States
| | - Duc Duy Nguyen
- Department of Mathematics, University of Kentucky, Lexington, Kentucky 40506, United States
| |
Collapse
|
8
|
Otte M, Netschitailo O, Weidtkamp-Peters S, Seidel CA, Beye M. Recognition of polymorphic Csd proteins determines sex in the honeybee. SCIENCE ADVANCES 2023; 9:eadg4239. [PMID: 37792946 PMCID: PMC10550236 DOI: 10.1126/sciadv.adg4239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 09/05/2023] [Indexed: 10/06/2023]
Abstract
Sex in honeybees, Apis mellifera, is genetically determined by heterozygous versus homo/hemizygous genotypes involving numerous alleles at the single complementary sex determination locus. The molecular mechanism of sex determination is however unknown because there are more than 4950 known possible allele combinations, but only two sexes in the species. We show how protein variants expressed from complementary sex determiner (csd) gene determine sex. In females, the amino acid differences between Csd variants at the potential-specifying domain (PSD) direct the selection of a conserved coiled-coil domain for binding and protein complexation. This recognition mechanism activates Csd proteins and, thus, the female pathway. In males, the absence of polymorphisms establishes other binding elements at PSD for binding and complexation of identical Csd proteins. This second recognition mechanism inactivates Csd proteins and commits male development via default pathway. Our results demonstrate that the recognition of different versus identical variants of a single protein is a mechanism to determine sex.
Collapse
Affiliation(s)
- Marianne Otte
- Institute of Evolutionary Genetics, Heinrich-Heine University, Düsseldorf, Germany
| | - Oksana Netschitailo
- Institute of Evolutionary Genetics, Heinrich-Heine University, Düsseldorf, Germany
| | | | - Claus A. M. Seidel
- Institut für Physikalische Chemie, Heinrich-Heine University, Düsseldorf, Germany
| | - Martin Beye
- Institute of Evolutionary Genetics, Heinrich-Heine University, Düsseldorf, Germany
| |
Collapse
|
9
|
Samantaray M, Pushan SS, Rajagopalan M, Ramaswamy A. Structural dynamics of the RNA dependent RNA polymerase of H1N1 strain affecting humans: a bioinformatics approach. J Biomol Struct Dyn 2023:1-14. [PMID: 37728538 DOI: 10.1080/07391102.2023.2259481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Accepted: 09/11/2023] [Indexed: 09/21/2023]
Abstract
The Influenza flu is a pandemic disease that renders the highest risk factor to the society due to its efficient ability of airborne transmission. Studies on the H1N1 strain gained significant focus, since its pandemic outbreak in 2009 and particularly the computational studies on its structural elements significantly aided in revealing their functional uniqueness. Among the 10 structural proteins of H1N1, the RNA-dependent RNA polymerase (RdRp) heterotrimeric protein complex, which is responsible for the synthesis of viral RNA (vRNA) from the negative-sense RNA genome of the virus, is the focus of the present study. This study aimed to investigate the structural dynamics of the RdRp complex with particular emphasis on the reported 17 mutations. The mutant strain is more stabilized by strong concerted residue-residue interactions at both intra- and inter- monomeric levels. In comparison, the mutant strain is structurally flexible with enhanced stabilizing interactions. The structural dynamics of RdRp are significantly governed by the dynamics of the (i) endonuclease domain of PA, (ii) RNA-entry region of PB1 and (iii) cap-binding region of PB2. Explicitly, the cap binding region of PB2 expresses (i) a concerted motion with the RNA-entry region, along with (ii) an anti-correlated motion with the endonuclease domain of the PA subunit, which further supports the stable dynamics of cap-binding towards RNA binding. These findings contribute to the understanding of the structural dynamics associated with the pandemic and mutant structures of RdRp and render a basic knowledge for further development of novel inhibitors towards influenza flu affected humans.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Mahesh Samantaray
- Department of Bioinformatics, School of Life Sciences, Pondicherry University, R. V. Nagar, Puducherry, India
| | - Shilpa Sri Pushan
- Department of Bioinformatics, School of Life Sciences, Pondicherry University, R. V. Nagar, Puducherry, India
| | - Muthukumaran Rajagopalan
- Department of Biological Sciences and Bioengineering, Indian Institute of Technology Kanpur, Kanpur, India
| | - Amutha Ramaswamy
- Department of Bioinformatics, School of Life Sciences, Pondicherry University, R. V. Nagar, Puducherry, India
| |
Collapse
|
10
|
Li J, Kang G, Wang J, Yuan H, Wu Y, Meng S, Wang P, Zhang M, Wang Y, Feng Y, Huang H, de Marco A. Affinity maturation of antibody fragments: A review encompassing the development from random approaches to computational rational optimization. Int J Biol Macromol 2023; 247:125733. [PMID: 37423452 DOI: 10.1016/j.ijbiomac.2023.125733] [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/03/2023] [Revised: 07/04/2023] [Accepted: 07/06/2023] [Indexed: 07/11/2023]
Abstract
Routinely screened antibody fragments usually require further in vitro maturation to achieve the desired biophysical properties. Blind in vitro strategies can produce improved ligands by introducing random mutations into the original sequences and selecting the resulting clones under more and more stringent conditions. Rational approaches exploit an alternative perspective that aims first at identifying the specific residues potentially involved in the control of biophysical mechanisms, such as affinity or stability, and then to evaluate what mutations could improve those characteristics. The understanding of the antigen-antibody interactions is instrumental to develop this process the reliability of which, consequently, strongly depends on the quality and completeness of the structural information. Recently, methods based on deep learning approaches critically improved the speed and accuracy of model building and are promising tools for accelerating the docking step. Here, we review the features of the available bioinformatic instruments and analyze the reports illustrating the result obtained with their application to optimize antibody fragments, and nanobodies in particular. Finally, the emerging trends and open questions are summarized.
Collapse
Affiliation(s)
- Jiaqi Li
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Frontiers Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin 300072, China
| | - Guangbo Kang
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Frontiers Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin 300072, China
| | - Jiewen Wang
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Frontiers Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin 300072, China
| | - Haibin Yuan
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Frontiers Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin 300072, China
| | - Yili Wu
- Zhejiang Provincial Clinical Research Center for Mental Disorders, School of Mental Health and the Affiliated Kangning Hospital, Institute of Aging, Key Laboratory of Alzheimer's Disease of Zhejiang Province, Wenzhou Medical University, Oujiang Laboratory, Wenzhou, Zhejiang 325035, China
| | - Shuxian Meng
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China
| | - Ping Wang
- New Technology R&D Department, Tianjin Modern Innovative TCM Technology Company Limited, Tianjin 300392, China
| | - Miao Zhang
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; China Resources Biopharmaceutical Company Limited, Beijing 100029, China
| | - Yuli Wang
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Tianjin Pharmaceutical Da Ren Tang Group Corporation Limited, Traditional Chinese Pharmacy Research Institute, Tianjin Key Laboratory of Quality Control in Chinese Medicine, Tianjin 300457, China; State Key Laboratory of Drug Delivery Technology and Pharmacokinetics, Tianjin Institute of Pharmaceutical Research, Tianjin 300193, China
| | - Yuanhang Feng
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China
| | - He Huang
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Frontiers Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin 300072, China.
| | - Ario de Marco
- Laboratory for Environmental and Life Sciences, University of Nova Gorica, Nova Gorica, Slovenia.
| |
Collapse
|
11
|
Hall MWJ, Shorthouse D, Alcraft R, Jones PH, Hall BA. Mutations observed in somatic evolution reveal underlying gene mechanisms. Commun Biol 2023; 6:753. [PMID: 37468606 PMCID: PMC10356810 DOI: 10.1038/s42003-023-05136-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Accepted: 07/11/2023] [Indexed: 07/21/2023] Open
Abstract
Highly sensitive DNA sequencing techniques have allowed the discovery of large numbers of somatic mutations in normal tissues. Some mutations confer a competitive advantage over wild-type cells, generating expanding clones that spread through the tissue. Competition between mutant clones leads to selection. This process can be considered a large scale, in vivo screen for mutations increasing cell fitness. It follows that somatic missense mutations may offer new insights into the relationship between protein structure, function and cell fitness. We present a flexible statistical method for exploring the selection of structural features in data sets of somatic mutants. We show how this approach can evidence selection of specific structural features in key drivers in aged tissues. Finally, we show how drivers may be classified as fitness-enhancing and fitness-suppressing through different patterns of mutation enrichment. This method offers a route to understanding the mechanism of protein function through in vivo mutant selection.
Collapse
Affiliation(s)
| | - David Shorthouse
- Department of Medical Physics and Biomedical Engineering, Malet Place Engineering Building, University College London, Gower Street, London, WC1E 6BT, UK
| | - Rachel Alcraft
- Advanced Research Computing, University College London, London, UK
| | - Philip H Jones
- Wellcome Sanger Institute, Hinxton, CB10 1SA, UK
- Department of Oncology, University of Cambridge, Cambridge, CB2 0XZ, UK
| | - Benjamin A Hall
- Department of Medical Physics and Biomedical Engineering, Malet Place Engineering Building, University College London, Gower Street, London, WC1E 6BT, UK.
| |
Collapse
|
12
|
Grassmann G, Di Rienzo L, Gosti G, Leonetti M, Ruocco G, Miotto M, Milanetti E. Electrostatic complementarity at the interface drives transient protein-protein interactions. Sci Rep 2023; 13:10207. [PMID: 37353566 PMCID: PMC10290103 DOI: 10.1038/s41598-023-37130-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Accepted: 06/16/2023] [Indexed: 06/25/2023] Open
Abstract
Understanding the mechanisms driving bio-molecules binding and determining the resulting complexes' stability is fundamental for the prediction of binding regions, which is the starting point for drug-ability and design. Characteristics like the preferentially hydrophobic composition of the binding interfaces, the role of van der Waals interactions, and the consequent shape complementarity between the interacting molecular surfaces are well established. However, no consensus has yet been reached on the role of electrostatic. Here, we perform extensive analyses on a large dataset of protein complexes for which both experimental binding affinity and pH data were available. Probing the amino acid composition, the disposition of the charges, and the electrostatic potential they generated on the protein molecular surfaces, we found that (i) although different classes of dimers do not present marked differences in the amino acid composition and charges disposition in the binding region, (ii) homodimers with identical binding region show higher electrostatic compatibility with respect to both homodimers with non-identical binding region and heterodimers. Interestingly, (iii) shape and electrostatic complementarity, for patches defined on short-range interactions, behave oppositely when one stratifies the complexes by their binding affinity: complexes with higher binding affinity present high values of shape complementarity (the role of the Lennard-Jones potential predominates) while electrostatic tends to be randomly distributed. Conversely, complexes with low values of binding affinity exploit Coulombic complementarity to acquire specificity, suggesting that electrostatic complementarity may play a greater role in transient (or less stable) complexes. In light of these results, (iv) we provide a novel, fast, and efficient method, based on the 2D Zernike polynomial formalism, to measure electrostatic complementarity without the need of knowing the complex structure. Expanding the electrostatic potential on a basis of 2D orthogonal polynomials, we can discriminate between transient and permanent protein complexes with an AUC of the ROC of [Formula: see text] 0.8. Ultimately, our work helps shedding light on the non-trivial relationship between the hydrophobic and electrostatic contributions in the binding interfaces, thus favoring the development of new predictive methods for binding affinity characterization.
Collapse
Affiliation(s)
- Greta Grassmann
- Department of Biochemical Sciences "Alessandro Rossi Fanelli", Sapienza University of Rome, Piazzale Aldo Moro 5, 00185, Rome, Italy
- Center for Life Nano & Neuro Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161, Rome, Italy
| | - Lorenzo Di Rienzo
- Center for Life Nano & Neuro Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161, Rome, Italy
| | - Giorgio Gosti
- Center for Life Nano & Neuro Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161, Rome, Italy
- Soft and Living Matter Laboratory, Institute of Nanotechnology, Consiglio Nazionale delle Ricerche, 00185, Rome, Italy
| | - Marco Leonetti
- Center for Life Nano & Neuro Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161, Rome, Italy
- Soft and Living Matter Laboratory, Institute of Nanotechnology, Consiglio Nazionale delle Ricerche, 00185, Rome, Italy
| | - Giancarlo Ruocco
- Center for Life Nano & Neuro Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161, Rome, Italy
- Department of Physics, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185, Rome, Italy
| | - Mattia Miotto
- Center for Life Nano & Neuro Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161, Rome, Italy.
| | - Edoardo Milanetti
- Center for Life Nano & Neuro Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161, Rome, Italy.
- Department of Physics, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185, Rome, Italy.
| |
Collapse
|
13
|
Kiemel K, Gurke M, Paraskevopoulou S, Havenstein K, Weithoff G, Tiedemann R. Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species. Sci Rep 2022; 12:22626. [PMID: 36587065 PMCID: PMC9805463 DOI: 10.1038/s41598-022-27137-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2022] [Accepted: 12/26/2022] [Indexed: 01/01/2023] Open
Abstract
Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.
Collapse
Affiliation(s)
- K. Kiemel
- grid.11348.3f0000 0001 0942 1117Unit of Evolutionary Biology/Systematic Zoology, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht Straße 24-25, 14476 Potsdam, Germany
| | - M. Gurke
- grid.422371.10000 0001 2293 9957Museum für Naturkunde – Leibniz Institute for Evolution and Biodiversity Science, Invalidenstraße 43, 10115 Berlin, Germany ,grid.7468.d0000 0001 2248 7639Department of Biology, Humboldt-University, Invalidenstraße 42, 10115 Berlin, Germany
| | - S. Paraskevopoulou
- grid.4514.40000 0001 0930 2361Department of Biology, Lund University, Microbiology Group, Sölvegatan 35, 223 62 Lund, Sweden
| | - K. Havenstein
- grid.11348.3f0000 0001 0942 1117Unit of Evolutionary Biology/Systematic Zoology, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht Straße 24-25, 14476 Potsdam, Germany
| | - G. Weithoff
- grid.11348.3f0000 0001 0942 1117Unit of Ecology and Ecosystem Modelling, Institute for Biochemistry and Biology, University of Potsdam, Am Neuen Palais 10, 14469 Potsdam, Germany
| | - R. Tiedemann
- grid.11348.3f0000 0001 0942 1117Unit of Evolutionary Biology/Systematic Zoology, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht Straße 24-25, 14476 Potsdam, Germany
| |
Collapse
|
14
|
Tiwari S, Pandey VP, Yadav K, Dwivedi UN. Modulation of interaction of BRCA1-RAD51 and BRCA1-AURKA protein complexes by natural metabolites using as possible therapeutic intervention toward cardiotoxic effects of cancer drugs: an in-silico approach. J Biomol Struct Dyn 2022; 40:12863-12879. [PMID: 34632941 DOI: 10.1080/07391102.2021.1976278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Breast cancer type 1 susceptibility protein (BRCA1) plays an important role in maintaining genome stability and is known to interact with several proteins involved in cellular pathways, gene transcription regulation and DNA damage response. More than 40% of inherited breast cancer cases are due to BRCA1 mutation. It is also a prognostic marker in non-small cell lung cancer patients as well as a gatekeeper of cardiac function. Interaction of mutant BRCA1 with other proteins is known to disrupt the tumor suppression mechanism. Two directly interacting proteins with BRCA1 namely, DNA repair protein RAD51 (RAD51) and Aurora kinase A (AURKA), known to regulate homologous recombination (HR) and G/M cell cycle transition, respectively, form protein complex with both wild and mutant BRCA1. To analyze the interactions, protein-protein complexes were generated for each pair of proteins. In order to combat the cardiotoxic effects of cancer drugs, pharmacokinetically screened natural metabolites derived from plant, marine and bacterial sources and along with FDA-approved cancer drugs as control, were subjected to molecular docking. Piperoleine B and dihydrocircumin were the best docked natural metabolites in both RAD51 and AURKA complexes, respectively. Molecular dynamics simulation (MDS) analysis and binding free energy calculations for the best docked natural metabolite and drug for both the mutant BRCA1 complexes suggested better stability for the natural metabolites piperolein B and dihydrocurcumin as compared to drug. Thus, both natural metabolites could be further analyzed for their role against the cardiotoxic effects of cancer drugs through wet lab experiments.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sameeksha Tiwari
- Department of Biochemistry, University of Lucknow, Lucknow, India
| | - Veda P Pandey
- Department of Biochemistry, University of Lucknow, Lucknow, India
| | - Kusum Yadav
- Department of Biochemistry, University of Lucknow, Lucknow, India
| | - Upendra N Dwivedi
- Department of Biochemistry, University of Lucknow, Lucknow, India.,Institute for Development of Advanced Computing, ONGC Centre for Advanced Studies, University of Lucknow, Lucknow, India
| |
Collapse
|
15
|
Liu J, Xia KL, Wu J, Yau SST, Wei GW. Biomolecular Topology: Modelling and Analysis. ACTA MATHEMATICA SINICA, ENGLISH SERIES 2022; 38:1901-1938. [PMID: 36407804 PMCID: PMC9640850 DOI: 10.1007/s10114-022-2326-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Accepted: 07/12/2022] [Indexed: 05/25/2023]
Abstract
With the great advancement of experimental tools, a tremendous amount of biomolecular data has been generated and accumulated in various databases. The high dimensionality, structural complexity, the nonlinearity, and entanglements of biomolecular data, ranging from DNA knots, RNA secondary structures, protein folding configurations, chromosomes, DNA origami, molecular assembly, to others at the macromolecular level, pose a severe challenge in their analysis and characterization. In the past few decades, mathematical concepts, models, algorithms, and tools from algebraic topology, combinatorial topology, computational topology, and topological data analysis, have demonstrated great power and begun to play an essential role in tackling the biomolecular data challenge. In this work, we introduce biomolecular topology, which concerns the topological problems and models originated from the biomolecular systems. More specifically, the biomolecular topology encompasses topological structures, properties and relations that are emerged from biomolecular structures, dynamics, interactions, and functions. We discuss the various types of biomolecular topology from structures (of proteins, DNAs, and RNAs), protein folding, and protein assembly. A brief discussion of databanks (and databases), theoretical models, and computational algorithms, is presented. Further, we systematically review related topological models, including graphs, simplicial complexes, persistent homology, persistent Laplacians, de Rham-Hodge theory, Yau-Hausdorff distance, and the topology-based machine learning models.
Collapse
Affiliation(s)
- Jian Liu
- School of Mathematical Sciences, Hebei Normal University, Shijiazhuang, 050024 P. R. China
- Yanqi Lake Beijing Institute of Mathematical Sciences and Applications, Beijing, 101408 P. R. China
| | - Ke-Lin Xia
- School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore, 639798 Singapore
| | - Jie Wu
- Yanqi Lake Beijing Institute of Mathematical Sciences and Applications, Beijing, 101408 P. R. China
- Department of Mathematical Sciences, Tsinghua University, Beijing, 100084 P. R. China
| | - Stephen Shing-Toung Yau
- Yanqi Lake Beijing Institute of Mathematical Sciences and Applications, Beijing, 101408 P. R. China
- Department of Mathematical Sciences, Tsinghua University, Beijing, 100084 P. R. China
| | - Guo-Wei Wei
- Department of Mathematics & Department of Biochemistry and Molecular Biology & Department of Electrical and Computer Engineering, Michigan State University, Wells Hall 619 Red Cedar Road, East Lansing, MI 48824-1027 USA
| |
Collapse
|
16
|
Mishra PM, Anjum F, Uversky VN, Nandi CK. SARS-CoV-2 Spike mutations modify the interaction between virus Spike and human ACE2 receptors. Biochem Biophys Res Commun 2022; 620:8-14. [PMID: 35772213 PMCID: PMC9221686 DOI: 10.1016/j.bbrc.2022.06.064] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Accepted: 06/21/2022] [Indexed: 11/20/2022]
Abstract
The high mutability of the SARS-CoV-2 virus is a growing concern among scientific communities and health professionals since it brings the effectiveness of repurposed drugs and vaccines for COVID-19 into question. Although the mutational investigation of the Spike protein of the SARS-CoV-2 virus has been confirmed by many different researchers, there is no thorough investigation carried out at the interacting region to reveal the mutational status and its associated severity. All the energetically favorable mutations and their detailed analytical features that could impact the infection severity of the SARS-CoV-2 virus need to be identified. Therefore, we have thoroughly investigated the most important site of the SARS-CoV-2 virus, which is the interface region (Residue 417–505) of the virus Spike that interacts with the human ACE2 receptor. Further, we have utilized molecular dynamic simulation to observe the relative stability of the Spike protein with partner ACE2, as a consequence of these mutations. In our study, we have identified 52 energetically favorable Spike mutations at the interface while binding to ACE2, of which only 36 significantly enhance the stabilization of the Spike-ACE2 complex. The stability order and molecular interactions of these mutations were also identified. The highest stabilizing mutation V503D confirmed in our study is also known for neutralization resistance.
Collapse
Affiliation(s)
- Pushpendra Mani Mishra
- School of Basic Sciences, Indian Institute of Technology, Mandi, HP, 175005, India; Advanced Material Research Centre, Indian Institute of Technology, Mandi, HP, 175005, India; Bio-X Centre, Indian Institute of Technology, Mandi, HP, 175005, India
| | - Farhan Anjum
- School of Basic Sciences, Indian Institute of Technology, Mandi, HP, 175005, India; Advanced Material Research Centre, Indian Institute of Technology, Mandi, HP, 175005, India; Bio-X Centre, Indian Institute of Technology, Mandi, HP, 175005, India
| | - Vladimir N Uversky
- Department of Molecular Medicine and Byrd Alzheimer's Research Institute, Morsani College of Medicine, University of South Florida, Tampa, FL, USA.
| | - Chayan Kanti Nandi
- School of Basic Sciences, Indian Institute of Technology, Mandi, HP, 175005, India; Advanced Material Research Centre, Indian Institute of Technology, Mandi, HP, 175005, India; Bio-X Centre, Indian Institute of Technology, Mandi, HP, 175005, India.
| |
Collapse
|
17
|
Liu X, Feng H, Wu J, Xia K. Hom-Complex-Based Machine Learning (HCML) for the Prediction of Protein-Protein Binding Affinity Changes upon Mutation. J Chem Inf Model 2022; 62:3961-3969. [PMID: 36040839 DOI: 10.1021/acs.jcim.2c00580] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Protein-protein interactions (PPIs) are involved in almost all biological processes in the cell. Understanding protein-protein interactions holds the key for the understanding of biological functions, diseases and the development of therapeutics. Recently, artificial intelligence (AI) models have demonstrated great power in PPIs. However, a key issue for all AI-based PPI models is efficient molecular representations and featurization. Here, we propose Hom-complex-based PPI representation, and Hom-complex-based machine learning models for the prediction of PPI binding affinity changes upon mutation, for the first time. In our model, various Hom complexes Hom(G1, G) can be generated for the graph representation G of protein-protein complex by using different graphs G1, which reveal G1-related inner connections within the graph representation G of protein-protein complex. Further, for a specific graph G1, a series of nested Hom complexes are generated to give a multiscale characterization of the PPIs. Its persistent homology and persistent Euler characteristic are used as molecular descriptors and further combined with the machine learning model, in particular, gradient boosting tree (GBT). We systematically test our model on the two most-commonly used data sets, that is, SKEMPI and AB-Bind. It has been found that our model outperforms all the existing models as far as we know, which demonstrates the great potential of our model for the analysis of PPIs. Our model can be used for the analysis and design of efficient antibodies for SARS-CoV-2.
Collapse
Affiliation(s)
- Xiang Liu
- Chern Institute of Mathematics and LPMC, Nankai University, Tianjin, China, 300071.,Division of Mathematical Sciences, School of Physical and Mathematical Sciences Nanyang Technological University, Singapore 637371
| | - Huitao Feng
- Division of Mathematical Sciences, School of Physical and Mathematical Sciences Nanyang Technological University, Singapore 637371.,Mathematical Science Research Center, Chongqing University of Technology, Chongqing, China, 400054
| | - Jie Wu
- Yanqi Lake Beijing Institute of Mathematical Sciences and Applications (BIMSA), Beijing, China,101408
| | - Kelin Xia
- Division of Mathematical Sciences, School of Physical and Mathematical Sciences Nanyang Technological University, Singapore 637371
| |
Collapse
|
18
|
Thirumal Kumar D, Udhaya Kumar S, Jain N, Sowmya B, Balsekar K, Siva R, Kamaraj B, Sidenna M, George Priya Doss C, Zayed H. Computational structural assessment of BReast CAncer type 1 susceptibility protein (BRCA1) and BRCA1-Associated Ring Domain protein 1 (BARD1) mutations on the protein-protein interface. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2022; 130:375-397. [PMID: 35534113 DOI: 10.1016/bs.apcsb.2022.02.003] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Breast cancer type 1 susceptibility protein (BRCA1) is closely related to the BRCA2 (breast cancer type 2 susceptibility protein) and BARD1 (BRCA1-associated RING domain-1) proteins. The homodimers were formed through their RING fingers; however they form more compact heterodimers preferentially, influencing BRCA1 residues 1-109 and BARD1 residues 26-119. We implemented an integrative computational pipeline to screen all the mutations in BRCA1 and identify the most significant mutations influencing the Protein-Protein Interactions (PPI) in the BRCA1-BARD1 protein complex. The amino acids involved in the PPI regions were identified from the PDBsum database with the PDB ID: 1JM7. We screened 2118 missense mutations in BRCA1 and none in BARD1 for pathogenicity and stability and analyzed the amino acid sequences for conserved residues. We identified the most significant mutations from these screenings as V11G, M18K, L22S, and T97R positioned in the PPI regions of the BRCA1-BARD1 protein complex. We further performed protein-protein docking using the ZDOCK server. The native protein-protein complex showed the highest binding score of 2118.613, and the V11G mutant protein complex showed the least binding score of 1992.949. The other three mutation protein complexes had binding scores between the native and V11G protein complexes. Finally, a molecular dynamics simulation study using GROMACS was performed to comprehend changes in the BRCA1-BARD1 complex's binding pattern due to the mutation. From the analysis, we observed the highest deviation with lowest compactness and a decrease in the intramolecular h-bonds in the BRCA1-BARD1 protein complex with the V11G mutation compared to the native complex or the complexes with other mutations.
Collapse
Affiliation(s)
- D Thirumal Kumar
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India; Meenakshi Academy of Higher Education and Research (Deemed to be University), Chennai, Tamil Nadu, India
| | - S Udhaya Kumar
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Nikita Jain
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Baviri Sowmya
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Kamakshi Balsekar
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - R Siva
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Balu Kamaraj
- Department of Neuroscience Technology, College of Applied Medical Sciences, Imam Abdulrahman Bin Faisal University, Jubail, Saudi Arabia
| | - Mariem Sidenna
- Department of Biomedical Sciences, College of Health and Sciences, QU Health, Qatar University, Doha, Qatar
| | - C George Priya Doss
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Hatem Zayed
- Department of Biomedical Sciences, College of Health and Sciences, QU Health, Qatar University, Doha, Qatar.
| |
Collapse
|
19
|
Grassmann G, Miotto M, Di Rienzo L, Gosti G, Ruocco G, Milanetti E. A novel computational strategy for defining the minimal protein molecular surface representation. PLoS One 2022; 17:e0266004. [PMID: 35421111 PMCID: PMC9009619 DOI: 10.1371/journal.pone.0266004] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Accepted: 03/12/2022] [Indexed: 11/18/2022] Open
Abstract
Most proteins perform their biological function by interacting with one or more molecular partners. In this respect, characterizing local features of the molecular surface, that can potentially be involved in the interaction with other molecules, represents a step forward in the investigation of the mechanisms of recognition and binding between molecules. Predictive methods often rely on extensive samplings of molecular patches with the aim to identify hot spots on the surface. In this framework, analysis of large proteins and/or many molecular dynamics frames is often unfeasible due to the high computational cost. Thus, finding optimal ways to reduce the number of points to be sampled maintaining the biological information (including the surface shape) carried by the molecular surface is pivotal. In this perspective, we here present a new theoretical and computational algorithm with the aim of defining a set of molecular surfaces composed of points not uniformly distributed in space, in such a way as to maximize the information of the overall shape of the molecule by minimizing the number of total points. We test our procedure’s ability in recognizing hot-spots by describing the local shape properties of portions of molecular surfaces through a recently developed method based on the formalism of 2D Zernike polynomials. The results of this work show the ability of the proposed algorithm to preserve the key information of the molecular surface using a reduced number of points compared to the complete surface, where all points of the surface are used for the description. In fact, the methodology shows a significant gain of the information stored in the sampling procedure compared to uniform random sampling.
Collapse
Affiliation(s)
| | - Mattia Miotto
- Center for Life Nano & Neuroscience, Italian Institute of Technology, Rome, Italy
| | - Lorenzo Di Rienzo
- Center for Life Nano & Neuroscience, Italian Institute of Technology, Rome, Italy
| | - Giorgio Gosti
- Center for Life Nano & Neuroscience, Italian Institute of Technology, Rome, Italy
| | - Giancarlo Ruocco
- Center for Life Nano & Neuroscience, Italian Institute of Technology, Rome, Italy
- Department of Physics, Sapienza University, Rome, Italy
| | - Edoardo Milanetti
- Center for Life Nano & Neuroscience, Italian Institute of Technology, Rome, Italy
- Department of Physics, Sapienza University, Rome, Italy
- * E-mail:
| |
Collapse
|
20
|
Wee J, Xia K. Persistent spectral based ensemble learning (PerSpect-EL) for protein-protein binding affinity prediction. Brief Bioinform 2022; 23:6533501. [PMID: 35189639 DOI: 10.1093/bib/bbac024] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Revised: 12/30/2021] [Accepted: 01/17/2022] [Indexed: 12/14/2022] Open
Abstract
Protein-protein interactions (PPIs) play a significant role in nearly all cellular and biological activities. Data-driven machine learning models have demonstrated great power in PPIs. However, the design of efficient molecular featurization poses a great challenge for all learning models for PPIs. Here, we propose persistent spectral (PerSpect) based PPI representation and featurization, and PerSpect-based ensemble learning (PerSpect-EL) models for PPI binding affinity prediction, for the first time. In our model, a sequence of Hodge (or combinatorial) Laplacian (HL) matrices at various different scales are generated from a specially designed filtration process. PerSpect attributes, which are statistical and combinatorial properties of spectrum information from these HL matrices, are used as features for PPI characterization. Each PerSpect attribute is input into a 1D convolutional neural network (CNN), and these CNN networks are stacked together in our PerSpect-based ensemble learning models. We systematically test our model on the two most commonly used datasets, i.e. SKEMPI and AB-Bind. It has been found that our model can achieve state-of-the-art results and outperform all existing models to the best of our knowledge.
Collapse
Affiliation(s)
- JunJie Wee
- Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| | - Kelin Xia
- Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| |
Collapse
|
21
|
Ahsan T, Sajib AA. Missense variants in the TNFA epitopes and their effects on interaction with therapeutic antibodies-in silico analysis. J Genet Eng Biotechnol 2022; 20:7. [PMID: 35006391 PMCID: PMC8748575 DOI: 10.1186/s43141-021-00288-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Accepted: 12/13/2021] [Indexed: 11/10/2022]
Abstract
BACKGROUND Tumor necrosis factor alpha (TNFA) is an important cytokine that influences multiple biological processes. It is one of the key mediators of acute and chronic systemic inflammatory reactions and plays a central role in several autoimmune diseases. A number of approved monoclonal antibodies (mAbs) are widely used to subside these autoimmune diseases. However, there is a high rate of non-responsiveness to treatments with these mAbs. Therefore, it is important to be able to predict responses of the patients in an individualistic manner to these therapeutic antibodies before administration. In the present study, we used in silico tools to explore the effects of missense variants in the respective epitopes of four therapeutic anti-TNFA mAbs-adalimumab (ADA), certolizumab pegol (CZP), golimumab (GLM), and infliximab (IFX)-on their interactions with TNFA. RESULTS The binding affinities of CZP and ADA to corresponding epitopes appear to be reduced by four (TNFAR131Q, TNFAE135G, TNFAR138Q, and TNFAR138W) and two (TNFAG66C and TNFAG66S) variants, respectively. The binding of GLM and IFX appears to be affected by TNFAR141S and TNFAR138W, respectively. TNFAG66C and TNFAG66S may be associated with autoimmune diseases, whereas TNFAE135G, TNFAR138W, and TNFAR141S may be pathogenic per se. CONCLUSION These variants may contribute to the observed inter-individual variability in response to anti-TNFA mAbs treatments and be used as markers to predict responses, and thus optimize therapeutic benefits to the patients.
Collapse
Affiliation(s)
- Tamim Ahsan
- Molecular Biotechnology Division, National Institute of Biotechnology, Dhaka, 1349 Bangladesh
| | - Abu Ashfaqur Sajib
- Department of Genetic Engineering & Biotechnology, University of Dhaka, Dhaka, 1000 Bangladesh
| |
Collapse
|
22
|
Dhusia K, Madrid C, Su Z, Wu Y. EXCESP: A Structure-Based Online Database for Extracellular Interactome of Cell Surface Proteins in Humans. J Proteome Res 2022; 21:349-359. [PMID: 34978816 DOI: 10.1021/acs.jproteome.1c00612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The interactions between ectodomains of cell surface proteins are vital players in many important cellular processes, such as regulating immune responses, coordinating cell differentiation, and shaping neural plasticity. However, while the construction of a large-scale protein interactome has been greatly facilitated by the development of high-throughput experimental techniques, little progress has been made to support the discovery of extracellular interactome for cell surface proteins. Harnessed by the recent advances in computational modeling of protein-protein interactions, here we present a structure-based online database for the extracellular interactome of cell surface proteins in humans, called EXCESP. The database contains both experimentally determined and computationally predicted interactions among all type-I transmembrane proteins in humans. All structural models for these interactions and their binding affinities were further computationally modeled. Moreover, information such as expression levels of each protein in different cell types and its relation to various signaling pathways from other online resources has also been integrated into the database. In summary, the database serves as a valuable addition to the existing online resources for the study of cell surface proteins. It can contribute to the understanding of the functions of cell surface proteins in the era of systems biology.
Collapse
Affiliation(s)
- Kalyani Dhusia
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, New York 10461, United States
| | - Carlos Madrid
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, New York 10461, United States.,Laboratory for Macromolecular Analysis and Proteomics, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, New York 10461, United States
| | - Zhaoqian Su
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, New York 10461, United States
| | - Yinghao Wu
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, New York 10461, United States
| |
Collapse
|
23
|
Ndagi U, Abdullahi M, Hamza AN, Magaji MG, Mhlongo NN, Babazhitsu M, Majiya H, Makun HA, Lawal MM. Impact of Drug Repurposing on SARS-Cov-2 Main Protease. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2022; 96. [PMCID: PMC10036164 DOI: 10.1134/s0036024423030299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/25/2023]
Abstract
The recent emergence of the severe acute respiratory disease caused by a novel coronavirus remains a concern posing many challenges to public health and the global economy. The resolved crystal structure of the main protease of SARS-CoV-2 or SCV2 (Mpro) has led to its identification as an attractive target for designing potent antiviral drugs. Herein, we provide a comparative molecular impact of hydroxychloroquine (HCQ), remdesivir, and β-D-N4-Hydroxycytidine (NHC) binding on SCV2 Mpro using various computational approaches like molecular docking and molecular dynamics (MD) simulation. Data analyses showed that HCQ, remdesivir, and NHC binding to SARS-CoV-2 Mpro decrease the protease loop capacity to fluctuate. These binding influences the drugs’ optimum orientation in the conformational space of SCV2 Mpro and produce noticeable steric effects on the interactive residues. An increased hydrogen bond formation was observed in SCV2 Mpro–NHC complex with a decreased receptor residence time during NHC binding. The binding mode of remdesivir to SCV2 Mpro differs from other drugs having van der Waals interaction as the force stabilizing protein–remdesivir complex. Electrostatic interaction dominates in the SCV2 Mpro−HCQ and SCV2 Mpro–NHC. Residue Glu166 was highly involved in the stability of remdesivir and NHC binding at the SCV2 Mpro active site, while Asp187 provides stability for HCQ binding.
Collapse
Affiliation(s)
- Umar Ndagi
- Africa Centre of Excellence for Mycotoxin and Food Safety, Federal University of Technology, Minna, Nigeria
| | - Maryam Abdullahi
- Faculty of Pharmaceutical Sciences, Ahmadu Bello University, Zaria, Kaduna State, Nigeria
| | - Asmau N. Hamza
- Faculty of Pharmaceutical Sciences, Ahmadu Bello University, Zaria, Kaduna State, Nigeria
| | - Mohd G. Magaji
- Faculty of Pharmaceutical Sciences, Ahmadu Bello University, Zaria, Kaduna State, Nigeria
| | - Ndumiso N. Mhlongo
- Department of Medical Biochemistry, School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, 4001 Durban, South Africa
| | - Makun Babazhitsu
- Department of Medical Microbiology and Parasitology, Faculty of Basic Clinical Sciences, College of Health Sciences, Usman Danfodio University, Sokoto, Nigeria
| | - Hussaini Majiya
- Department of Microbiology, Ibrahim Badamasi Babangida University, Lapai, Niger State, Nigeria
| | - Hussaini Anthony Makun
- Africa Centre of Excellence for Mycotoxin and Food Safety, Federal University of Technology, Minna, Nigeria
| | - Monsurat M. Lawal
- Department of Medical Biochemistry, School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, 4001 Durban, South Africa
| |
Collapse
|
24
|
Qiu L, Song J, Zhang JZH. Computational Alanine Scanning Reveals Common Features of TCR/pMHC Recognition in HLA-DQ8-Associated Celiac Disease. Methods Mol Biol 2022; 2385:293-312. [PMID: 34888725 DOI: 10.1007/978-1-0716-1767-0_13] [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] [Indexed: 06/13/2023]
Abstract
In HLA-DQ8-associated celiac disease, Gliadin-γ1 or Gliadin-α1 peptide is presented to the cell surface and recognized by several types of T-cell receptor (TCR), but it is still unclear how the TCR, peptide, and the major histocompatibility complex (MHC) act together to trigger celiac disease. For now, most of the analysis is based on static crystal structures. And the detailed information about these structures based on energetic interaction is still lacking. Here, we took four types of celiac disease-related MHC-peptide-TCR structures from three patients to perform computational alanine scanning calculations using the molecular mechanics generalized born surface area (MM/GBSA) approach combined with a recently developed interaction entropy (IE) method to identify the key residues on TCR, peptide, and MHC. Our study aims to shed some light on the interaction mechanism of this complex protein interaction system. Based on detailed computational analysis and mutational calculations, important binding interactions in these triple-interaction complexes are analyzed, and critical residues responsible for TCR/pMHC recognition pattern in HLA-DQ8-associated celiac disease are presented. These detailed analysis and computational result should help shed light on our understanding of the celiac disease and the development of the medical treatment.
Collapse
Affiliation(s)
- Linqiong Qiu
- Shanghai Key Laboratory of Green Chemistry and Chemical Processes, School of Chemistry and Molecular Engineering, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai, China
| | - Jianing Song
- NYU-ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai, China
| | - John Z H Zhang
- Shanghai Key Laboratory of Green Chemistry and Chemical Processes, School of Chemistry and Molecular Engineering, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai, China.
- NYU-ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai, China.
- Department of Chemistry, New York University, New York, NY, USA.
| |
Collapse
|
25
|
Cho JH, Ju WS, Seo SY, Kim BH, Kim JS, Kim JG, Park SJ, Choo YK. The Potential Role of Human NME1 in Neuronal Differentiation of Porcine Mesenchymal Stem Cells: Application of NB-hNME1 as a Human NME1 Suppressor. Int J Mol Sci 2021; 22:ijms222212194. [PMID: 34830075 PMCID: PMC8619003 DOI: 10.3390/ijms222212194] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Revised: 11/02/2021] [Accepted: 11/08/2021] [Indexed: 12/31/2022] Open
Abstract
This study aimed to investigate the effects of the human macrophage (MP) secretome in cellular xenograft rejection. The role of human nucleoside diphosphate kinase A (hNME1), from the secretome of MPs involved in the neuronal differentiation of miniature pig adipose tissue-derived mesenchymal stem cells (mp AD-MSCs), was evaluated by proteomic analysis. Herein, we first demonstrate that hNME1 strongly binds to porcine ST8 alpha-N-acetyl-neuraminide alpha-2,8-sialyltransferase 1 (pST8SIA1), which is a ganglioside GD3 synthase. When hNME1 binds with pST8SIA1, it induces degradation of pST8SIA1 in mp AD-MSCs, thereby inhibiting the expression of ganglioside GD3 followed by decreased neuronal differentiation of mp AD-MSCs. Therefore, we produced nanobodies (NBs) named NB-hNME1 that bind to hNME1 specifically, and the inhibitory effect of NB-hNME1 was evaluated for blocking the binding between hNME1 and pST8SIA1. Consequently, NB-hNME1 effectively blocked the binding of hNME1 to pST8SIA1, thereby recovering the expression of ganglioside GD3 and neuronal differentiation of mp AD-MSCs. Our findings suggest that mp AD-MSCs could be a potential candidate for use as an additive, such as an immunosuppressant, in stem cell transplantation.
Collapse
Affiliation(s)
- Jin Hyoung Cho
- Department of Biological Science, College of Natural Sciences, Wonkwang University, 460, Iksan-daero, Iksan-si 54538, Korea; (J.H.C.); (W.S.J.); (S.Y.S.); (J.-G.K.); (S.J.P.)
- GreenBio Corp. Central Research, 201-19, Bubaljungand-ro, Bubal-eup, Icheon-si 17321, Korea
| | - Won Seok Ju
- Department of Biological Science, College of Natural Sciences, Wonkwang University, 460, Iksan-daero, Iksan-si 54538, Korea; (J.H.C.); (W.S.J.); (S.Y.S.); (J.-G.K.); (S.J.P.)
- Institute for Glycoscience, Wonkwang University, 460, Iksan-daero, Iksan-si 54538, Korea
| | - Sang Young Seo
- Department of Biological Science, College of Natural Sciences, Wonkwang University, 460, Iksan-daero, Iksan-si 54538, Korea; (J.H.C.); (W.S.J.); (S.Y.S.); (J.-G.K.); (S.J.P.)
| | - Bo Hyun Kim
- CHA Fertility Center Bundang, 59, Yatap-ro, Bundang-gu, Seongnam-si 13496, Korea;
| | - Ji-Su Kim
- Primate Resources Center (PRC), Korea Research Institute of Bioscience and Biotechnology, 181, Ipsin-gil, Jeongeup-si 56216, Korea;
| | - Jong-Geol Kim
- Department of Biological Science, College of Natural Sciences, Wonkwang University, 460, Iksan-daero, Iksan-si 54538, Korea; (J.H.C.); (W.S.J.); (S.Y.S.); (J.-G.K.); (S.J.P.)
| | - Soon Ju Park
- Department of Biological Science, College of Natural Sciences, Wonkwang University, 460, Iksan-daero, Iksan-si 54538, Korea; (J.H.C.); (W.S.J.); (S.Y.S.); (J.-G.K.); (S.J.P.)
| | - Young-Kug Choo
- Department of Biological Science, College of Natural Sciences, Wonkwang University, 460, Iksan-daero, Iksan-si 54538, Korea; (J.H.C.); (W.S.J.); (S.Y.S.); (J.-G.K.); (S.J.P.)
- Institute for Glycoscience, Wonkwang University, 460, Iksan-daero, Iksan-si 54538, Korea
- Correspondence: ; Tel.: +82-63-850-6087; Fax: +82-63-857-8837
| |
Collapse
|
26
|
Flores SC, Alexiou A, Glaros A. Mining the Protein Data Bank to improve prediction of changes in protein-protein binding. PLoS One 2021; 16:e0257614. [PMID: 34727109 PMCID: PMC8562805 DOI: 10.1371/journal.pone.0257614] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2021] [Accepted: 09/05/2021] [Indexed: 12/23/2022] Open
Abstract
Predicting the effect of mutations on protein-protein interactions is important for relating structure to function, as well as for in silico affinity maturation. The effect of mutations on protein-protein binding energy (ΔΔG) can be predicted by a variety of atomic simulation methods involving full or limited flexibility, and explicit or implicit solvent. Methods which consider only limited flexibility are naturally more economical, and many of them are quite accurate, however results are dependent on the atomic coordinate set used. In this work we perform a sequence and structure based search of the Protein Data Bank to find additional coordinate sets and repeat the calculation on each. The method increases precision and Positive Predictive Value, and decreases Root Mean Square Error, compared to using single structures. Given the ongoing growth of near-redundant structures in the Protein Data Bank, our method will only increase in applicability and accuracy.
Collapse
Affiliation(s)
| | - Athanasios Alexiou
- Department of Computer Science and Biomedical Informatics, University of Thessaly, Volos, Greece
| | - Anastasios Glaros
- Eukaryotic Single Cell Genomics Facility, Science For Life Laboratory, Stockholm, Sweden
| |
Collapse
|
27
|
Abbasi WA, Abbas SA, Andleeb S. PANDA: Predicting the change in proteins binding affinity upon mutations by finding a signal in primary structures. J Bioinform Comput Biol 2021; 19:2150015. [PMID: 34126874 DOI: 10.1142/s0219720021500153] [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] [Indexed: 11/18/2022]
Abstract
Accurately determining a change in protein binding affinity upon mutations is important to find novel therapeutics and to assist mutagenesis studies. Determination of change in binding affinity upon mutations requires sophisticated, expensive, and time-consuming wet-lab experiments that can be supported with computational methods. Most of the available computational prediction techniques depend upon protein structures that bound their applicability to only protein complexes with recognized 3D structures. In this work, we explore the sequence-based prediction of change in protein binding affinity upon mutation and question the effectiveness of [Formula: see text]-fold cross-validation (CV) across mutations adopted in previous studies to assess the generalization ability of such predictors with no known mutation during training. We have used protein sequence information instead of protein structures along with machine learning techniques to accurately predict the change in protein binding affinity upon mutation. Our proposed sequence-based novel change in protein binding affinity predictor called PANDA performs comparably to the existing methods gauged through an appropriate CV scheme and an external independent test dataset. On an external test dataset, our proposed method gives a maximum Pearson correlation coefficient of 0.52 in comparison to the state-of-the-art existing protein structure-based method called MutaBind which gives a maximum Pearson correlation coefficient of 0.59. Our proposed protein sequence-based method, to predict a change in binding affinity upon mutations, has wide applicability and comparable performance in comparison to existing protein structure-based methods. We made PANDA easily accessible through a cloud-based webserver and python code available at https://sites.google.com/view/wajidarshad/software and https://github.com/wajidarshad/panda, respectively.
Collapse
Affiliation(s)
- Wajid Arshad Abbasi
- Computational Biology and Data Analysis Lab., Department of Computer Sciences & Information Technology, King Abdullah Campus, University of Azad Jammu & Kashmir, Muzaffarabad, AJ&K 13100, Pakistan
| | - Syed Ali Abbas
- Computational Biology and Data Analysis Lab., Department of Computer Sciences & Information Technology, King Abdullah Campus, University of Azad Jammu & Kashmir, Muzaffarabad, AJ&K 13100, Pakistan
| | - Saiqa Andleeb
- Biotechnology Lab., Department of Zoology, King Abdullah Campus, University of Azad Jammu & Kashmir, Muzaffarabad, AJ&K 13100, Pakistan
| |
Collapse
|
28
|
Sequeiros-Borja CE, Surpeta B, Brezovsky J. Recent advances in user-friendly computational tools to engineer protein function. Brief Bioinform 2021; 22:bbaa150. [PMID: 32743637 PMCID: PMC8138880 DOI: 10.1093/bib/bbaa150] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Revised: 06/03/2020] [Accepted: 06/16/2020] [Indexed: 12/14/2022] Open
Abstract
Progress in technology and algorithms throughout the past decade has transformed the field of protein design and engineering. Computational approaches have become well-engrained in the processes of tailoring proteins for various biotechnological applications. Many tools and methods are developed and upgraded each year to satisfy the increasing demands and challenges of protein engineering. To help protein engineers and bioinformaticians navigate this emerging wave of dedicated software, we have critically evaluated recent additions to the toolbox regarding their application for semi-rational and rational protein engineering. These newly developed tools identify and prioritize hotspots and analyze the effects of mutations for a variety of properties, comprising ligand binding, protein-protein and protein-nucleic acid interactions, and electrostatic potential. We also discuss notable progress to target elusive protein dynamics and associated properties like ligand-transport processes and allosteric communication. Finally, we discuss several challenges these tools face and provide our perspectives on the further development of readily applicable methods to guide protein engineering efforts.
Collapse
Affiliation(s)
- Carlos Eduardo Sequeiros-Borja
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw, Warsaw, Poland
| | - Bartłomiej Surpeta
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw, Warsaw, Poland
| | - Jan Brezovsky
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw
| |
Collapse
|
29
|
Arines FM, Hamlin AJ, Yang X, Liu YYJ, Li M. A selective transmembrane recognition mechanism by a membrane-anchored ubiquitin ligase adaptor. J Cell Biol 2021; 220:211632. [PMID: 33351099 PMCID: PMC7759299 DOI: 10.1083/jcb.202001116] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2020] [Revised: 09/10/2020] [Accepted: 10/30/2020] [Indexed: 12/30/2022] Open
Abstract
While it is well-known that E3 ubiquitin ligases can selectively ubiquitinate membrane proteins in response to specific environmental cues, the underlying mechanisms for the selectivity are poorly understood. In particular, the role of transmembrane regions, if any, in target recognition remains an open question. Here, we describe how Ssh4, a yeast E3 ligase adaptor, recognizes the PQ-loop lysine transporter Ypq1 only after lysine starvation. We show evidence of an interaction between two transmembrane helices of Ypq1 (TM5 and TM7) and the single transmembrane helix of Ssh4. This interaction is regulated by the conserved PQ motif. Strikingly, recent structural studies of the PQ-loop family have suggested that TM5 and TM7 undergo major conformational changes during substrate transport, implying that transport-associated conformational changes may determine the selectivity. These findings thus provide critical information concerning the regulatory mechanism through which transmembrane domains can be specifically recognized in response to changing environmental conditions.
Collapse
Affiliation(s)
- Felichi Mae Arines
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI
| | - Aaron Jeremy Hamlin
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI
| | - Xi Yang
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI
| | - Yun-Yu Jennifer Liu
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI
| | - Ming Li
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI
| |
Collapse
|
30
|
Kandeel M, Yamamoto M, Tani H, Kobayashi A, Gohda J, Kawaguchi Y, Park BK, Kwon HJ, Inoue JI, Alkattan A. Discovery of New Fusion Inhibitor Peptides against SARS-CoV-2 by Targeting the Spike S2 Subunit. Biomol Ther (Seoul) 2021; 29:282-289. [PMID: 33424013 PMCID: PMC8094075 DOI: 10.4062/biomolther.2020.201] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Revised: 12/08/2020] [Accepted: 12/10/2020] [Indexed: 12/31/2022] Open
Abstract
A novel coronavirus, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), caused a worldwide pandemic. Our aim in this study is to produce new fusion inhibitors against SARS-CoV-2, which can be the basis for developing new antiviral drugs. The fusion core comprising the heptad repeat domains (HR1 and HR2) of SARS-CoV-2 spike (S) were used to design the peptides. A total of twelve peptides were generated, comprising a short or truncated 24-mer (peptide #1), a long 36-mer peptide (peptide #2), and ten peptide #2 analogs. In contrast to SARS-CoV, SARS-CoV-2 S-mediated cell-cell fusion cannot be inhibited with a minimal length, 24-mer peptide. Peptide #2 demonstrated potent inhibition of SARS-CoV-2 S-mediated cell-cell fusion at 1 µM concentration. Three peptide #2 analogs showed IC50 values in the low micromolar range (4.7-9.8 µM). Peptide #2 inhibited the SARS-CoV-2 pseudovirus assay at IC50=1.49 µM. Given their potent inhibition of viral activity and safety and lack of cytotoxicity, these peptides provide an attractive avenue for the development of new prophylactic and therapeutic agents against SARS-CoV-2.
Collapse
Affiliation(s)
- Mahmoud Kandeel
- Department of Biomedical Sciences, College of Veterinary Medicine, King Faisal University, Al-Ahsa 31982, Saudi Arabia.,Department of Pharmacology, Faculty of Veterinary Medicine, Kafrelsheikh University, Kafrelsheikh 33516, Egypt
| | - Mizuki Yamamoto
- Research Center for Asian Infectious Diseases, Institute of Medical Science, The University of Tokyo, Tokyo 108-8639, Japan
| | - Hideki Tani
- Department of Virology, Toyama Institute of Health, Toyama 939-0363, Japan
| | - Ayako Kobayashi
- Research Center for Asian Infectious Diseases, Institute of Medical Science, The University of Tokyo, Tokyo 108-8639, Japan
| | - Jin Gohda
- Research Center for Asian Infectious Diseases, Institute of Medical Science, The University of Tokyo, Tokyo 108-8639, Japan
| | - Yasushi Kawaguchi
- Research Center for Asian Infectious Diseases, Institute of Medical Science, The University of Tokyo, Tokyo 108-8639, Japan.,Division of Molecular Virology, Department of Microbiology and Immunology, The Institute of Medical Science, The University of Tokyo, Tokyo 108-8639, Japan
| | - Byoung Kwon Park
- Department of Microbiology, Hallym University College of Medicine, Chuncheon 24252, Republic of Korea
| | - Hyung-Joo Kwon
- Department of Microbiology, Hallym University College of Medicine, Chuncheon 24252, Republic of Korea
| | - Jun-Ichiro Inoue
- Senior Professor Office, The University of Tokyo, Tokyo 108-8639, Japan
| | - Abdallah Alkattan
- Department of Biomedical Sciences, College of Veterinary Medicine, King Faisal University, Al-Ahsa 31982, Saudi Arabia
| |
Collapse
|
31
|
Wang B, Su Z, Wu Y. Computational Assessment of Protein-Protein Binding Affinity by Reverse Engineering the Energetics in Protein Complexes. GENOMICS PROTEOMICS & BIOINFORMATICS 2021; 19:1012-1022. [PMID: 33838354 PMCID: PMC9403033 DOI: 10.1016/j.gpb.2021.03.004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/10/2018] [Revised: 03/07/2019] [Accepted: 05/17/2019] [Indexed: 11/29/2022]
Abstract
The cellular functions of proteins are maintained by forming diverse complexes. The stability of these complexes is quantified by the measurement of binding affinity, and mutations that alter the binding affinity can cause various diseases such as cancer and diabetes. As a result, accurate estimation of the binding stability and the effects of mutations on changes of binding affinity is a crucial step to understanding the biological functions of proteins and their dysfunctional consequences. It has been hypothesized that the stability of a protein complex is dependent not only on the residues at its binding interface by pairwise interactions but also on all other remaining residues that do not appear at the binding interface. Here, we computationally reconstruct the binding affinity by decomposing it into the contributions of interfacial residues and other non-interfacial residues in a protein complex. We further assume that the contributions of both interfacial and non-interfacial residues to the binding affinity depend on their local structural environments such as solvent-accessible surfaces and secondary structural types. The weights of all corresponding parameters are optimized by Monte-Carlo simulations. After cross-validation against a large-scale dataset, we show that the model not only shows a strong correlation between the absolute values of the experimental and calculated binding affinities, but can also be an effective approach to predict the relative changes of binding affinity from mutations. Moreover, we have found that the optimized weights of many parameters can capture the first-principle chemical and physical features of molecular recognition, therefore reversely engineering the energetics of protein complexes. These results suggest that our method can serve as a useful addition to current computational approaches for predicting binding affinity and understanding the molecular mechanism of protein–protein interactions.
Collapse
Affiliation(s)
- Bo Wang
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, NY 10461, USA
| | - Zhaoqian Su
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, NY 10461, USA
| | - Yinghao Wu
- Department of Systems and Computational Biology, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, NY 10461, USA.
| |
Collapse
|
32
|
Ajadi MB, Soremekun OS, Adewumi AT, Kumalo HM, Soliman MES. Functional Analysis of Single Nucleotide Polymorphism in ZUFSP Protein and Implication in Pathogenesis. Protein J 2021; 40:28-40. [PMID: 33512633 DOI: 10.1007/s10930-021-09962-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/04/2021] [Indexed: 11/25/2022]
Abstract
Researches have revealed that functional non-synonymous Single Nucleotide Polymorphism (nsSNPs) present in the Zinc-finger with UFM1-Specific Peptidase domain protein (ZUFSP) may be involved in genetic instability and carcinogenesis. For the first time, we employed in-silico approach using predictive tools to identify and validate potential nsSNPs that could be pathogenic. Our result revealed that 8 nsSNPs (rs 112738382, rs 140094037, rs 201652589, rs 201847265, rs 202076827, rs 373634906, rs 375114528, rs 772591104) are pathogenic after being subjected to rigorous filtering process. The structural impact of the nsSNPs on ZUFSP structure indicated that the nsSNPs affect the stability of the protein by lowering ZUFSP protein stability. Furthermore, conservation analysis showed that rs 201652589, rs 140094037, rs 201847265, and rs 772591104 were highly conserved. Interestingly, the protein-protein affinity between ZUFSP and Ubiquitin was altered rs 201652589, rs 140094037, rs 201847265, and rs 772591104 had a binding affinity of - 0.46, - 0.83, - 1.62, and - 1.12 kcal/mol respectively. Our study has been able to identify potential nsSNPs that could be used as genetic biomarkers for some diseases arising as a result of aberration in the ZUFSP structure, however, being a predictive study, the identified nsSNPs need to be experimentally investigated.
Collapse
Affiliation(s)
- Mary B Ajadi
- Department of Medical Biochemistry, School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Howard Campus, Durban, 4000, South Africa
- Chemical Pathology Department, Faculty of Basic Medical Sciences, College of Health Sciences, Ladoke Akintola University of Technology, PMB 4400, Osogbo, Nigeria
| | - Opeyemi S Soremekun
- Molecular Bio-Computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban, 4001, South Africa
| | - Adeniyi T Adewumi
- Molecular Bio-Computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban, 4001, South Africa
| | - Hezekiel M Kumalo
- Department of Medical Biochemistry, School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Howard Campus, Durban, 4000, South Africa
| | - Mahmoud E S Soliman
- Molecular Bio-Computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban, 4001, South Africa.
| |
Collapse
|
33
|
Gonzalez TR, Martin KP, Barnes JE, Patel JS, Ytreberg FM. Assessment of software methods for estimating protein-protein relative binding affinities. PLoS One 2020; 15:e0240573. [PMID: 33347442 PMCID: PMC7751979 DOI: 10.1371/journal.pone.0240573] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2020] [Accepted: 12/07/2020] [Indexed: 11/19/2022] Open
Abstract
A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83-98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.
Collapse
Affiliation(s)
- Tawny R. Gonzalez
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
| | - Kyle P. Martin
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| | - Jonathan E. Barnes
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| | - Jagdish Suresh Patel
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Biological Sciences, University of Idaho, Moscow, Idaho, United States of America
| | - F. Marty Ytreberg
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| |
Collapse
|
34
|
Su Z, Wang B, Almo SC, Wu Y. Understanding the Targeting Mechanisms of Multi-Specific Biologics in Immunotherapy with Multiscale Modeling. iScience 2020; 23:101835. [PMID: 33305190 PMCID: PMC7710644 DOI: 10.1016/j.isci.2020.101835] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Revised: 09/29/2020] [Accepted: 11/17/2020] [Indexed: 11/30/2022] Open
Abstract
Immunotherapeutics are frequently associated with adverse side effects due to the elicitation of global immune modulation. To lower the risk of these side effects, recombinant DNA technology is employed to enhance the selectivity of cell targeting by genetically fusing different biomolecules, yielding new species referred to as multi-specific biologics. The design of new multi-specific biologics is a central challenge for the realization of new immunotherapies. To understand the molecular determinants responsible for regulating the binding between multi-specific biologics and surface-bound membrane receptors, we developed a multiscale computational framework that integrates various simulation approaches covering different timescales and spatial resolutions. Our model system of multi-specific biologics contains two natural ligands of immune receptors, which are covalently tethered by a peptide linker. Using this method, a number of interesting features of multi-specific biologics were identified. Our study therefore provides an important strategy to design the next-generation biologics for immunotherapy. Two proteins are connected by different linkers as a model of bispecific biologics Conformational dynamics of biologics are captured by microsecond MD simulations Coarse-grained simulations are used to test binding between biologics and receptors Biologics with long and flexible linkers are more efficient in targeting receptors
Collapse
Affiliation(s)
- Zhaoqian Su
- Department of Systems and Computational Biology, Albert Einstein College of Medicine of Yeshiva University, 1300 Morris Park Avenue, Bronx, NY 10461, USA
| | - Bo Wang
- Department of Systems and Computational Biology, Albert Einstein College of Medicine of Yeshiva University, 1300 Morris Park Avenue, Bronx, NY 10461, USA
| | - Steven C Almo
- Department of Biochemistry, Albert Einstein College of Medicine of Yeshiva University, 1300 Morris Park Avenue, Bronx, NY 10461, USA.,Department of Physiology and Biophysics, Albert Einstein College of Medicine of Yeshiva University, 1300 Morris Park Avenue, Bronx, NY 10461, USA
| | - Yinghao Wu
- Department of Systems and Computational Biology, Albert Einstein College of Medicine of Yeshiva University, 1300 Morris Park Avenue, Bronx, NY 10461, USA
| |
Collapse
|
35
|
de Groot NS, Torrent Burgas M. Bacteria use structural imperfect mimicry to hijack the host interactome. PLoS Comput Biol 2020; 16:e1008395. [PMID: 33275611 PMCID: PMC7744059 DOI: 10.1371/journal.pcbi.1008395] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Revised: 12/16/2020] [Accepted: 09/23/2020] [Indexed: 12/25/2022] Open
Abstract
Bacteria use protein-protein interactions to infect their hosts and hijack fundamental pathways, which ensures their survival and proliferation. Hence, the infectious capacity of the pathogen is closely related to its ability to interact with host proteins. Here, we show that hubs in the host-pathogen interactome are isolated in the pathogen network by adapting the geometry of the interacting interfaces. An imperfect mimicry of the eukaryotic interfaces allows pathogen proteins to actively bind to the host's target while preventing deleterious effects on the pathogen interactome. Understanding how bacteria recognize eukaryotic proteins may pave the way for the rational design of new antibiotic molecules.
Collapse
Affiliation(s)
- Natalia Sanchez de Groot
- Gene Function and Evolution Lab, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, Barcelona, Spain
- * E-mail: (NSdG); (MTB)
| | - Marc Torrent Burgas
- Systems Biology of Infection Lab, Department of Biochemistry and Molecular Biology, Biosciences Faculty, Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Spain
- * E-mail: (NSdG); (MTB)
| |
Collapse
|
36
|
Milanetti E, Miotto M, Di Rienzo L, Monti M, Gosti G, Ruocco G. 2D Zernike polynomial expansion: Finding the protein-protein binding regions. Comput Struct Biotechnol J 2020; 19:29-36. [PMID: 33363707 PMCID: PMC7750141 DOI: 10.1016/j.csbj.2020.11.051] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Revised: 11/26/2020] [Accepted: 11/28/2020] [Indexed: 01/26/2023] Open
Abstract
We present a method for efficiently and effectively assessing whether and where two proteins can interact with each other to form a complex. This is still largely an open problem, even for those relatively few cases where the 3D structure of both proteins is known. In fact, even if much of the information about the interaction is encoded in the chemical and geometric features of the structures, the set of possible contact patches and of their relative orientations are too large to be computationally affordable in a reasonable time, thus preventing the compilation of reliable interactome. Our method is able to rapidly and quantitatively measure the geometrical shape complementarity between interacting proteins, comparing their molecular iso-electron density surfaces expanding the surface patches in term of 2D Zernike polynomials. We first test the method against the real binding region of a large dataset of known protein complexes, reaching a success rate of 0.72. We then apply the method for the blind recognition of binding sites, identifying the real region of interaction in about 60% of the analyzed cases. Finally, we investigate how the efficiency in finding the right binding region depends on the surface roughness as a function of the expansion order.
Collapse
Affiliation(s)
- Edoardo Milanetti
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy.,Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| | - Mattia Miotto
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy.,Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| | - Lorenzo Di Rienzo
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| | - Michele Monti
- Centre for Genomic Regulation (CRG), The Barcelona Institute for Science and Technology, Dr. Aiguader 88, 08003 Barcelona, Spain.,RNA System Biology Lab, Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Giorgio Gosti
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| | - Giancarlo Ruocco
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy.,Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| |
Collapse
|
37
|
Vihinen M. Functional effects of protein variants. Biochimie 2020; 180:104-120. [PMID: 33164889 DOI: 10.1016/j.biochi.2020.10.009] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Revised: 10/15/2020] [Accepted: 10/19/2020] [Indexed: 12/11/2022]
Abstract
Genetic and other variations frequently affect protein functions. Scientific articles can contain confusing descriptions about which function or property is affected, and in many cases the statements are pure speculation without any experimental evidence. To clarify functional effects of protein variations of genetic or non-genetic origin, a systematic conceptualisation and framework are introduced. This framework describes protein functional effects on abundance, activity, specificity and affinity, along with countermeasures, which allow cells, tissues and organisms to tolerate, avoid, repair, attenuate or resist (TARAR) the effects. Effects on abundance discussed include gene dosage, restricted expression, mis-localisation and degradation. Enzymopathies, effects on kinetics, allostery and regulation of protein activity are subtopics for the effects of variants on activity. Variation outcomes on specificity and affinity comprise promiscuity, specificity, affinity and moonlighting. TARAR mechanisms redress variations with active and passive processes including chaperones, redundancy, robustness, canalisation and metabolic and signalling rewiring. A framework for pragmatic protein function analysis and presentation is introduced. All of the mechanisms and effects are described along with representative examples, most often in relation to diseases. In addition, protein function is discussed from evolutionary point of view. Application of the presented framework facilitates unambiguous, detailed and specific description of functional effects and their systematic study.
Collapse
Affiliation(s)
- Mauno Vihinen
- Department of Experimental Medical Science, BMC B13, Lund University, SE-22 184, Lund, Sweden.
| |
Collapse
|
38
|
Huang X, Zheng W, Pearce R, Zhang Y. SSIPe: accurately estimating protein-protein binding affinity change upon mutations using evolutionary profiles in combination with an optimized physical energy function. Bioinformatics 2020; 36:2429-2437. [PMID: 31830252 DOI: 10.1093/bioinformatics/btz926] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2019] [Revised: 11/08/2019] [Accepted: 12/09/2019] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION Most proteins perform their biological functions through interactions with other proteins in cells. Amino acid mutations, especially those occurring at protein interfaces, can change the stability of protein-protein interactions (PPIs) and impact their functions, which may cause various human diseases. Quantitative estimation of the binding affinity changes (ΔΔGbind) caused by mutations can provide critical information for protein function annotation and genetic disease diagnoses. RESULTS We present SSIPe, which combines protein interface profiles, collected from structural and sequence homology searches, with a physics-based energy function for accurate ΔΔGbind estimation. To offset the statistical limits of the PPI structure and sequence databases, amino acid-specific pseudocounts were introduced to enhance the profile accuracy. SSIPe was evaluated on large-scale experimental data containing 2204 mutations from 177 proteins, where training and test datasets were stringently separated with the sequence identity between proteins from the two datasets below 30%. The Pearson correlation coefficient between estimated and experimental ΔΔGbind was 0.61 with a root-mean-square-error of 1.93 kcal/mol, which was significantly better than the other methods. Detailed data analyses revealed that the major advantage of SSIPe over other traditional approaches lies in the novel combination of the physical energy function with the new knowledge-based interface profile. SSIPe also considerably outperformed a former profile-based method (BindProfX) due to the newly introduced sequence profiles and optimized pseudocount technique that allows for consideration of amino acid-specific prior mutation probabilities. AVAILABILITY AND IMPLEMENTATION Web-server/standalone program, source code and datasets are freely available at https://zhanglab.ccmb.med.umich.edu/SSIPe and https://github.com/tommyhuangthu/SSIPe. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | - Wei Zheng
- Department of Computational Medicine and Bioinformatics
| | - Robin Pearce
- Department of Computational Medicine and Bioinformatics
| | - Yang Zhang
- Department of Computational Medicine and Bioinformatics.,Department of Biological Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
| |
Collapse
|
39
|
Di Rienzo L, Milanetti E, Testi C, Montemiglio LC, Baiocco P, Boffi A, Ruocco G. A novel strategy for molecular interfaces optimization: The case of Ferritin-Transferrin receptor interaction. Comput Struct Biotechnol J 2020; 18:2678-2686. [PMID: 33101606 PMCID: PMC7548301 DOI: 10.1016/j.csbj.2020.09.020] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Revised: 09/10/2020] [Accepted: 09/11/2020] [Indexed: 11/24/2022] Open
Abstract
Protein-protein interactions regulate almost all cellular functions and rely on a fine tune of surface amino acids properties involved on both molecular partners. The disruption of a molecular association can be caused even by a single residue mutation, often leading to a pathological modification of a biochemical pathway. Therefore the evaluation of the effects of amino acid substitutions on binding, and the ad hoc design of protein-protein interfaces, is one of the biggest challenges in computational biology. Here, we present a novel strategy for computational mutation and optimization of protein-protein interfaces. Modeling the interaction surface properties using the Zernike polynomials, we describe the shape and electrostatics of binding sites with an ordered set of descriptors, making possible the evaluation of complementarity between interacting surfaces. With a Monte Carlo approach, we obtain protein mutants with controlled molecular complementarities. Applying this strategy to the relevant case of the interaction between Ferritin and Transferrin Receptor, we obtain a set of Ferritin mutants with increased or decreased complementarity. The extensive molecular dynamics validation of the method results confirms its efficacy, showing that this strategy represents a very promising approach in designing correct molecular interfaces.
Collapse
Affiliation(s)
- Lorenzo Di Rienzo
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| | - Edoardo Milanetti
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Claudia Testi
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| | | | - Paola Baiocco
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
- Department of Biochemical Sciences ‘A. Rossi Fanelli’ Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Alberto Boffi
- Department of Biochemical Sciences ‘A. Rossi Fanelli’ Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Giancarlo Ruocco
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| |
Collapse
|
40
|
Huang X, Pearce R, Zhang Y. EvoEF2: accurate and fast energy function for computational protein design. Bioinformatics 2020; 36:1135-1142. [PMID: 31588495 DOI: 10.1093/bioinformatics/btz740] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2019] [Revised: 09/19/2019] [Accepted: 09/25/2019] [Indexed: 01/26/2023] Open
Abstract
MOTIVATION The accuracy and success rate of de novo protein design remain limited, mainly due to the parameter over-fitting of current energy functions and their inability to discriminate incorrect designs from correct designs. RESULTS We developed an extended energy function, EvoEF2, for efficient de novo protein sequence design, based on a previously proposed physical energy function, EvoEF. Remarkably, EvoEF2 recovered 32.5%, 47.9% and 22.3% of all, core and surface residues for 148 test monomers, and was generally applicable to protein-protein interaction design, as it recapitulated 30.9%, 42.4%, 31.3% and 21.4% of all, core, interface and surface residues for 88 test dimers, significantly outperforming EvoEF on the native sequence recapitulation. We further used I-TASSER to evaluate the foldability of the 148 designed monomer sequences, where all of them were predicted to fold into structures with high fold- and atomic-level similarity to their corresponding native structures, as demonstrated by the fact that 87.8% of the predicted structures shared a root-mean-square-deviation less than 2 Å to their native counterparts. The study also demonstrated that the usefulness of physical energy functions is highly correlated with the parameter optimization processes, and EvoEF2, with parameters optimized using sequence recapitulation, is more suitable for computational protein sequence design than EvoEF, which was optimized on thermodynamic mutation data. AVAILABILITY AND IMPLEMENTATION The source code of EvoEF2 and the benchmark datasets are freely available at https://zhanglab.ccmb.med.umich.edu/EvoEF. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Xiaoqiang Huang
- Department of Computational Medicine and Bioinformatics, MI 48109, USA
| | - Robin Pearce
- Department of Computational Medicine and Bioinformatics, MI 48109, USA
| | - Yang Zhang
- Department of Computational Medicine and Bioinformatics, MI 48109, USA.,Department of Biological Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
| |
Collapse
|
41
|
Zhang N, Lu H, Chen Y, Zhu Z, Yang Q, Wang S, Li M. PremPRI: Predicting the Effects of Missense Mutations on Protein-RNA Interactions. Int J Mol Sci 2020; 21:ijms21155560. [PMID: 32756481 PMCID: PMC7432928 DOI: 10.3390/ijms21155560] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Revised: 07/28/2020] [Accepted: 07/30/2020] [Indexed: 12/23/2022] Open
Abstract
Protein–RNA interactions are crucial for many cellular processes, such as protein synthesis and regulation of gene expression. Missense mutations that alter protein–RNA interaction may contribute to the pathogenesis of many diseases. Here, we introduce a new computational method PremPRI, which predicts the effects of single mutations occurring in RNA binding proteins on the protein–RNA interactions by calculating the binding affinity changes quantitatively. The multiple linear regression scoring function of PremPRI is composed of three sequence- and eight structure-based features, and is parameterized on 248 mutations from 50 protein–RNA complexes. Our model shows a good agreement between calculated and experimental values of binding affinity changes with a Pearson correlation coefficient of 0.72 and the corresponding root-mean-square error of 0.76 kcal·mol−1, outperforming three other available methods. PremPRI can be used for finding functionally important variants, understanding the molecular mechanisms, and designing new protein–RNA interaction inhibitors.
Collapse
|
42
|
Beg MA, Athar F. Molecular modeling and in silico characterization of mycobacterial Rv3101c and Rv3102c proteins: prerequisite molecular target in cell division. ACTA ACUST UNITED AC 2020. [DOI: 10.15406/ppij.2020.08.00300] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
43
|
Shringari SR, Giannakoulias S, Ferrie JJ, Petersson EJ. Rosetta custom score functions accurately predict ΔΔG of mutations at protein-protein interfaces using machine learning. Chem Commun (Camb) 2020; 56:6774-6777. [PMID: 32441721 DOI: 10.1039/d0cc01959c] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Protein-protein interfaces play essential roles in a variety of biological processes and many therapeutic molecules are targeted at these interfaces. However, accurate predictions of the effects of interfacial mutations to identify "hotspots" have remained elusive despite the myriad of modeling and machine learning methods tested. Here, for the first time, we demonstrate that nonlinear reweighting of energy terms from Rosetta, through the use of machine learning, exhibits improved predictability of ΔΔG values associated with interfacial mutations.
Collapse
Affiliation(s)
- Sumant R Shringari
- Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, PA 19104, USA.
| | | | | | | |
Collapse
|
44
|
Pahari S, Li G, Murthy AK, Liang S, Fragoza R, Yu H, Alexov E. SAAMBE-3D: Predicting Effect of Mutations on Protein-Protein Interactions. Int J Mol Sci 2020; 21:E2563. [PMID: 32272725 PMCID: PMC7177817 DOI: 10.3390/ijms21072563] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Revised: 04/04/2020] [Accepted: 04/05/2020] [Indexed: 12/26/2022] Open
Abstract
Maintaining wild type protein-protein interactions is essential for the normal function of cell and any mutation that alter their characteristics can cause disease. Therefore, the ability to correctly and quickly predict the effect of amino acid mutations is crucial for understanding disease effects and to be able to carry out genome-wide studies. Here, we report a new development of the SAAMBE method, SAAMBE-3D, which is a machine learning-based approach, resulting in accurate predictions and is extremely fast. It achieves the Pearson correlation coefficient ranging from 0.78 to 0.82 depending on the training protocol in benchmarking five-fold validation test against the SKEMPI v2.0 database and outperforms currently existing algorithms on various blind-tests. Furthermore, optimized and tested via five-fold cross-validation on the Cornell University dataset, the SAAMBE-3D achieves AUC of 1.0 and 0.96 on a homo and hereto-dimer test datasets. Another important feature of SAAMBE-3D is that it is very fast, it takes less than a fraction of a second to complete a prediction. SAAMBE-3D is available as a web server and as well as a stand-alone code, the last one being another important feature allowing other researchers to directly download the code and run it on their local computer. Combined all together, SAAMBE-3D is an accurate and fast software applicable for genome-wide studies to assess the effect of amino acid mutations on protein-protein interactions. The webserver and the stand-alone codes (SAAMBE-3D for predicting the change of binding free energy and SAAMBE-3D-DN for predicting if the mutation is disruptive or non-disruptive) are available.
Collapse
Affiliation(s)
- Swagata Pahari
- Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA; (S.P.); (G.L.); (A.K.M.)
| | - Gen Li
- Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA; (S.P.); (G.L.); (A.K.M.)
| | - Adithya Krishna Murthy
- Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA; (S.P.); (G.L.); (A.K.M.)
| | - Siqi Liang
- Department of Computational Biology, Cornell University, Ithaca, NY 14850, USA; (S.L.); (R.F.); (H.Y.)
| | - Robert Fragoza
- Department of Computational Biology, Cornell University, Ithaca, NY 14850, USA; (S.L.); (R.F.); (H.Y.)
| | - Haiyuan Yu
- Department of Computational Biology, Cornell University, Ithaca, NY 14850, USA; (S.L.); (R.F.); (H.Y.)
| | - Emil Alexov
- Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA; (S.P.); (G.L.); (A.K.M.)
| |
Collapse
|
45
|
Zhou G, Chen M, Ju CJT, Wang Z, Jiang JY, Wang W. Mutation effect estimation on protein-protein interactions using deep contextualized representation learning. NAR Genom Bioinform 2020; 2:lqaa015. [PMID: 32166223 PMCID: PMC7059401 DOI: 10.1093/nargab/lqaa015] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2019] [Revised: 01/20/2020] [Accepted: 02/24/2020] [Indexed: 12/14/2022] Open
Abstract
The functional impact of protein mutations is reflected on the alteration of conformation and thermodynamics of protein–protein interactions (PPIs). Quantifying the changes of two interacting proteins upon mutations is commonly carried out by computational approaches. Hence, extensive research efforts have been put to the extraction of energetic or structural features on proteins, followed by statistical learning methods to estimate the effects of mutations on PPI properties. Nonetheless, such features require extensive human labors and expert knowledge to obtain, and have limited abilities to reflect point mutations. We present an end-to-end deep learning framework, MuPIPR (Mutation Effects in Protein–protein Interaction PRediction Using Contextualized Representations), to estimate the effects of mutations on PPIs. MuPIPR incorporates a contextualized representation mechanism of amino acids to propagate the effects of a point mutation to surrounding amino acid representations, therefore amplifying the subtle change in a long protein sequence. On top of that, MuPIPR leverages a Siamese residual recurrent convolutional neural encoder to encode a wild-type protein pair and its mutation pair. Multi-layer perceptron regressors are applied to the protein pair representations to predict the quantifiable changes of PPI properties upon mutations. Experimental evaluations show that, with only sequence information, MuPIPR outperforms various state-of-the-art systems on estimating the changes of binding affinity for SKEMPI v1, and offers comparable performance on SKEMPI v2. Meanwhile, MuPIPR also demonstrates state-of-the-art performance on estimating the changes of buried surface areas. The software implementation is available at https://github.com/guangyu-zhou/MuPIPR.
Collapse
Affiliation(s)
- Guangyu Zhou
- Department of Computer Science, University of California, Los Angeles, CA 90095, USA
| | - Muhao Chen
- Department of Computer Science, University of California, Los Angeles, CA 90095, USA.,Department of Computer and Information Science, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Chelsea J T Ju
- Department of Computer Science, University of California, Los Angeles, CA 90095, USA
| | - Zheng Wang
- Department of Computer Science, University of California, Los Angeles, CA 90095, USA
| | - Jyun-Yu Jiang
- Department of Computer Science, University of California, Los Angeles, CA 90095, USA
| | - Wei Wang
- Department of Computer Science, University of California, Los Angeles, CA 90095, USA
| |
Collapse
|
46
|
Zhang N, Chen Y, Lu H, Zhao F, Alvarez RV, Goncearenco A, Panchenko AR, Li M. MutaBind2: Predicting the Impacts of Single and Multiple Mutations on Protein-Protein Interactions. iScience 2020; 23:100939. [PMID: 32169820 PMCID: PMC7068639 DOI: 10.1016/j.isci.2020.100939] [Citation(s) in RCA: 80] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2019] [Revised: 11/21/2019] [Accepted: 02/20/2020] [Indexed: 01/17/2023] Open
Abstract
Missense mutations may affect proteostasis by destabilizing or over-stabilizing protein complexes and changing the pathway flux. Predicting the effects of stabilizing mutations on protein-protein interactions is notoriously difficult because existing experimental sets are skewed toward mutations reducing protein-protein binding affinity and many computational methods fail to correctly evaluate their effects. To address this issue, we developed a method MutaBind2, which estimates the impacts of single as well as multiple mutations on protein-protein interactions. MutaBind2 employs only seven features, and the most important of them describe interactions of proteins with the solvent, evolutionary conservation of the site, and thermodynamic stability of the complex and each monomer. This approach shows a distinct improvement especially in evaluating the effects of mutations increasing binding affinity. MutaBind2 can be used for finding disease driver mutations, designing stable protein complexes, and discovering new protein-protein interaction inhibitors.
Collapse
Affiliation(s)
- Ning Zhang
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Yuting Chen
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Haoyu Lu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Feiyang Zhao
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Roberto Vera Alvarez
- National Center for Biotechnology Information, National Institutes of Health, Bethesda, MD 20894, USA
| | - Alexander Goncearenco
- National Center for Biotechnology Information, National Institutes of Health, Bethesda, MD 20894, USA
| | - Anna R Panchenko
- National Center for Biotechnology Information, National Institutes of Health, Bethesda, MD 20894, USA.
| | - Minghui Li
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China.
| |
Collapse
|
47
|
Ndagi U, Abdullahi M, Hamza AN, Soliman ME. An analogue of a kinase inhibitor exhibits subjective characteristics that contribute to its inhibitory activities as a potential anti-cancer candidate: insights through computational biomolecular modelling of UM-164 binding with lyn protein. RSC Adv 2020; 10:145-161. [PMID: 35492550 PMCID: PMC9047091 DOI: 10.1039/c9ra07204g] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2019] [Accepted: 12/12/2019] [Indexed: 12/31/2022] Open
Abstract
The recent emergence of lyn kinase as a driver of aggressive behaviour in triple-negative breast cancer (TNBC) remains a major concern posing a burden for people living with breast cancer and drug development. The binding of UM-164 to lyn protein has been noted to impact the conformational dynamics required for drug fitness. Herein, we provide the first account of the molecular impact of an experimental drug, UM-164 binding on lyn protein using various computational approaches including molecular docking and molecular dynamics simulation. These computational modelling methods enabled us to analyse parameters, for example principal component analysis (PCA), dynamics cross-correlation matrices (DCCM) analysis, hydrogen bond occupancy, thermodynamics calculation and ligand–residue interaction. Findings from these analyses revealed that UM-164 exhibited a higher binding affinity of −9.9 kcal mol−1 with lyn protein than Dasatinib, with a binding affinity of −8.3 kcal mol−1 on docking. It was observed that the binding of UM-164 to lyn protein decreases the capacity of its loop to fluctuate, influences the ligand optimum orientation on the conformational space of lyn protein, and increases the hydrogen bond formation in the lyn-UM-164 system. Also, an increase in drug binding energy of UM-164 was recorded with increasing residue correlation in the lyn-UM-164 system. It is quite informative to note that Met85 was a key stabilising factor in the binding of UM-164 to lyn protein. These findings can provide important insights that will potentially serve as a baseline in the design of novel lyn inhibitors. It could also stimulate further research into multidimensional approaches required to curb the influence of lyn protein in TNBC. This study provides the first account of the molecular impact of UM-164 binding on lyn protein using various computational approaches.![]()
Collapse
Affiliation(s)
- Umar Ndagi
- Faculty of Natural Sciences
- Ibrahim Badamasi Babangida University
- Nigeria
| | - Maryam Abdullahi
- Molecular Bio-Computation and Drug Design Research Group
- School of Health Sciences
- University of KwaZulu-Natal
- Durban 4000
- South Africa
| | - Asmau N. Hamza
- Faculty of Pharmaceutical Sciences
- Ahmadu Bello University
- Zaria
- Nigeria
| | - Mahmoud E. Soliman
- Molecular Bio-Computation and Drug Design Research Group
- School of Health Sciences
- University of KwaZulu-Natal
- Durban 4000
- South Africa
| |
Collapse
|
48
|
Robust Prediction of Single and Multiple Point Protein Mutations Stability Changes. Biomolecules 2019; 10:biom10010067. [PMID: 31906171 PMCID: PMC7023245 DOI: 10.3390/biom10010067] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2019] [Revised: 12/19/2019] [Accepted: 12/20/2019] [Indexed: 11/16/2022] Open
Abstract
Accurate prediction of protein stability changes resulting from amino acid substitutions is of utmost importance in medicine to better understand which mutations are deleterious, leading to diseases, and which are neutral. Since conducting wet lab experiments to get a better understanding of protein mutations is costly and time consuming, and because of huge number of possible mutations the need of computational methods that could accurately predict effects of amino acid mutations is of greatest importance. In this research, we present a robust methodology to predict the energy changes of a proteins upon mutations. The proposed prediction scheme is based on two step algorithm that is a Holdout Random Sampler followed by a neural network model for regression. The Holdout Random Sampler is utilized to analysis the energy change, the corresponding uncertainty, and to obtain a set of admissible energy changes, expressed as a cumulative distribution function. These values are further utilized to train a simple neural network model that can predict the energy changes. Results were blindly tested (validated) against experimental energy changes, giving Pearson correlation coefficients of 0.66 for Single Point Mutations and 0.77 for Multiple Point Mutations. These results confirm the successfulness of our method, since it outperforms majority of previous studies in this field.
Collapse
|
49
|
Kulandaisamy A, Zaucha J, Sakthivel R, Frishman D, Michael Gromiha M. Pred‐MutHTP: Prediction of disease‐causing and neutral mutations in human transmembrane proteins. Hum Mutat 2019; 41:581-590. [DOI: 10.1002/humu.23961] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2019] [Revised: 11/05/2019] [Accepted: 11/20/2019] [Indexed: 12/24/2022]
Affiliation(s)
- A. Kulandaisamy
- Department of Biotechnology, Bhupat and Jyoti Mehta School of BioSciencesIndian Institute of Technology MadrasChennai Tamilnadu India
| | - Jan Zaucha
- Department of Bioinformatics, Wissenschaftszentrum WeihenstephanTechnische Universität MünchenFreising Germany
| | - Ramasamy Sakthivel
- Department of Biotechnology, Bhupat and Jyoti Mehta School of BioSciencesIndian Institute of Technology MadrasChennai Tamilnadu India
| | - Dmitrij Frishman
- Department of Bioinformatics, Wissenschaftszentrum WeihenstephanTechnische Universität MünchenFreising Germany
- Department of BioinformaticsPeter the Great St. Petersburg Polytechnic UniversitySt. Petersburg Russian Federation
| | - M. Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of BioSciencesIndian Institute of Technology MadrasChennai Tamilnadu India
- Advanced Computational Drug Discovery Unit, Tokyo Tech World Research Hub Initiative (WRHI), Institute of Innovative ResearchTokyo Institute of TechnologyYokohama Japan
| |
Collapse
|
50
|
Jayaraman M, Rajendra SK, Ramadas K. Structural insight into conformational dynamics of non-active site mutations in KasA: A Mycobacterium tuberculosis target protein. Gene 2019; 720:144082. [DOI: 10.1016/j.gene.2019.144082] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2018] [Revised: 08/23/2019] [Accepted: 08/26/2019] [Indexed: 10/26/2022]
|