1
|
Kozome D, Sljoka A, Laurino P. Remote loop evolution reveals a complex biological function for chitinase enzymes beyond the active site. Nat Commun 2024; 15:3227. [PMID: 38622119 PMCID: PMC11018821 DOI: 10.1038/s41467-024-47588-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2024] [Accepted: 04/08/2024] [Indexed: 04/17/2024] Open
Abstract
Loops are small secondary structural elements that play a crucial role in the emergence of new enzyme functions. However, the evolutionary molecular mechanisms how proteins acquire these loop elements and obtain new function is poorly understood. To address this question, we study glycoside hydrolase family 19 (GH19) chitinase-an essential enzyme family for pathogen degradation in plants. By revealing the evolutionary history and loops appearance of GH19 chitinase, we discover that one loop which is remote from the catalytic site, is necessary to acquire the new antifungal activity. We demonstrate that this remote loop directly accesses the fungal cell wall, and surprisingly, it needs to adopt a defined structure supported by long-range intramolecular interactions to perform its function. Our findings prove that nature applies this strategy at the molecular level to achieve a complex biological function while maintaining the original activity in the catalytic pocket, suggesting an alternative way to design new enzyme function.
Collapse
Affiliation(s)
- Dan Kozome
- Protein Engineering and Evolution Unit, Okinawa Institute of Science and Technology Graduate University (OIST), Okinawa, 904-0495, Japan
| | - Adnan Sljoka
- Center for Advanced Intelligence Project, RIKEN, Tokyo, 103-0027, Japan
- Department of Chemistry, York University, Toronto, ON, M3J 1P3, Canada
| | - Paola Laurino
- Protein Engineering and Evolution Unit, Okinawa Institute of Science and Technology Graduate University (OIST), Okinawa, 904-0495, Japan.
- Institute for Protein Research, Osaka University, Suita, Japan.
| |
Collapse
|
2
|
Hsu MF, Sriramoju MK, Lai CH, Chen YR, Huang JS, Ko TP, Huang KF, Hsu STD. Structure, dynamics, and stability of the smallest and most complex 7 1 protein knot. J Biol Chem 2024; 300:105553. [PMID: 38072060 PMCID: PMC10840475 DOI: 10.1016/j.jbc.2023.105553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Revised: 11/21/2023] [Accepted: 12/04/2023] [Indexed: 12/29/2023] Open
Abstract
Proteins can spontaneously tie a variety of intricate topological knots through twisting and threading of the polypeptide chains. Recently developed artificial intelligence algorithms have predicted several new classes of topological knotted proteins, but the predictions remain to be authenticated experimentally. Here, we showed by X-ray crystallography and solution-state NMR spectroscopy that Q9PR55, an 89-residue protein from Ureaplasma urealyticum, possesses a novel 71 knotted topology that is accurately predicted by AlphaFold 2, except for the flexible N terminus. Q9PR55 is monomeric in solution, making it the smallest and most complex knotted protein known to date. In addition to its exceptional chemical stability against urea-induced unfolding, Q9PR55 is remarkably robust to resist the mechanical unfolding-coupled proteolysis by a bacterial proteasome, ClpXP. Our results suggest that the mechanical resistance against pulling-induced unfolding is determined by the complexity of the knotted topology rather than the size of the molecule.
Collapse
Affiliation(s)
- Min-Feng Hsu
- Institute of Biological Chemistry, Academia Sinica, Taipei, Taiwan
| | | | - Chih-Hsuan Lai
- Institute of Biological Chemistry, Academia Sinica, Taipei, Taiwan
| | - Yun-Ru Chen
- Academia Sinica Protein Clinic, Academia Sinica, Taipei, Taiwan
| | - Jing-Siou Huang
- Academia Sinica Protein Clinic, Academia Sinica, Taipei, Taiwan
| | - Tzu-Ping Ko
- Institute of Biological Chemistry, Academia Sinica, Taipei, Taiwan
| | - Kai-Fa Huang
- Institute of Biological Chemistry, Academia Sinica, Taipei, Taiwan; Academia Sinica Protein Clinic, Academia Sinica, Taipei, Taiwan
| | - Shang-Te Danny Hsu
- Institute of Biological Chemistry, Academia Sinica, Taipei, Taiwan; Academia Sinica Protein Clinic, Academia Sinica, Taipei, Taiwan; Institute of Biochemical Sciences, National Taiwan University, Taipei, Taiwan; International Institute for Sustainability with Knotted Chiral Meta Matter (WPI-SKCM(2)), Hiroshima University, Higashihiroshima, Japan.
| |
Collapse
|
3
|
Chen J, Kuhn LA, Raschka S. Techniques for Developing Reliable Machine Learning Classifiers Applied to Understanding and Predicting Protein:Protein Interaction Hot Spots. Methods Mol Biol 2024; 2714:235-268. [PMID: 37676603 DOI: 10.1007/978-1-0716-3441-7_14] [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: 09/08/2023]
Abstract
With machine learning now transforming the sciences, successful prediction of biological structure or activity is mainly limited by the extent and quality of data available for training, the astute choice of features for prediction, and thorough assessment of the robustness of prediction on a variety of new cases. In this chapter, we address these issues while developing and sharing protocols to build a robust dataset and rigorously compare several predictive classifiers using the open-source Python machine learning library, scikit-learn. We show how to evaluate whether enough data has been used for training and whether the classifier has been overfit to training data. The most telling experiment is 500-fold repartitioning of the training and test sets, followed by prediction, which gives a good indication of whether a classifier performs consistently well on different datasets. An intuitive method is used to quantify which features are most important for correct prediction.The resulting well-trained classifier, hotspotter, can robustly predict the small subset of amino acid residues on the surface of a protein that are energetically most important for binding a protein partner: the interaction hot spots. Hotspotter has been trained and tested here on a curated dataset assembled from 1046 non-redundant alanine scanning mutation sites with experimentally measured change in binding free energy values from 97 different protein complexes; this dataset is available to download. The accessible surface area of the wild-type residue at a given site and its degree of evolutionary conservation proved the most important features to identify hot spots. A variant classifier was trained and validated for proteins where only the amino acid sequence is available, augmented by secondary structure assignment. This version of hotspotter requiring fewer features is almost as robust as the structure-based classifier. Application to the ACE2 (angiotensin converting enzyme 2) receptor, which mediates COVID-19 virus entry into human cells, identified the critical hot spot triad of ACE2 residues at the center of the small interface with the CoV-2 spike protein. Hotspotter results can be used to guide the strategic design of protein interfaces and ligands and also to identify likely interfacial residues for protein:protein docking.
Collapse
Affiliation(s)
- Jiaxing Chen
- Bioinformatics and Genomics Graduate Program, Pennsylvania State University, University Park, PA, USA
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA
| | - Leslie A Kuhn
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA.
| | - Sebastian Raschka
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA
- Department of Statistics, University of Wisconsin-Madison, Madison, WI, USA
| |
Collapse
|
4
|
Pražnikar J. Using graphlet degree vectors to predict atomic displacement parameters in protein structures. Acta Crystallogr D Struct Biol 2023; 79:1109-1119. [PMID: 37987168 PMCID: PMC10833351 DOI: 10.1107/s2059798323009142] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Accepted: 10/17/2023] [Indexed: 11/22/2023] Open
Abstract
In structural biology, atomic displacement parameters, commonly used in the form of B values, describe uncertainties in atomic positions. Their distribution over the structure can provide hints on local structural reliability and mobility. A spatial macromolecular model can be represented by a graph whose nodes are atoms and whose edges correspond to all interatomic contacts within a certain distance. Small connected subgraphs, called graphlets, provide information about the wiring of a particular atom. The multiple linear regression approach based on this information aims to predict a distribution of values of isotropic atomic displacement parameters (B values) within a protein structure, given the atomic coordinates and molecular packing. By modeling the dynamic component of atomic uncertainties, this method allows the B values obtained from experimental crystallographic or cryo-electron microscopy studies to be reproduced relatively well.
Collapse
Affiliation(s)
- Jure Pražnikar
- Faculty of Mathematics, Natural Sciences and Information Technologies, University of Primorska, Glagoljaška 8, Koper, Slovenia
- Department of Biochemistry, Molecular and Structural Biology, Institute Jožef Stefan, Jamova 39, Ljubljana, Slovenia
| |
Collapse
|
5
|
Tyler S, Laforge C, Guzzo A, Nicolaï A, Maisuradze GG, Senet P. Einstein Model of a Graph to Characterize Protein Folded/Unfolded States. Molecules 2023; 28:6659. [PMID: 37764437 PMCID: PMC10536427 DOI: 10.3390/molecules28186659] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Revised: 09/11/2023] [Accepted: 09/14/2023] [Indexed: 09/29/2023] Open
Abstract
The folded structures of proteins can be accurately predicted by deep learning algorithms from their amino-acid sequences. By contrast, in spite of decades of research studies, the prediction of folding pathways and the unfolded and misfolded states of proteins, which are intimately related to diseases, remains challenging. A two-state (folded/unfolded) description of protein folding dynamics hides the complexity of the unfolded and misfolded microstates. Here, we focus on the development of simplified order parameters to decipher the complexity of disordered protein structures. First, we show that any connected, undirected, and simple graph can be associated with a linear chain of atoms in thermal equilibrium. This analogy provides an interpretation of the usual topological descriptors of a graph, namely the Kirchhoff index and Randić resistance, in terms of effective force constants of a linear chain. We derive an exact relation between the Kirchhoff index and the average shortest path length for a linear graph and define the free energies of a graph using an Einstein model. Second, we represent the three-dimensional protein structures by connected, undirected, and simple graphs. As a proof of concept, we compute the topological descriptors and the graph free energies for an all-atom molecular dynamics trajectory of folding/unfolding events of the proteins Trp-cage and HP-36 and for the ensemble of experimental NMR models of Trp-cage. The present work shows that the local, nonlocal, and global force constants and free energies of a graph are promising tools to quantify unfolded/disordered protein states and folding/unfolding dynamics. In particular, they allow the detection of transient misfolded rigid states.
Collapse
Affiliation(s)
- Steve Tyler
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université de Bourgogne, 21078 Dijon CEDEX, France
| | - Christophe Laforge
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université de Bourgogne, 21078 Dijon CEDEX, France
| | - Adrien Guzzo
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université de Bourgogne, 21078 Dijon CEDEX, France
| | - Adrien Nicolaï
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université de Bourgogne, 21078 Dijon CEDEX, France
| | - Gia G. Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853, USA
| | - Patrick Senet
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université de Bourgogne, 21078 Dijon CEDEX, France
| |
Collapse
|
6
|
Abstract
A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.
Collapse
Affiliation(s)
- Ke-Wei Chen
- Lab of Computional Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Tian-Yu Sun
- Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Yun-Dong Wu
- Lab of Computional Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
- Shenzhen Bay Laboratory, Shenzhen 518132, China
| |
Collapse
|
7
|
Fowler NJ, Albalwi MF, Lee S, Hounslow AM, Williamson MP. Improved methodology for protein NMR structure calculation using hydrogen bond restraints and ANSURR validation: The SH2 domain of SH2B1. Structure 2023; 31:975-986.e3. [PMID: 37311460 DOI: 10.1016/j.str.2023.05.012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Revised: 05/02/2023] [Accepted: 05/18/2023] [Indexed: 06/15/2023]
Abstract
Protein structures calculated using NMR data are less accurate and less well-defined than they could be. Here we use the program ANSURR to show that this deficiency is at least in part due to a lack of hydrogen bond restraints. We describe a protocol to introduce hydrogen bond restraints into the structure calculation of the SH2 domain from SH2B1 in a systematic and transparent way and show that the structures generated are more accurate and better defined as a result. We also show that ANSURR can be used as a guide to know when the structure calculation is good enough to stop.
Collapse
Affiliation(s)
- Nicholas J Fowler
- School of Biosciences, University of Sheffield, S10 2TN Sheffield, UK.
| | - Marym F Albalwi
- School of Biosciences, University of Sheffield, S10 2TN Sheffield, UK
| | - Subin Lee
- School of Biosciences, University of Sheffield, S10 2TN Sheffield, UK
| | - Andrea M Hounslow
- School of Biosciences, University of Sheffield, S10 2TN Sheffield, UK
| | - Mike P Williamson
- School of Biosciences, University of Sheffield, S10 2TN Sheffield, UK.
| |
Collapse
|
8
|
Wang G. Thermal Ring-Based Heat Switches in Hyperthermophilic Class II Bacterial Fructose Aldolase. ACS OMEGA 2023; 8:24624-24634. [PMID: 37457467 PMCID: PMC10339327 DOI: 10.1021/acsomega.3c03001] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/01/2023] [Accepted: 05/24/2023] [Indexed: 07/18/2023]
Abstract
Both thermophilic and hyperthermophilic enzymes in bacterial and archaeal species are activated above a specific temperature threshold but inactivated at another higher temperature. However, the underlying structural basis for these two heat switches is still unresolved. Here, graph theory was used to test if the temperature-dependent noncovalent interactions and metal bridges as identified in a series of crystal structures of the class II bacterial fructose 1,6-bisphosphate aldolase homodimer or homotetramer with or without natural substrates and products bound could form systematic fluidic grid-like mesh networks with topological grids as thermal rings to regulate their structural thermostability and functional thermoactivity. The results indicated that the second biggest grid in the Thermus aquaticus fructose 1,6-diphosphate aldolase dimer may control the specific temperature threshold to release the swapping flexible active sites at the dimeric interface for heat-evoked activation. Meanwhile, the third biggest grid may serve as a necessary structural motif against heat inactivation. Finally, the smallest grid may act as a stiff thermostable anchor. Its dissociation at the maximal melting temperature threshold may stop the catalytic activity. Taken as a whole, this computational study may render the structural motifs for the optimal growth temperature and the extreme heat stability of hyperthermophiles.
Collapse
|
9
|
Kumar P, Vyas P, Faisal SM, Chang YF, Akif M. Crystal structure of a variable region segment of Leptospira host-interacting outer surface protein, LigA, reveals the orientation of Ig-like domains. Int J Biol Macromol 2023:125445. [PMID: 37336372 DOI: 10.1016/j.ijbiomac.2023.125445] [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: 02/07/2023] [Revised: 06/13/2023] [Accepted: 06/14/2023] [Indexed: 06/21/2023]
Abstract
Leptospiral immunoglobulin-like (Lig) protein family is a surface-exposed protein from the pathogenic Leptospira. The Lig protein family has been identified as an essential virulence factor of L. interrogan. One of the family members, LigA, contains 13 homologous tandem repeats of bacterial Ig-like (Big) domains in its extracellular portion. It is crucial in binding with the host's Extracellular matrices (ECM) and complement factors. However, its vital role in the invasion and evasion of pathogenic Leptospira, structural details, and domain organization of the extracellular portion of this protein are not explored thoroughly. Here, we described the first high-resolution crystal structure of a variable region segment (LigA8-9) of LigA at 1.87 Å resolution. The structure showed some remarkably distinctive aspects compared with the most closely related Immunoglobulin superfamily (IgSF) members. The structure illustrated the relative orientation of two domains and highlighted the role of the linker region in the domain orientation. We also observed an apparent electron density of Ca2+ ions coordinated with a proper interacting geometry within the protein. Molecular dynamic simulations demonstrated the involvement of a linker salt bridge in providing rigidity between the two domains. Our study proposes an overall arrangement of Ig-like domains in the LigA protein. The structural understanding of the extracellular portion of LigA and its interaction with the ECM provides insight into developing new therapeutics directed toward leptospirosis.
Collapse
Affiliation(s)
- Pankaj Kumar
- Laboratory of Structural Biology, Department of Biochemistry, School of Life Sciences, University of Hyderabad, Gachibowli, Hyderabad, India
| | - Pallavi Vyas
- Laboratory of Vaccine Immunology, National Institute of Animal Biotechnology, Gachibowli, Hyderabad, Telangana, India
| | - Syed M Faisal
- Laboratory of Vaccine Immunology, National Institute of Animal Biotechnology, Gachibowli, Hyderabad, Telangana, India
| | - Yung-Fu Chang
- Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY, USA
| | - Mohd Akif
- Laboratory of Structural Biology, Department of Biochemistry, School of Life Sciences, University of Hyderabad, Gachibowli, Hyderabad, India.
| |
Collapse
|
10
|
Jones RD, Jones AM. Model of ligand-triggered information transmission in G-protein coupled receptor complexes. Front Endocrinol (Lausanne) 2023; 14:1111594. [PMID: 37361529 PMCID: PMC10286511 DOI: 10.3389/fendo.2023.1111594] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 03/21/2023] [Indexed: 06/28/2023] Open
Abstract
We present a model for the effects of ligands on information transmission in G-Protein Coupled Receptor (GPCR) complexes. The model is built ab initio entirely on principles of statistical mechanics and tenets of information transmission theory and was validated in part using agonist-induced effector activity and signaling bias for the angiotensin- and adrenergic-mediated signaling pathways, with in vitro observations of phosphorylation sites on the C tail of the GPCR complex, and single-cell information-transmission experiments. The model extends traditional kinetic models that form the basis for many existing models of GPCR signaling. It is based on maximizing the rates of entropy production and information transmission through the GPCR complex. The model predicts that (1) phosphatase-catalyzed reactions, as opposed to kinase-catalyzed reactions, on the C-tail and internal loops of the GPCR are responsible for controlling the signaling activity, (2) signaling favors the statistical balance of the number of switches in the ON state and the number in the OFF state, and (3) biased-signaling response depends discontinuously on ligand concentration.
Collapse
Affiliation(s)
- Roger D. Jones
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
- European Centre for Living Technology, Ca’ Foscari University of Venice, Venice, Italy
- Systems Engineering and Research Center, Stevens Institute of Technology, Hoboken, NJ, United States
| | - Alan M. Jones
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
- Department of Pharmacology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
| |
Collapse
|
11
|
Huang SK, Picard LP, Rahmatullah RSM, Pandey A, Van Eps N, Sunahara RK, Ernst OP, Sljoka A, Prosser RS. Mapping the conformational landscape of the stimulatory heterotrimeric G protein. Nat Struct Mol Biol 2023; 30:502-511. [PMID: 36997760 DOI: 10.1038/s41594-023-00957-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Accepted: 02/24/2023] [Indexed: 04/01/2023]
Abstract
Heterotrimeric G proteins serve as membrane-associated signaling hubs, in concert with their cognate G-protein-coupled receptors. Fluorine nuclear magnetic resonance spectroscopy was employed to monitor the conformational equilibria of the human stimulatory G-protein α subunit (Gsα) alone, in the intact Gsαβ1γ2 heterotrimer or in complex with membrane-embedded human adenosine A2A receptor (A2AR). The results reveal a concerted equilibrium that is strongly affected by nucleotide and interactions with the βγ subunit, the lipid bilayer and A2AR. The α1 helix of Gsα exhibits significant intermediate timescale dynamics. The α4β6 loop and α5 helix undergo membrane/receptor interactions and order-disorder transitions respectively, associated with G-protein activation. The αN helix adopts a key functional state that serves as an allosteric conduit between the βγ subunit and receptor, while a significant fraction of the ensemble remains tethered to the membrane and receptor upon activation.
Collapse
Affiliation(s)
- Shuya Kate Huang
- Department of Chemistry, University of Toronto, UTM, Mississauga, Ontario, Canada
| | | | - Rima S M Rahmatullah
- Department of Chemistry, University of Toronto, UTM, Mississauga, Ontario, Canada
| | - Aditya Pandey
- Department of Chemistry, University of Toronto, UTM, Mississauga, Ontario, Canada
| | - Ned Van Eps
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
| | - Roger K Sunahara
- Department of Pharmacology, University of California San Diego School of Medicine, La Jolla, CA, USA
| | - Oliver P Ernst
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada
| | - Adnan Sljoka
- RIKEN Center for Advanced Intelligence Project, RIKEN, Tokyo, Japan.
| | - R Scott Prosser
- Department of Chemistry, University of Toronto, UTM, Mississauga, Ontario, Canada.
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada.
| |
Collapse
|
12
|
Yang Z, Zeng X, Zhao Y, Chen R. AlphaFold2 and its applications in the fields of biology and medicine. Signal Transduct Target Ther 2023; 8:115. [PMID: 36918529 PMCID: PMC10011802 DOI: 10.1038/s41392-023-01381-z] [Citation(s) in RCA: 51] [Impact Index Per Article: 51.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2022] [Revised: 12/27/2022] [Accepted: 02/16/2023] [Indexed: 03/16/2023] Open
Abstract
AlphaFold2 (AF2) is an artificial intelligence (AI) system developed by DeepMind that can predict three-dimensional (3D) structures of proteins from amino acid sequences with atomic-level accuracy. Protein structure prediction is one of the most challenging problems in computational biology and chemistry, and has puzzled scientists for 50 years. The advent of AF2 presents an unprecedented progress in protein structure prediction and has attracted much attention. Subsequent release of structures of more than 200 million proteins predicted by AF2 further aroused great enthusiasm in the science community, especially in the fields of biology and medicine. AF2 is thought to have a significant impact on structural biology and research areas that need protein structure information, such as drug discovery, protein design, prediction of protein function, et al. Though the time is not long since AF2 was developed, there are already quite a few application studies of AF2 in the fields of biology and medicine, with many of them having preliminarily proved the potential of AF2. To better understand AF2 and promote its applications, we will in this article summarize the principle and system architecture of AF2 as well as the recipe of its success, and particularly focus on reviewing its applications in the fields of biology and medicine. Limitations of current AF2 prediction will also be discussed.
Collapse
Affiliation(s)
- Zhenyu Yang
- West China Biomedical Big Data Center, West China Hospital, Sichuan University, Chengdu, 610041, China
| | - Xiaoxi Zeng
- West China Biomedical Big Data Center, West China Hospital, Sichuan University, Chengdu, 610041, China.
| | - Yi Zhao
- West China Biomedical Big Data Center, West China Hospital, Sichuan University, Chengdu, 610041, China.
- Key Laboratory of Intelligent Information Processing, Advanced Computer Research Center, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, 100190, China.
| | - Runsheng Chen
- West China Biomedical Big Data Center, West China Hospital, Sichuan University, Chengdu, 610041, China.
- Key Laboratory of RNA Biology, Center for Big Data Research in Health, Institute of Biophysics, Chinese Academy of Sciences, Beijing, 100101, China.
- Pingshan Translational Medicine Center, Shenzhen Bay Laboratory, Shenzhen, 518118, China.
| |
Collapse
|
13
|
Marasco R, Fusi M, Coscolín C, Barozzi A, Almendral D, Bargiela R, Nutschel CGN, Pfleger C, Dittrich J, Gohlke H, Matesanz R, Sanchez-Carrillo S, Mapelli F, Chernikova TN, Golyshin PN, Ferrer M, Daffonchio D. Enzyme adaptation to habitat thermal legacy shapes the thermal plasticity of marine microbiomes. Nat Commun 2023; 14:1045. [PMID: 36828822 PMCID: PMC9958047 DOI: 10.1038/s41467-023-36610-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2021] [Accepted: 02/08/2023] [Indexed: 02/26/2023] Open
Abstract
Microbial communities respond to temperature with physiological adaptation and compositional turnover. Whether thermal selection of enzymes explains marine microbiome plasticity in response to temperature remains unresolved. By quantifying the thermal behaviour of seven functionally-independent enzyme classes (esterase, extradiol dioxygenase, phosphatase, beta-galactosidase, nuclease, transaminase, and aldo-keto reductase) in native proteomes of marine sediment microbiomes from the Irish Sea to the southern Red Sea, we record a significant effect of the mean annual temperature (MAT) on enzyme response in all cases. Activity and stability profiles of 228 esterases and 5 extradiol dioxygenases from sediment and seawater across 70 locations worldwide validate this thermal pattern. Modelling the esterase phase transition temperature as a measure of structural flexibility confirms the observed relationship with MAT. Furthermore, when considering temperature variability in sites with non-significantly different MATs, the broadest range of enzyme thermal behaviour and the highest growth plasticity of the enriched heterotrophic bacteria occur in samples with the widest annual thermal variability. These results indicate that temperature-driven enzyme selection shapes microbiome thermal plasticity and that thermal variability finely tunes such processes and should be considered alongside MAT in forecasting microbial community thermal response.
Collapse
Affiliation(s)
- Ramona Marasco
- Biological and Environmental Sciences and Engineering Division (BESE), Red Sea Research Centre (RSRC), King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia
| | - Marco Fusi
- Biological and Environmental Sciences and Engineering Division (BESE), Red Sea Research Centre (RSRC), King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia
- Centre for Conservation and Restoration Science, Edinburgh Napier University Sighthill Campus, Edinburgh, UK
| | | | - Alan Barozzi
- Biological and Environmental Sciences and Engineering Division (BESE), Red Sea Research Centre (RSRC), King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia
| | - David Almendral
- Instituto de Catalisis y Petroleoquimica (ICP), CSIC, Madrid, Spain
| | - Rafael Bargiela
- Centre for Environmental Biotechnology, School of Natural Sciences, Bangor University, Deiniol Rd, Bangor, UK
| | | | - Christopher Pfleger
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Jonas Dittrich
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Holger Gohlke
- Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, Jülich, Germany
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
- John von Neumann Institute for Computing (NIC) and Jülich Supercomputing Centre (JSC), Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Ruth Matesanz
- Spectroscopy Laboratory, Centro de Investigaciones Biologicas Margarita Salas (CIB), CSIC, Madrid, Spain
| | - Sergio Sanchez-Carrillo
- Instituto de Catalisis y Petroleoquimica (ICP), CSIC, Madrid, Spain
- Centro de Biologia Molecular Severo Ochoa (CBM), CSIC-UAM, Madrid, Spain
| | - Francesca Mapelli
- Department of Food Environmental and Nutritional Sciences, University of Milan, Milan, Italy
| | - Tatyana N Chernikova
- Centre for Environmental Biotechnology, School of Natural Sciences, Bangor University, Deiniol Rd, Bangor, UK
| | - Peter N Golyshin
- Centre for Environmental Biotechnology, School of Natural Sciences, Bangor University, Deiniol Rd, Bangor, UK
| | - Manuel Ferrer
- Instituto de Catalisis y Petroleoquimica (ICP), CSIC, Madrid, Spain.
| | - Daniele Daffonchio
- Biological and Environmental Sciences and Engineering Division (BESE), Red Sea Research Centre (RSRC), King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia.
| |
Collapse
|
14
|
Tucs A, Tsuda K, Sljoka A. Probing Conformational Dynamics of Antibodies with Geometric Simulations. Methods Mol Biol 2023; 2552:125-139. [PMID: 36346589 DOI: 10.1007/978-1-0716-2609-2_6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
This chapter describes the application of constrained geometric simulations for prediction of antibody structural dynamics. We utilize constrained geometric simulations method FRODAN, which is a low computational complexity alternative to molecular dynamics (MD) simulations that can rapidly explore flexible motions in protein structures. FRODAN is highly suited for conformational dynamics analysis of large proteins, complexes, intrinsically disordered proteins, and dynamics that occurs on longer biologically relevant time scales that are normally inaccessible to classical MD simulations. This approach predicts protein dynamics at an all-atom scale while retaining realistic covalent bonding, maintaining dihedral angles in energetically good conformations while avoiding steric clashes in addition to performing other geometric and stereochemical criteria checks. In this chapter, we apply FRODAN to showcase its applicability for probing functionally relevant dynamics of IgG2a, including large-amplitude domain-domain motions and motions of complementarity determining region (CDR) loops. As was suggested in previous experimental studies, our simulations show that antibodies can explore a large range of conformational space.
Collapse
Affiliation(s)
- Andrejs Tucs
- Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa, Chiba, Japan
| | - Koji Tsuda
- Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa, Chiba, Japan
- RIKEN Center for Advanced Intelligence Project, Tokyo, Japan
- Research and Services Division of Materials Data and Integrated System, National Institute for Materials Science, Tsukuba, Ibaraki, Japan
| | - Adnan Sljoka
- RIKEN Center for Advanced Intelligence Project, Tokyo, Japan.
- Department of Chemistry, York University, Toronto, Canada.
| |
Collapse
|
15
|
Baksh KA, Augustine J, Sljoka A, Prosser RS, Zamble DB. Mechanistic insights into the nickel-dependent allosteric response of the Helicobacter pylori NikR transcription factor. J Biol Chem 2022; 299:102785. [PMID: 36502919 PMCID: PMC9860126 DOI: 10.1016/j.jbc.2022.102785] [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/03/2022] [Revised: 11/30/2022] [Accepted: 12/03/2022] [Indexed: 12/13/2022] Open
Abstract
In Helicobacter pylori, the nickel-responsive NikR transcription factor plays a key role in regulating intracellular nickel concentrations, which is an essential process for survival of this pathogen in the acidic human stomach. Nickel binding to H. pylori NikR (HpNikR) allosterically activates DNA binding to target promoters encoding genes involved in nickel homeostasis and acid adaptation, to either activate or repress their transcription. We previously showed that HpNikR adopts an equilibrium between an open conformation and DNA-binding competent cis and trans states. Nickel binding slows down conformational exchange between these states and shifts the equilibrium toward the binding-competent states. The protein then becomes stabilized in a cis conformation upon binding the ureA promoter. Here, we investigate how nickel binding creates this response and how it is transmitted to the DNA-binding domains. Through mutagenesis, DNA-binding studies, and computational methods, the allosteric response to nickel was found to be propagated from the nickel-binding sites to the DNA-binding domains via the β-sheets of the metal-binding domain and a network of residues at the inter-domain interface. Our computational results suggest that nickel binding increases protein rigidity to slow down the conformational exchange. A thymine base in the ureA promoter sequence, known to be critical for high affinity DNA binding by HpNikR, was also found to be important for the allosteric response, while a modified version of this promoter further highlighted the importance of the DNA sequence in modulating the response. Collectively, our results provide insights into regulation of a key protein for H. pylori survival.
Collapse
Affiliation(s)
- Karina A. Baksh
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
| | - Jerry Augustine
- Department of Chemistry, University of Toronto, Toronto, Ontario, Canada
| | - Adnan Sljoka
- RIKEN Center for Advanced Intelligence Project, RIKEN, Chuo-ku, Tokyo, Japan,For correspondence: R. Scott Prosser; Adnan Sljoka
| | - R. Scott Prosser
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada,Department of Chemistry, University of Toronto, Toronto, Ontario, Canada,For correspondence: R. Scott Prosser; Adnan Sljoka
| | - Deborah B. Zamble
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada,Department of Chemistry, University of Toronto, Toronto, Ontario, Canada
| |
Collapse
|
16
|
Tanaka H, Seki Y, Ueno S, Shibutani Y. Conformational deformation of a multi-jointed elastic loop. Sci Rep 2022; 12:19984. [PMID: 36411324 PMCID: PMC9678883 DOI: 10.1038/s41598-022-24355-7] [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: 06/21/2022] [Accepted: 11/14/2022] [Indexed: 11/23/2022] Open
Abstract
A new class of deformation is presented for a planar loop structure made up of slender elastic bodies and joints. In demonstrating the circumferential shortening of the multi-jointed elastic loop, diverse three-dimensional (3D) deformations emerge through piecewise deflections and discrete rotations. These 3D morphologies correspond to conformations of molecular ring systems. Through image processing, the 3D reconstructions of the deformed structures are characterized by number, geometry, and initial imperfections of the body segments. We elucidate from measurements that the conformational deformation without self-stress results from a cyclical assembly of compressive bending of elastic bodies with high shear rigidity. The mechanical insights gained may apply in controlling the polymorphism exhibited by the cyclical structures across scales.
Collapse
Affiliation(s)
- Hiro Tanaka
- grid.136593.b0000 0004 0373 3971Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871 Japan
| | - Yuji Seki
- grid.136593.b0000 0004 0373 3971Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871 Japan
| | - Shohei Ueno
- grid.136593.b0000 0004 0373 3971Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871 Japan
| | - Yoji Shibutani
- grid.136593.b0000 0004 0373 3971Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871 Japan ,grid.267852.c0000 0004 0637 2083Nanotechnology Program, VNU Vietnam Japan University, Luu Huu Phuoc Street, My Dinh 1 Ward, Nam Tu Liem District, Hanoi, Vietnam
| |
Collapse
|
17
|
Berry M, Limberg D, Lee-Trimble ME, Hayward R, Santangelo CD. Controlling the configuration space topology of mechanical structures. Phys Rev E 2022; 106:055002. [PMID: 36559440 DOI: 10.1103/physreve.106.055002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2022] [Accepted: 09/15/2022] [Indexed: 06/17/2023]
Abstract
Linkages are mechanical devices constructed from rigid bars and freely rotating joints studied both for their utility in engineering and as mathematical idealizations in a number of physical systems. Recently, there has been a resurgence of interest in designing linkages in the physics community due to the concurrent developments of mechanical metamaterials, topological mechanics, and the discovery of anomalous rigidity in fiber networks and vertex models. These developments raise a natural question: to what extent can the motion of a linkage or mechanical structure be designed? Here, we describe a method to design the topology of the configuration space of a linkage by first identifying the manifold of critical points, then perturbing around such critical configurations. Unlike other methods, our methods are tractable and provide a simple visual toolkit for mechanism design. We demonstrate our procedure by designing a mechanism to gate the propagation of a soliton in a Kane-Lubensky chain of interconnected rotors.
Collapse
Affiliation(s)
- M Berry
- Department of Physics, Syracuse University, Syracuse, New York 13244, USA
| | - David Limberg
- Department of Polymer Science and Engineering, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - M E Lee-Trimble
- Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Ryan Hayward
- Department of Chemical and Biological Engineering, University of Colorado, Boulder, Colorado 80309, USA
| | - C D Santangelo
- Department of Physics, Syracuse University, Syracuse, New York 13244, USA
| |
Collapse
|
18
|
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
|
19
|
Pauly T, Bolakhrif N, Kaiser J, Nagel-Steger L, Gremer L, Gohlke H, Willbold D. Met/Val129 polymorphism of the full-length human prion protein dictates distinct pathways of amyloid formation. J Biol Chem 2022; 298:102430. [PMID: 36037966 PMCID: PMC9513279 DOI: 10.1016/j.jbc.2022.102430] [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: 05/17/2022] [Revised: 08/18/2022] [Accepted: 08/22/2022] [Indexed: 11/25/2022] Open
Abstract
Methionine/valine polymorphism at position 129 of the human prion protein, huPrP, is tightly associated with the pathogenic phenotype, disease progress, and age of onset of neurodegenerative diseases such as Creutzfeldt–Jakob disease or Fatal Familial Insomnia. This raises the question of whether and how the amino acid type at position 129 influences the structural properties of huPrP, affecting its folding, stability, and amyloid formation behavior. Here, our detailed biophysical characterization of the 129M and 129V variants of recombinant full-length huPrP(23–230) by amyloid formation kinetics, CD spectroscopy, molecular dynamics simulations, and sedimentation velocity analysis reveals differences in their aggregation propensity and oligomer content, leading to deviating pathways for the conversion into amyloid at acidic pH. We determined that the 129M variant exhibits less secondary structure content before amyloid formation and higher resistance to thermal denaturation compared to the 129V variant, whereas the amyloid conformation of both variants shows similar thermal stability. Additionally, our molecular dynamics simulations and rigidity analyses at the atomistic level identify intramolecular interactions responsible for the enhanced monomer stability of the 129M variant, involving more frequent minimum distances between E196 and R156, forming a salt bridge. Removal of the N-terminal half of the 129M full-length variant diminishes its differences compared to the 129V full-length variant and highlights the relevance of the flexible N terminus in huPrP. Taken together, our findings provide insight into structural properties of huPrP and the effects of the amino acid identity at position 129 on amyloid formation behavior.
Collapse
Affiliation(s)
- Thomas Pauly
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry) and JuStruct: Jülich Center of Structural Biology, Forschungszentrum Jülich, 52425 Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
| | - Najoua Bolakhrif
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry) and JuStruct: Jülich Center of Structural Biology, Forschungszentrum Jülich, 52425 Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
| | - Jesko Kaiser
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
| | - Luitgard Nagel-Steger
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry) and JuStruct: Jülich Center of Structural Biology, Forschungszentrum Jülich, 52425 Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
| | - Lothar Gremer
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry) and JuStruct: Jülich Center of Structural Biology, Forschungszentrum Jülich, 52425 Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany; Jülich Supercomputing Centre (JSC), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany; Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany; John von Neumann Institute for Computing (NIC), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Dieter Willbold
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry) and JuStruct: Jülich Center of Structural Biology, Forschungszentrum Jülich, 52425 Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany.
| |
Collapse
|
20
|
Yüksel S, Bonus M, Schwabe T, Pfleger C, Zimmer T, Enke U, Saß I, Gohlke H, Benndorf K, Kusch J. Uncoupling of Voltage- and Ligand-Induced Activation in HCN2 Channels by Glycine Inserts. Front Physiol 2022; 13:895324. [PMID: 36091400 PMCID: PMC9452628 DOI: 10.3389/fphys.2022.895324] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 06/20/2022] [Indexed: 11/24/2022] Open
Abstract
Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are tetramers that generate electrical rhythmicity in special brain neurons and cardiomyocytes. The channels are activated by membrane hyperpolarization. The binding of cAMP to the four available cyclic nucleotide-binding domains (CNBD) enhances channel activation. We analyzed in the present study the mechanism of how the effect of cAMP binding is transmitted to the pore domain. Our strategy was to uncouple the C-linker (CL) from the channel core by inserting one to five glycine residues between the S6 gate and the A′-helix (constructs 1G to 5G). We quantified in full-length HCN2 channels the resulting functional effects of the inserted glycines by current activation as well as the structural dynamics and statics using molecular dynamics simulations and Constraint Network Analysis. We show functionally that already in 1G the cAMP effect on activation is lost and that with the exception of 3G and 5G the concentration-activation relationships are shifted to depolarized voltages with respect to HCN2. The strongest effect was found for 4G. Accordingly, the activation kinetics were accelerated by all constructs, again with the strongest effect in 4G. The simulations reveal that the average residue mobility of the CL and CNBD domains is increased in all constructs and that the junction between the S6 and A′-helix is turned into a flexible hinge, resulting in a destabilized gate in all constructs. Moreover, for 3G and 4G, there is a stronger downward displacement of the CL-CNBD than in HCN2 and the other constructs, resulting in an increased kink angle between S6 and A′-helix, which in turn loosens contacts between the S4-helix and the CL. This is suggested to promote a downward movement of the S4-helix, similar to the effect of hyperpolarization. In addition, exclusively in 4G, the selectivity filter in the upper pore region and parts of the S4-helix are destabilized. The results provide new insights into the intricate activation of HCN2 channels.
Collapse
Affiliation(s)
- Sezin Yüksel
- Universitätsklinikum Jena, Institut für Physiologie II, Jena, Germany
| | - Michele Bonus
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Tina Schwabe
- Universitätsklinikum Jena, Institut für Physiologie II, Jena, Germany
| | - Christopher Pfleger
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Thomas Zimmer
- Universitätsklinikum Jena, Institut für Physiologie II, Jena, Germany
| | - Uta Enke
- Universitätsklinikum Jena, Institut für Physiologie II, Jena, Germany
| | - Inga Saß
- Universitätsklinikum Jena, Institut für Physiologie II, Jena, Germany
| | - Holger Gohlke
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry) and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, Jülich, Germany
- *Correspondence: Holger Gohlke, ; Klaus Benndorf, ; Jana Kusch,
| | - Klaus Benndorf
- Universitätsklinikum Jena, Institut für Physiologie II, Jena, Germany
- *Correspondence: Holger Gohlke, ; Klaus Benndorf, ; Jana Kusch,
| | - Jana Kusch
- Universitätsklinikum Jena, Institut für Physiologie II, Jena, Germany
- *Correspondence: Holger Gohlke, ; Klaus Benndorf, ; Jana Kusch,
| |
Collapse
|
21
|
Low pH structure of heliorhodopsin reveals chloride binding site and intramolecular signaling pathway. Sci Rep 2022; 12:13955. [PMID: 35977989 PMCID: PMC9385722 DOI: 10.1038/s41598-022-17716-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Accepted: 07/29/2022] [Indexed: 11/11/2022] Open
Abstract
Within the microbial rhodopsin family, heliorhodopsins (HeRs) form a phylogenetically distinct group of light-harvesting retinal proteins with largely unknown functions. We have determined the 1.97 Å resolution X-ray crystal structure of Thermoplasmatales archaeon SG8-52-1 heliorhodopsin (TaHeR) in the presence of NaCl under acidic conditions (pH 4.5), which complements the known 2.4 Å TaHeR structure acquired at pH 8.0. The low pH structure revealed that the hydrophilic Schiff base cavity (SBC) accommodates a chloride anion to stabilize the protonated retinal Schiff base when its primary counterion (Glu-108) is neutralized. Comparison of the two structures at different pH revealed conformational changes connecting the SBC and the extracellular loop linking helices A–B. We corroborated this intramolecular signaling transduction pathway with computational studies, which revealed allosteric network changes propagating from the perturbed SBC to the intracellular and extracellular space, suggesting TaHeR may function as a sensory rhodopsin. This intramolecular signaling mechanism may be conserved among HeRs, as similar changes were observed for HeR 48C12 between its pH 8.8 and pH 4.3 structures. We additionally performed DEER experiments, which suggests that TaHeR forms possible dimer-of-dimer associations which may be integral to its putative functionality as a light sensor in binding a transducer protein.
Collapse
|
22
|
Chen S, Giardina F, Choi GPT, Mahadevan L. Modular representation and control of floppy networks. Proc Math Phys Eng Sci 2022. [DOI: 10.1098/rspa.2022.0082] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Geometric graph models of systems as diverse as proteins, DNA assemblies, architected materials and robot swarms are useful abstract representations of these objects that also unify ways to study their properties and control them in space and time. While much work has been done in the context of characterizing the behaviour of these networks close to critical points associated with bond and rigidity percolation, isostaticity, etc., much less is known about floppy, underconstrained networks that are far more common in nature and technology. Here, we combine geometric rigidity and algebraic sparsity to provide a framework for identifying the zero energy floppy modes via a representation that illuminates the underlying hierarchy and modularity of the network and thence the control of its nestedness and locality. Our framework allows us to demonstrate a range of applications of this approach that include robotic reaching tasks with motion primitives, and predicting the linear and nonlinear response of elastic networks based solely on infinitesimal rigidity and sparsity, which we test using physical experiments. Our approach is thus likely to be of use broadly in dissecting the geometrical properties of floppy networks using algebraic sparsity to optimize their function and performance.
Collapse
Affiliation(s)
- Siheng Chen
- School of Engineering and Applied Sciences, Harvard University, Cambridge MA 02138, USA
| | - Fabio Giardina
- School of Engineering and Applied Sciences, Harvard University, Cambridge MA 02138, USA
| | - Gary P. T. Choi
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge MA 02139, USA
| | - L. Mahadevan
- School of Engineering and Applied Sciences, Harvard University, Cambridge MA 02138, USA
- Department of Physics, Harvard University, Cambridge MA 02138, USA
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge MA 02138, USA
| |
Collapse
|
23
|
Sieg J, Sandmeier CC, Lieske J, Meents A, Lemmen C, Streit WR, Rarey M. Analyzing structural features of proteins from deep-sea organisms. Proteins 2022; 90:1521-1537. [PMID: 35313380 DOI: 10.1002/prot.26337] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 03/10/2022] [Accepted: 03/15/2022] [Indexed: 12/31/2022]
Abstract
Protein adaptations to extreme environmental conditions are drivers in biotechnological process optimization and essential to unravel the molecular limits of life. Most proteins with such desirable adaptations are found in extremophilic organisms inhabiting extreme environments. The deep sea is such an environment and a promising resource that poses multiple extremes on its inhabitants. Conditions like high hydrostatic pressure and high or low temperature are prevalent and many deep-sea organisms tolerate multiple of these extremes. While molecular adaptations to high temperature are comparatively good described, adaptations to other extremes like high pressure are not well-understood yet. To fully unravel the molecular mechanisms of individual adaptations it is probably necessary to disentangle multifactorial adaptations. In this study, we evaluate differences of protein structures from deep-sea organisms and their respective related proteins from nondeep-sea organisms. We created a data collection of 1281 experimental protein structures from 25 deep-sea organisms and paired them with orthologous proteins. We exhaustively evaluate differences between the protein pairs with machine learning and Shapley values to determine characteristic differences in sequence and structure. The results show a reasonable discrimination of deep-sea and nondeep-sea proteins from which we distinguish correlations previously attributed to thermal stability from other signals potentially describing adaptions to high pressure. While some distinct correlations can be observed the overall picture appears intricate.
Collapse
Affiliation(s)
- Jochen Sieg
- Universität Hamburg, ZBH - Center for Bioinformatics, Hamburg, Germany
| | | | - Julia Lieske
- Deutsches Elektronen-Synchrotron DESY, Center for Free-Electron Laser Science, Hamburg, Germany
| | - Alke Meents
- Deutsches Elektronen-Synchrotron DESY, Center for Free-Electron Laser Science, Hamburg, Germany
| | | | - Wolfgang R Streit
- Universität Hamburg, Department of Microbiology and Biotechnology, Hamburg, Germany
| | - Matthias Rarey
- Universität Hamburg, ZBH - Center for Bioinformatics, Hamburg, Germany
| |
Collapse
|
24
|
Heimsch KC, Gertzen CGW, Schuh AK, Nietzel T, Rahlfs S, Przyborski JM, Gohlke H, Schwarzländer M, Becker K, Fritz-Wolf K. Structure and Function of Redox-Sensitive Superfolder Green Fluorescent Protein Variant. Antioxid Redox Signal 2022; 37:1-18. [PMID: 35072524 PMCID: PMC9293687 DOI: 10.1089/ars.2021.0234] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Aims: Genetically encoded green fluorescent protein (GFP)-based redox biosensors are widely used to monitor specific and dynamic redox processes in living cells. Over the last few years, various biosensors for a variety of applications were engineered and enhanced to match the organism and cellular environments, which should be investigated. In this context, the unicellular intraerythrocytic parasite Plasmodium, the causative agent of malaria, represents a challenge, as the small size of the organism results in weak fluorescence signals that complicate precise measurements, especially for cell compartment-specific observations. To address this, we have functionally and structurally characterized an enhanced redox biosensor superfolder roGFP2 (sfroGFP2). Results: SfroGFP2 retains roGFP2-like behavior, yet with improved fluorescence intensity (FI) in cellulo. SfroGFP2-based redox biosensors are pH insensitive in a physiological pH range and show midpoint potentials comparable with roGFP2-based redox biosensors. Using crystallography and rigidity theory, we identified the superfolding mutations as being responsible for improved structural stability of the biosensor in a redox-sensitive environment, thus explaining the improved FI in cellulo. Innovation: This work provides insight into the structure and function of GFP-based redox biosensors. It describes an improved redox biosensor (sfroGFP2) suitable for measuring oxidizing effects within small cells where applicability of other redox sensor variants is limited. Conclusion: Improved structural stability of sfroGFP2 gives rise to increased FI in cellulo. Fusion to hGrx1 (human glutaredoxin-1) provides the hitherto most suitable biosensor for measuring oxidizing effects in Plasmodium. This sensor is of major interest for studying glutathione redox changes in small cells, as well as subcellular compartments in general. Antioxid. Redox Signal. 37, 1-18.
Collapse
Affiliation(s)
- Kim C Heimsch
- Biochemistry and Molecular Biology, Interdisciplinary Research Center, Justus Liebig University Giessen, Giessen, Germany
| | - Christoph G W Gertzen
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany.,Center for Structural Studies (CSS), Heinrich Heine University Düsseldorf, Düsseldorf, Germany
| | - Anna Katharina Schuh
- Biochemistry and Molecular Biology, Interdisciplinary Research Center, Justus Liebig University Giessen, Giessen, Germany
| | - Thomas Nietzel
- Institute of Plant Biology and Biotechnology, University of Münster, Münster, Germany
| | - Stefan Rahlfs
- Biochemistry and Molecular Biology, Interdisciplinary Research Center, Justus Liebig University Giessen, Giessen, Germany
| | - Jude M Przyborski
- Biochemistry and Molecular Biology, Interdisciplinary Research Center, Justus Liebig University Giessen, Giessen, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany.,John von Neumann Institute of Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Markus Schwarzländer
- Institute of Plant Biology and Biotechnology, University of Münster, Münster, Germany
| | - Katja Becker
- Biochemistry and Molecular Biology, Interdisciplinary Research Center, Justus Liebig University Giessen, Giessen, Germany
| | - Karin Fritz-Wolf
- Biochemistry and Molecular Biology, Interdisciplinary Research Center, Justus Liebig University Giessen, Giessen, Germany.,Max-Planck Institute of Medical Research, Heidelberg, Germany
| |
Collapse
|
25
|
Ullanat V, Kasukurthi N, Viswanath S. PrISM: Precision for Integrative Structural Models. Bioinformatics 2022; 38:3837-3839. [PMID: 35723541 DOI: 10.1093/bioinformatics/btac400] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2022] [Revised: 06/04/2022] [Accepted: 06/16/2022] [Indexed: 11/14/2022] Open
Abstract
MOTIVATION A single precision value is currently reported for an integrative model. However, precision may vary for different regions of an integrative model owing to varying amounts of input information. RESULTS We develop PrISM (Precision for Integrative Structural Models), to efficiently identify high and low-precision regions for integrative models. AVAILABILITY PrISM is written in Python and available under the GNU General Public License v3.0 at https://github.com/isblab/prism; benchmark data used in this paper is available at doi:10.5281/zenodo.6241200. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Varun Ullanat
- National Center for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India
| | - Nikhil Kasukurthi
- National Center for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India
| | - Shruthi Viswanath
- National Center for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India
| |
Collapse
|
26
|
Malik FK, Guo JT. Insights into protein-DNA interactions from hydrogen bond energy-based comparative protein-ligand analyses. Proteins 2022; 90:1303-1314. [PMID: 35122321 PMCID: PMC9018545 DOI: 10.1002/prot.26313] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2021] [Revised: 01/17/2022] [Accepted: 01/31/2022] [Indexed: 01/18/2023]
Abstract
Hydrogen bonds play important roles in protein folding and protein-ligand interactions, particularly in specific protein-DNA recognition. However, the distributions of hydrogen bonds, especially hydrogen bond energy (HBE) in different types of protein-ligand complexes, is unknown. Here we performed a comparative analysis of hydrogen bonds among three non-redundant datasets of protein-protein, protein-peptide, and protein-DNA complexes. Besides comparing the number of hydrogen bonds in terms of types and locations, we investigated the distributions of HBE. Our results indicate that while there is no significant difference of hydrogen bonds within protein chains among the three types of complexes, interfacial hydrogen bonds are significantly more prevalent in protein-DNA complexes. More importantly, the interfacial hydrogen bonds in protein-DNA complexes displayed a unique energy distribution of strong and weak hydrogen bonds whereas majority of the interfacial hydrogen bonds in protein-protein and protein-peptide complexes are of predominantly high strength with low energy. Moreover, there is a significant difference in the energy distributions of minor groove hydrogen bonds between protein-DNA complexes with different binding specificity. Highly specific protein-DNA complexes contain more strong hydrogen bonds in the minor groove than multi-specific complexes, suggesting important role of minor groove in specific protein-DNA recognition. These results can help better understand protein-DNA interactions and have important implications in improving quality assessments of protein-DNA complex models.
Collapse
Affiliation(s)
- Fareeha K Malik
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, North Carolina, USA.,Research Center of Modeling and Simulation, National University of Science and Technology, Islamabad, Pakistan
| | - Jun-Tao Guo
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, North Carolina, USA
| |
Collapse
|
27
|
Lata S, Akif M. Probing structural basis for enhanced binding of SARS-CoV-2 P.1 variant spike protein with the human ACE2 receptor. J Cell Biochem 2022; 123:1207-1221. [PMID: 35620980 PMCID: PMC9347910 DOI: 10.1002/jcb.30276] [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: 01/21/2022] [Revised: 04/22/2022] [Accepted: 05/10/2022] [Indexed: 11/25/2022]
Abstract
The initial step of infection by severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) involves the binding of receptor binding domain (RBD) of the spike protein to the angiotensin converting enzyme 2 (ACE2) receptor. Each successive wave of SARS‐CoV‐2 reports emergence of many new variants, which is associated with mutations in the RBD as well as other parts of the spike protein. These mutations are reported to have enhanced affinity towards the ACE2 receptor as well as are also crucial for the virus transmission. Many computational and experimental studies have demonstrated the effect of individual mutation on the RBD‐ACE2 binding. However, the cumulative effect of mutations on the RBD and away from the RBD was not investigated in detail. We report here a comparative analysis on the structural communication and dynamics of the RBD and truncated S1 domain of spike protein in complex with the ACE2 receptor from SARS‐CoV‐2 wild type and its P.1 variant. Our integrative network and dynamics approaches highlighted a subtle conformational changes in the RBD as well as truncated S1 domain of spike protein at the protein contact level, responsible for the increased affinity with the ACE2 receptor. Moreover, our study also identified the commonalities and differences in the dynamics of the interactions between spike protein of SARS‐CoV‐2 wild type and its P.1 variant with the ACE2 receptor. Further, our investigation yielded an understanding towards identification of the unique RBD residues crucial for the interaction with the ACE2 host receptor. Overall, the study provides an insight for designing better therapeutics against the circulating P.1 variants as well as other future variants.
Collapse
Affiliation(s)
- Surabhi Lata
- Laboratory of Structural Biology, Department of Biochemistry, School of Life Sciences, University of Hyderabad, Hyderabad, Telangana, India
| | - Mohd Akif
- Laboratory of Structural Biology, Department of Biochemistry, School of Life Sciences, University of Hyderabad, Hyderabad, Telangana, India
| |
Collapse
|
28
|
The accuracy of protein structures in solution determined by AlphaFold and NMR. Structure 2022; 30:925-933.e2. [DOI: 10.1016/j.str.2022.04.005] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2022] [Revised: 03/18/2022] [Accepted: 04/13/2022] [Indexed: 02/05/2023]
|
29
|
Light- and pH-dependent structural changes in cyanobacteriochrome AnPixJg2. Photochem Photobiol Sci 2022; 21:447-469. [PMID: 35394641 DOI: 10.1007/s43630-022-00204-4] [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: 12/22/2021] [Accepted: 03/03/2022] [Indexed: 10/18/2022]
Abstract
Cyanobacteriochromes (CBCRs) are phytochrome-related photosensory proteins that play an essential role in regulating phototaxis, chromatic acclimation, and cell aggregation in cyanobacteria. Here, we apply solid-state NMR spectroscopy to the red/green GAF2 domain of the CBCR AnPixJ assembled in vitro with a uniformly 13C- and 15N-labeled bilin chromophore, tracking changes in electronic structure, geometry, and structural heterogeneity of the chromophore as well as intimate contacts between the chromophore and protein residues in the photocycle. Our data confirm that the bilin ring D is strongly twisted with respect to the B-C plane in both dark and photoproduct states. We also identify a greater structural heterogeneity of the bilin chromophore in the photoproduct than in the dark state. In addition, the binding pocket is more hydrated in the photoproduct. Observation of interfacial 1H contacts of the photoproduct chromophore, together with quantum mechanics/molecular mechanics (QM/MM)-based structural models for this photoproduct, clearly suggests the presence of a biprotonated (cationic) imidazolium side-chain for a conserved histidine residue (322) at a distance of ~2.7 Å, generalizing the recent theoretical findings that explicitly link the structural heterogeneity of the dark-state chromophore to the protonation of this specific residue. Moreover, we examine pH effects on this in vitro assembled holoprotein, showing a substantially altered electronic structure and protonation of the photoproduct chromophore even with a small pH drop from 7.8 to 7.2. Our studies provide further information regarding the light- and pH-induced changes of the chromophore and the rearrangements of the hydrogen-bonding and electrostatic interaction network around it. Possible correlations between structural heterogeneity of the chromophore, protonation of the histidine residue nearby, and hydration of the pocket in both photostates are discussed.
Collapse
|
30
|
Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex. Molecules 2022; 27:molecules27030799. [PMID: 35164065 PMCID: PMC8839245 DOI: 10.3390/molecules27030799] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 01/20/2022] [Accepted: 01/21/2022] [Indexed: 12/02/2022] Open
Abstract
The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.
Collapse
|
31
|
Hüdig M, Tronconi MA, Zubimendi JP, Sage TL, Poschmann G, Bickel D, Gohlke H, Maurino VG. Respiratory and C4-photosynthetic NAD-malic enzyme coexist in bundle sheath cell mitochondria and evolved via association of differentially adapted subunits. THE PLANT CELL 2022; 34:597-615. [PMID: 34734993 PMCID: PMC8773993 DOI: 10.1093/plcell/koab265] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Accepted: 10/26/2021] [Indexed: 05/29/2023]
Abstract
In plant mitochondria, nicotinamide adenine dinucleotide-malic enzyme (NAD-ME) has a housekeeping function in malate respiration. In different plant lineages, NAD-ME was independently co-opted in C4 photosynthesis. In the C4 Cleome species, Gynandropsis gynandra and Cleome angustifolia, all NAD-ME genes (NAD-MEα, NAD-MEβ1, and NAD-MEβ2) were affected by C4 evolution and are expressed at higher levels than their orthologs in the C3 species Tarenaya hassleriana. In T. hassleriana, the NAD-ME housekeeping function is performed by two heteromers, NAD-MEα/β1 and NAD-MEα/β2, with similar biochemical properties. In both C4 species, this role is restricted to NAD-MEα/β2. In the C4 species, NAD-MEα/β1 is exclusively present in the leaves, where it accounts for most of the enzymatic activity. Gynandropsis gynandra NAD-MEα/β1 (GgNAD-MEα/β1) exhibits high catalytic efficiency and is differentially activated by the C4 intermediate aspartate, confirming its role as the C4-decarboxylase. During C4 evolution, NAD-MEβ1 lost its catalytic activity; its contribution to the enzymatic activity results from a stabilizing effect on the associated α-subunit and the acquisition of regulatory properties. We conclude that in bundle sheath cell mitochondria of C4 species, the functions of NAD-ME as C4 photosynthetic decarboxylase and as a housekeeping enzyme coexist and are performed by isoforms that combine the same α-subunit with differentially adapted β-subunits.
Collapse
Affiliation(s)
- Meike Hüdig
- Molekulare Pflanzenphysiologie, Rheinische Friedrich-Wilhelms-Universität Bonn, Kirschallee, Bonn 53115, Germany
| | - Marcos A Tronconi
- Centro de Estudios Fotosintéticos y Bioquímicos, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, Suipacha 531, Rosario 2000, Argentina
| | - Juan P Zubimendi
- Centro de Estudios Fotosintéticos y Bioquímicos, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, Suipacha 531, Rosario 2000, Argentina
| | - Tammy L Sage
- Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada
| | - Gereon Poschmann
- Molecular Proteomics Laboratory, Biomedical Research Centre (BMFZ) & Institute of Molecular Medicine, Proteome Research, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf, Universitätsstraße 1, Düsseldorf 40225, Germany
| | - David Bickel
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Universitätsstraße 1, Düsseldorf 40225, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Universitätsstraße 1, Düsseldorf 40225, Germany
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry) & Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, Jülich 52425, Germany
| | - Veronica G Maurino
- Molekulare Pflanzenphysiologie, Rheinische Friedrich-Wilhelms-Universität Bonn, Kirschallee, Bonn 53115, Germany
| |
Collapse
|
32
|
Huang SK, Almurad O, Pejana RJ, Morrison ZA, Pandey A, Picard LP, Nitz M, Sljoka A, Prosser RS. Allosteric modulation of the adenosine A 2A receptor by cholesterol. eLife 2022; 11:e73901. [PMID: 34986091 PMCID: PMC8730723 DOI: 10.7554/elife.73901] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Accepted: 12/23/2021] [Indexed: 12/18/2022] Open
Abstract
Cholesterol is a major component of the cell membrane and commonly regulates membrane protein function. Here, we investigate how cholesterol modulates the conformational equilibria and signaling of the adenosine A2A receptor (A2AR) in reconstituted phospholipid nanodiscs. This model system conveniently excludes possible effects arising from cholesterol-induced phase separation or receptor oligomerization and focuses on the question of allostery. GTP hydrolysis assays show that cholesterol weakly enhances the basal signaling of A2AR while decreasing the agonist EC50. Fluorine nuclear magnetic resonance (19F NMR) spectroscopy shows that this enhancement arises from an increase in the receptor's active state population and a G-protein-bound precoupled state. 19F NMR of fluorinated cholesterol analogs reveals transient interactions with A2AR, indicating a lack of high-affinity binding or direct allosteric modulation. The combined results suggest that the observed allosteric effects are largely indirect and originate from cholesterol-mediated changes in membrane properties, as shown by membrane fluidity measurements and high-pressure NMR.
Collapse
Affiliation(s)
- Shuya Kate Huang
- Department of Chemistry, University of TorontoTorontoCanada
- Department of Chemical and Physical Sciences, University of Toronto MississaugaMississaugaCanada
| | - Omar Almurad
- Department of Chemistry, University of TorontoTorontoCanada
- Department of Chemical and Physical Sciences, University of Toronto MississaugaMississaugaCanada
| | - Reizel J Pejana
- Department of Chemistry, University of TorontoTorontoCanada
- Department of Chemical and Physical Sciences, University of Toronto MississaugaMississaugaCanada
| | | | - Aditya Pandey
- Department of Chemistry, University of TorontoTorontoCanada
- Department of Chemical and Physical Sciences, University of Toronto MississaugaMississaugaCanada
| | - Louis-Philippe Picard
- Department of Chemistry, University of TorontoTorontoCanada
- Department of Chemical and Physical Sciences, University of Toronto MississaugaMississaugaCanada
| | - Mark Nitz
- Department of Chemistry, University of TorontoTorontoCanada
| | - Adnan Sljoka
- RIKEN Center for Advanced Intelligence ProjectTokyoJapan
- York University, Department of ChemistryTorontoCanada
| | - R Scott Prosser
- Department of Chemistry, University of TorontoTorontoCanada
- Department of Chemical and Physical Sciences, University of Toronto MississaugaMississaugaCanada
- Department of Biochemistry, University of TorontoTorontoCanada
| |
Collapse
|
33
|
Lai J, Yang J, Gamsiz Uzun ED, Rubenstein BM, Sarkar IN. LYRUS: a machine learning model for predicting the pathogenicity of missense variants. BIOINFORMATICS ADVANCES 2021; 2:vbab045. [PMID: 35036922 PMCID: PMC8754197 DOI: 10.1093/bioadv/vbab045] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/08/2021] [Revised: 12/08/2021] [Accepted: 12/21/2021] [Indexed: 01/27/2023]
Abstract
SUMMARY Single amino acid variations (SAVs) are a primary contributor to variations in the human genome. Identifying pathogenic SAVs can provide insights to the genetic architecture of complex diseases. Most approaches for predicting the functional effects or pathogenicity of SAVs rely on either sequence or structural information. This study presents 〈Lai Yang Rubenstein Uzun Sarkar〉 (LYRUS), a machine learning method that uses an XGBoost classifier to predict the pathogenicity of SAVs. LYRUS incorporates five sequence-based, six structure-based and four dynamics-based features. Uniquely, LYRUS includes a newly proposed sequence co-evolution feature called the variation number. LYRUS was trained using a dataset that contains 4363 protein structures corresponding to 22 639 SAVs from the ClinVar database, and tested using the VariBench testing dataset. Performance analysis showed that LYRUS achieved comparable performance to current variant effect predictors. LYRUS's performance was also benchmarked against six Deep Mutational Scanning datasets for PTEN and TP53. AVAILABILITY AND IMPLEMENTATION LYRUS is freely available and the source code can be found at https://github.com/jiaying2508/LYRUS. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics Advances online.
Collapse
Affiliation(s)
- Jiaying Lai
- Center for Biomedical Informatics, Brown University, Providence, RI 02903, USA,Center for Computational Molecular Biology, Brown University, Providence, RI 02906, USA
| | - Jordan Yang
- Department of Chemistry, Brown University, Providence, RI 02906, USA
| | - Ece D Gamsiz Uzun
- Center for Computational Molecular Biology, Brown University, Providence, RI 02906, USA,Department of Pathology and Laboratory Medicine, Brown University Alpert Medical School, Providence, RI 02903, USA,Department of Pathology, Rhode Island Hospital, Providence, RI 02903, USA
| | - Brenda M Rubenstein
- Center for Computational Molecular Biology, Brown University, Providence, RI 02906, USA,Department of Chemistry, Brown University, Providence, RI 02906, USA,To whom correspondence should be addressed. and
| | - Indra Neil Sarkar
- Center for Biomedical Informatics, Brown University, Providence, RI 02903, USA,Rhode Island Quality Institute, Providence, RI 02908, USA,To whom correspondence should be addressed. and
| |
Collapse
|
34
|
Schwarz D, Georges G, Kelm S, Shi J, Vangone A, Deane CM. Co-evolutionary distance predictions contain flexibility information. Bioinformatics 2021; 38:65-72. [PMID: 34383892 DOI: 10.1093/bioinformatics/btab562] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Revised: 06/19/2021] [Accepted: 08/10/2021] [Indexed: 02/03/2023] Open
Abstract
MOTIVATION Co-evolution analysis can be used to accurately predict residue-residue contacts from multiple sequence alignments. The introduction of machine-learning techniques has enabled substantial improvements in precision and a shift from predicting binary contacts to predict distances between pairs of residues. These developments have significantly improved the accuracy of de novo prediction of static protein structures. With AlphaFold2 lifting the accuracy of some predicted protein models close to experimental levels, structure prediction research will move on to other challenges. One of those areas is the prediction of more than one conformation of a protein. Here, we examine the potential of residue-residue distance predictions to be informative of protein flexibility rather than simply static structure. RESULTS We used DMPfold to predict distance distributions for every residue pair in a set of proteins that showed both rigid and flexible behaviour. Residue pairs that were in contact in at least one reference structure were classified as rigid, flexible or neither. The predicted distance distribution of each residue pair was analysed for local maxima of probability indicating the most likely distance or distances between a pair of residues. We found that rigid residue pairs tended to have only a single local maximum in their predicted distance distributions while flexible residue pairs more often had multiple local maxima. These results suggest that the shape of predicted distance distributions contains information on the rigidity or flexibility of a protein and its constituent residues. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Dominik Schwarz
- Department of Statistics, University of Oxford, Oxford OX1 3LB, UK
| | - Guy Georges
- Department of Computational Engineering and Data Science, Large Molecule Research, Penzberg 82377, Germany
| | - Sebastian Kelm
- Computer-Aided Drug Design, UCB Pharma, Slough SL1 3WE, UK
| | - Jiye Shi
- Computer-Aided Drug Design, UCB Pharma, Slough SL1 3WE, UK
| | - Anna Vangone
- Department of Computational Engineering and Data Science, Large Molecule Research, Penzberg 82377, Germany
| | | |
Collapse
|
35
|
Abstract
In this paper, we offer an overview of a number of results on the static rigidity and infinitesimal rigidity of discrete structures which are embedded in projective geometric reasoning, representations, and transformations. Part I considers the fundamental case of a bar–joint framework in projective d-space and places particular emphasis on the projective invariance of infinitesimal rigidity, coning between dimensions, transfer to the spherical metric, slide joints and pure conditions for singular configurations. Part II extends the results, tools and concepts from Part I to additional types of rigid structures including body-bar, body–hinge and rod-bar frameworks, all drawing on projective representations, transformations and insights. Part III widens the lens to include the closely related cofactor matroids arising from multivariate splines, which also exhibit the projective invariance. These are another fundamental example of abstract rigidity matroids with deep analogies to rigidity. We conclude in Part IV with commentary on some nearby areas.
Collapse
|
36
|
Noby N, Auhim HS, Winter S, Worthy HL, Embaby AM, Saeed H, Hussein A, Pudney CR, Rizkallah PJ, Wells SA, Jones DD. Structure and in silico simulations of a cold-active esterase reveals its prime cold-adaptation mechanism. Open Biol 2021; 11:210182. [PMID: 34847772 PMCID: PMC8633780 DOI: 10.1098/rsob.210182] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Here we determined the structure of a cold active family IV esterase (EstN7) cloned from Bacillus cohnii strain N1. EstN7 is a dimer with a classical α/β hydrolase fold. It has an acidic surface that is thought to play a role in cold-adaption by retaining solvation under changed water solvent entropy at lower temperatures. The conformation of the functionally important cap region is significantly different to EstN7's closest relatives, forming a bridge-like structure with reduced helical content providing greater access to the active site through more than one substrate access tunnel. However, dynamics do not appear to play a major role in cold adaption. Molecular dynamics at different temperatures, rigidity analysis, normal mode analysis and geometric simulations of motion confirm the flexibility of the cap region but suggest that the rest of the protein is largely rigid. Rigidity analysis indicates the distribution of hydrophobic tethers is appropriate to colder conditions, where the hydrophobic effect is weaker than in mesophilic conditions due to reduced water entropy. Thus, it is likely that increased substrate accessibility and tolerance to changes in water entropy are important for of EstN7's cold adaptation rather than changes in dynamics.
Collapse
Affiliation(s)
- Nehad Noby
- Department of Biotechnology, Institute of Graduate Studies and Research, Alexandria University, Alexandria, Egypt,School of Biosciences, Molecular Biosciences Division, Cardiff University, Cardiff CF10 3AX, UK
| | - Husam Sabah Auhim
- School of Biosciences, Molecular Biosciences Division, Cardiff University, Cardiff CF10 3AX, UK,Department of Biology, College of Science, University of Baghdad, Baghdad, Iraq
| | - Samuel Winter
- Department of Biology and Biochemistry, University of Bath, Bath BA2 7AY, UK
| | - Harley L. Worthy
- School of Biosciences, Molecular Biosciences Division, Cardiff University, Cardiff CF10 3AX, UK
| | - Amira M. Embaby
- Department of Biotechnology, Institute of Graduate Studies and Research, Alexandria University, Alexandria, Egypt
| | - Hesham Saeed
- Department of Biotechnology, Institute of Graduate Studies and Research, Alexandria University, Alexandria, Egypt
| | - Ahmed Hussein
- Department of Biotechnology, Institute of Graduate Studies and Research, Alexandria University, Alexandria, Egypt
| | | | | | | | - D. Dafydd Jones
- School of Biosciences, Molecular Biosciences Division, Cardiff University, Cardiff CF10 3AX, UK
| |
Collapse
|
37
|
Hauptstein N, Pouyan P, Kehrein J, Dirauf M, Driessen MD, Raschig M, Licha K, Gottschaldt M, Schubert US, Haag R, Meinel L, Sotriffer C, Lühmann T. Molecular Insights into Site-Specific Interferon-α2a Bioconjugates Originated from PEG, LPG, and PEtOx. Biomacromolecules 2021; 22:4521-4534. [PMID: 34643378 DOI: 10.1021/acs.biomac.1c00775] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Conjugation of biologics with polymers modulates their pharmacokinetics, with polyethylene glycol (PEG) as the gold standard. We compared alternative polymers and two types of cyclooctyne linkers (BCN/DBCO) for bioconjugation of interferon-α2a (IFN-α2a) using 10 kDa polymers including linear mPEG, poly(2-ethyl-2-oxazoline) (PEtOx), and linear polyglycerol (LPG). IFN-α2a was azide functionalized via amber codon expansion and bioorthogonally conjugated to all cyclooctyne linked polymers. Polymer conjugation did not impact IFN-α2a's secondary structure and only marginally reduced IFN-α2a's bioactivity. In comparison to PEtOx, the LPG polymer attached via the less rigid cyclooctyne linker BCN was found to stabilize IFN-α2a against thermal stress. These findings were further detailed by molecular modeling studies which showed a modulation of protein flexibility upon PEtOx conjugation and a reduced amount of protein native contacts as compared to PEG and LPG originated bioconjugates. Polymer interactions with IFN-α2a were further assessed via a limited proteolysis (LIP) assay, which resulted in comparable proteolytic cleavage patterns suggesting weak interactions with the protein's surface. In conclusion, both PEtOx and LPG bioconjugates resulted in a similar biological outcome and may become promising PEG alternatives for bioconjugation.
Collapse
Affiliation(s)
- Niklas Hauptstein
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Paria Pouyan
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
| | - Josef Kehrein
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Michael Dirauf
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstr. 10, 07743 Jena, Germany.,Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
| | - Marc D Driessen
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Martina Raschig
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Kai Licha
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
| | - Michael Gottschaldt
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstr. 10, 07743 Jena, Germany.,Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
| | - Ulrich S Schubert
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstr. 10, 07743 Jena, Germany.,Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
| | - Rainer Haag
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
| | - Lorenz Meinel
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany.,Helmholtz Institute for RNA-Based Infection Research (HIRI), Helmholtz Center for Infection Research (HZI), 97080 Würzburg, Germany
| | - Christoph Sotriffer
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| | - Tessa Lühmann
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany
| |
Collapse
|
38
|
Perez-Mellor AF, Spezia R. Determination of kinetic properties in unimolecular dissociation of complex systems from graph theory based analysis of an ensemble of reactive trajectories. J Chem Phys 2021; 155:124103. [PMID: 34598552 DOI: 10.1063/5.0058382] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
In this paper, we report how graph theory can be used to analyze an ensemble of independent molecular trajectories, which can react during the simulation time-length, and obtain structural and kinetic information. This method is totally general and here is applied to the prototypical case of gas phase fragmentation of protonated cyclo-di-glycine. This methodology allows us to analyze the whole set of trajectories in an automatic computer-based way without the need of visual inspection but by getting all the needed information. In particular, we not only determine the appearance of different products and intermediates but also characterize the corresponding kinetics. The use of colored graph and canonical labeling allows for the correct characterization of the chemical species involved. In the present case, the simulations consist of an ensemble of unimolecular fragmentation trajectories at constant energy such that from the rate constants at different energies, the threshold energy can also be obtained for both global and specific pathways. This approach allows for the characterization of ion-molecule complexes, likely through a roaming mechanism, by properly taking into account the elusive nature of such species. Finally, it is possible to directly obtain the theoretical mass spectrum of the fragmenting species if the reacting system is an ion as in the specific example.
Collapse
Affiliation(s)
- Ariel F Perez-Mellor
- LAMBE UMR8587, Université d'Evry Val d'Essonne, CNRS, CEA, Université Paris-Saclay, Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, 91025 Evry, France
| | - Riccardo Spezia
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| |
Collapse
|
39
|
Koehl P, Orland H, Delarue M. Parameterizing elastic network models to capture the dynamics of proteins. J Comput Chem 2021; 42:1643-1661. [PMID: 34117647 DOI: 10.1002/jcc.26701] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Revised: 12/14/2020] [Accepted: 05/23/2021] [Indexed: 11/09/2022]
Abstract
Coarse-grained normal mode analyses of protein dynamics rely on the idea that the geometry of a protein structure contains enough information for computing its fluctuations around its equilibrium conformation. This geometry is captured in the form of an elastic network (EN), namely a network of edges between its residues. The normal modes of a protein are then identified with the normal modes of its EN. Different approaches have been proposed to construct ENs, focusing on the choice of the edges that they are comprised of, and on their parameterizations by the force constants associated with those edges. Here we propose new tools to guide choices on these two facets of EN. We study first different geometric models for ENs. We compare cutoff-based ENs, whose edges have lengths that are smaller than a cutoff distance, with Delaunay-based ENs and find that the latter provide better representations of the geometry of protein structures. We then derive an analytical method for the parameterization of the EN such that its dynamics leads to atomic fluctuations that agree with experimental B-factors. To limit overfitting, we attach a parameter referred to as flexibility constant to each atom instead of to each edge in the EN. The parameterization is expressed as a non-linear optimization problem whose parameters describe both rigid-body and internal motions. We show that this parameterization leads to improved ENs, whose dynamics mimic MD simulations better than ENs with uniform force constants, and reduces the number of normal modes needed to reproduce functional conformational changes.
Collapse
Affiliation(s)
- Patrice Koehl
- Department of Computer Sciences and Genome Center, University of California, Davis, California, USA
| | - Henri Orland
- Institut de Physique Théorique, Université Paris-Saclay, Gif sur Yvette, France
| | - Marc Delarue
- Unité de Dynamique Structurale des Macromolécules, Institut Pasteur, UMR 3528 du CNRS, Paris, France
| |
Collapse
|
40
|
Wei H, Wang B, Yang J, Gao J. RNA Flexibility Prediction With Sequence Profile and Predicted Solvent Accessibility. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2021; 18:2017-2022. [PMID: 31794403 DOI: 10.1109/tcbb.2019.2956496] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Structural flexibility plays an essential role in many biological processes. B-factor is an important indicator to measure the flexibility of protein or RNA structures. Many methods were developed to predict protein B-factors, but few studies have been done for RNA B-factor prediction. In this paper, we proposed a new method RNAbval to predict RNA B-factors using random forest. The method was developed using a comprehensive set of features, including the sequence profile and predicted solvent accessibility. RNAbval achieved an improvement of 9.2-20.5 percent over the state-of-the-art method on two benchmark test datasets. The proposed method is available at http://yanglab.nankai.edu.cn/RNAbval/.
Collapse
|
41
|
Becker D, Bharatam PV, Gohlke H. F/G Region Rigidity is Inversely Correlated to Substrate Promiscuity of Human CYP Isoforms Involved in Metabolism. J Chem Inf Model 2021; 61:4023-4030. [PMID: 34370479 DOI: 10.1021/acs.jcim.1c00558] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Of 57 human cytochrome P450 (CYP) enzymes, 12 metabolize 90% of xenobiotics. To our knowledge, no study has addressed the relation between enzyme dynamics and substrate promiscuity for more than three CYPs. Here, we show by constraint dilution simulations with the Constraint Network Analysis for the 12 isoforms that structural rigidity of the F/G region is significantly inversely correlated to the enzymes' substrate promiscuity. This highlights the functional importance of structural dynamics of the substrate tunnel.
Collapse
Affiliation(s)
- Daniel Becker
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Prasad V Bharatam
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research (NIPER), Sector 67, Sahibzada Ajit Singh Nagar, Mohali 160062, Punjab, India
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany.,John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| |
Collapse
|
42
|
Fowler NJ, Sljoka A, Williamson MP. The accuracy of NMR protein structures in the Protein Data Bank. Structure 2021; 29:1430-1439.e2. [PMID: 34331857 DOI: 10.1016/j.str.2021.07.001] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Revised: 06/18/2021] [Accepted: 07/14/2021] [Indexed: 11/18/2022]
Abstract
The program ANSURR measures the accuracy of NMR structures by comparing rigidity obtained from experimental backbone chemical shifts and from structures. We report on ANSURR analysis of 7,000 PDB NMR ensembles within the Protein Data Bank, which can be found at ansurr.com. The accuracy of NMR structures progressively improved up until 2005, but since then, it has plateaued. Most structures have accurate secondary structure, but are generally too floppy, particularly in loops. Thus, there is a need for more experimental restraints in loops. Currently, the best predictors of accuracy are Ramachandran distribution and the number of NOE restraints per residue. The precision of structures within the ensemble correlates well with accuracy, as does the number of hydrogen bond restraints per residue. Structure accuracy is improved when other components (such as additional polypeptide chains or ligands) are included.
Collapse
Affiliation(s)
- Nicholas J Fowler
- Department of Molecular Biology and Biotechnology, University of Sheffield, Sheffield S10 2TN, UK
| | - Adnan Sljoka
- RIKEN Center for Advanced Intelligence Project, RIKEN, 1-4-1 Nihombashi, Chuo-ku, Tokyo 103-0027 Japan; Department of Chemistry, University of Toronto, UTM, 3359 Mississauga Road North, Mississauga, ON L5L 1C6, Canada.
| | - Mike P Williamson
- Department of Molecular Biology and Biotechnology, University of Sheffield, Sheffield S10 2TN, UK.
| |
Collapse
|
43
|
El Harrar T, Frieg B, Davari MD, Jaeger KE, Schwaneberg U, Gohlke H. Aqueous ionic liquids redistribute local enzyme stability via long-range perturbation pathways. Comput Struct Biotechnol J 2021; 19:4248-4264. [PMID: 34429845 PMCID: PMC8355836 DOI: 10.1016/j.csbj.2021.07.001] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Revised: 07/05/2021] [Accepted: 07/06/2021] [Indexed: 01/25/2023] Open
Abstract
Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for biocatalysis due to their unique properties. On the other hand, the incubation of enzymes in IL or aIL often reduces enzyme activity. Recent studies proposed various aIL-induced effects to explain the reduction, classified as direct effects, e.g., local dehydration or competitive inhibition, and indirect effects, e.g., structural perturbations or disturbed catalytic site integrity. However, the molecular origin of indirect effects has largely remained elusive. Here we show by multi-μs long molecular dynamics simulations, free energy computations, and rigidity analyses that aIL favorably interact with specific residues of Bacillus subtilis Lipase A (BsLipA) and modify the local structural stability of this model enzyme by inducing long-range perturbations of noncovalent interactions. The perturbations percolate over neighboring residues and eventually affect the catalytic site and the buried protein core. Validation against a complete experimental site saturation mutagenesis library of BsLipA (3620 variants) reveals that the residues of the perturbation pathways are distinguished sequence positions where substitutions highly likely yield significantly improved residual activity. Our results demonstrate that identifying these perturbation pathways and specific IL ion-residue interactions there effectively predicts focused variant libraries with improved aIL tolerance.
Collapse
Affiliation(s)
- Till El Harrar
- Institute of Biotechnology, RWTH Aachen University, 52074 Aachen, Germany
- John-von-Neumann-Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Benedikt Frieg
- John-von-Neumann-Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Mehdi D. Davari
- Institute of Biotechnology, RWTH Aachen University, 52074 Aachen, Germany
| | - Karl-Erich Jaeger
- Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, 52428 Jülich, Germany
- Institute of Bio- and Geosciences IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Ulrich Schwaneberg
- Institute of Biotechnology, RWTH Aachen University, 52074 Aachen, Germany
- DWI – Leibniz Institute for Interactive Materials e.V., 52074 Aachen, Germany
| | - Holger Gohlke
- John-von-Neumann-Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
44
|
Lata S, Akif M. Comparative protein structure network analysis on 3CL pro from SARS-CoV-1 and SARS-CoV-2. Proteins 2021; 89:1216-1225. [PMID: 33983654 PMCID: PMC8242809 DOI: 10.1002/prot.26143] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2020] [Revised: 04/10/2021] [Accepted: 05/05/2021] [Indexed: 12/29/2022]
Abstract
The main protease Mpro, 3CLpro is an important target from coronaviruses. In spite of having 96% sequence identity among Mpros from SARS‐CoV‐1 and SARS‐CoV‐2; the inhibitors used to block the activity of SARS‐CoV‐1 Mpro so far, were found to have differential inhibitory effect on Mpro of SARS‐CoV‐2. The possible reason could be due to the difference of few amino acids among the peptidases. Since, overall 3‐D crystallographic structure of Mpro from SARS‐CoV‐1 and SARS‐CoV‐2 is quite similar and mapping a subtle structural variation is seemingly impossible. Hence, we have attempted to study a structural comparison of SARS‐CoV‐1 and SARS‐CoV‐2 Mpro in apo and inhibitor bound states using protein structure network (PSN) based approach at contacts level. The comparative PSNs analysis of apo Mpros from SARS‐CoV‐1 and SARS‐CoV‐2 uncovers small but significant local changes occurring near the active site region and distributed throughout the structure. Additionally, we have shown how inhibitor binding perturbs the PSG and the communication pathways in Mpros. Moreover, we have also investigated the network connectivity on the quaternary structure of Mpro and identified critical residue pairs for complex formation using three centrality measurement parameters along with the modularity analysis. Taken together, these results on the comparative PSN provide an insight into conformational changes that may be used as an additional guidance towards specific drug development.
Collapse
Affiliation(s)
- Surabhi Lata
- Department of Biochemistry, School of Life Sciences, University of Hyderabad, Hyderabad, India
| | - Mohd Akif
- Department of Biochemistry, School of Life Sciences, University of Hyderabad, Hyderabad, India
| |
Collapse
|
45
|
Nutschel C, Coscolín C, David B, Mulnaes D, Ferrer M, Jaeger KE, Gohlke H. Promiscuous Esterases Counterintuitively Are Less Flexible than Specific Ones. J Chem Inf Model 2021; 61:2383-2395. [PMID: 33949194 DOI: 10.1021/acs.jcim.1c00152] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Understanding mechanisms of promiscuity is increasingly important from a fundamental and application point of view. As to enzyme structural dynamics, more promiscuous enzymes generally have been recognized to also be more flexible. However, examples for the opposite received much less attention. Here, we exploit comprehensive experimental information on the substrate promiscuity of 147 esterases tested against 96 esters together with computationally efficient rigidity analyses to understand the molecular origin of the observed promiscuity range. Unexpectedly, our data reveal that promiscuous esterases are significantly less flexible than specific ones, are significantly more thermostable, and have a significantly increased specific activity. These results may be reconciled with a model according to which structural flexibility in the case of specific esterases serves for conformational proofreading. Our results signify that an esterase sequence space can be screened by rigidity analyses for promiscuous esterases as starting points for further exploration in biotechnology and synthetic chemistry.
Collapse
Affiliation(s)
- Christina Nutschel
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Cristina Coscolín
- Institute of Catalysis, Consejo Superior de Investigaciones Científicas, 28049 Madrid, Spain
| | - Benoit David
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Daniel Mulnaes
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Manuel Ferrer
- Institute of Catalysis, Consejo Superior de Investigaciones Científicas, 28049 Madrid, Spain
| | - Karl-Erich Jaeger
- Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, 52425 Jülich, Germany.,Institute of Bio- and Geosciences IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Holger Gohlke
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany.,Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
46
|
|
47
|
Huang SK, Pandey A, Tran DP, Villanueva NL, Kitao A, Sunahara RK, Sljoka A, Prosser RS. Delineating the conformational landscape of the adenosine A 2A receptor during G protein coupling. Cell 2021; 184:1884-1894.e14. [PMID: 33743210 DOI: 10.1016/j.cell.2021.02.041] [Citation(s) in RCA: 81] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Revised: 12/02/2020] [Accepted: 02/17/2021] [Indexed: 12/17/2022]
Abstract
G-protein-coupled receptors (GPCRs) represent a ubiquitous membrane protein family and are important drug targets. Their diverse signaling pathways are driven by complex pharmacology arising from a conformational ensemble rarely captured by structural methods. Here, fluorine nuclear magnetic resonance spectroscopy (19F NMR) is used to delineate key functional states of the adenosine A2A receptor (A2AR) complexed with heterotrimeric G protein (Gαsβ1γ2) in a phospholipid membrane milieu. Analysis of A2AR spectra as a function of ligand, G protein, and nucleotide identifies an ensemble represented by inactive states, a G-protein-bound activation intermediate, and distinct nucleotide-free states associated with either partial- or full-agonist-driven activation. The Gβγ subunit is found to be critical in facilitating ligand-dependent allosteric transmission, as shown by 19F NMR, biochemical, and computational studies. The results provide a mechanistic basis for understanding basal signaling, efficacy, precoupling, and allostery in GPCRs.
Collapse
Affiliation(s)
- Shuya Kate Huang
- Department of Chemistry, University of Toronto, UTM, 3359 Mississauga Road North, Mississauga, Ontario L5L 1C6, Canada
| | - Aditya Pandey
- Department of Chemistry, University of Toronto, UTM, 3359 Mississauga Road North, Mississauga, Ontario L5L 1C6, Canada; Department of Biochemistry, University of Toronto, 1 King's College Circle, Toronto, Ontario M5S 1A8, Canada
| | - Duy Phuoc Tran
- School of Life Science and Technology, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8550, Japan
| | - Nicolas L Villanueva
- Department of Pharmacology, University of California San Diego School of Medicine, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Akio Kitao
- School of Life Science and Technology, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8550, Japan
| | - Roger K Sunahara
- Department of Pharmacology, University of California San Diego School of Medicine, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Adnan Sljoka
- Department of Chemistry, University of Toronto, UTM, 3359 Mississauga Road North, Mississauga, Ontario L5L 1C6, Canada; RIKEN Center for Advanced Intelligence Project, RIKEN, 1-4-1 Nihombashi, Chuo-ku, Tokyo 103-0027, Japan.
| | - R Scott Prosser
- Department of Chemistry, University of Toronto, UTM, 3359 Mississauga Road North, Mississauga, Ontario L5L 1C6, Canada; Department of Biochemistry, University of Toronto, 1 King's College Circle, Toronto, Ontario M5S 1A8, Canada.
| |
Collapse
|
48
|
Monoclonal antibody stability can be usefully monitored using the excitation-energy-dependent fluorescence edge-shift. Biochem J 2021; 477:3599-3612. [PMID: 32869839 PMCID: PMC7527260 DOI: 10.1042/bcj20200580] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2020] [Revised: 08/27/2020] [Accepted: 09/01/2020] [Indexed: 11/17/2022]
Abstract
Among the major challenges in the development of biopharmaceuticals are structural heterogeneity and aggregation. The development of a successful therapeutic monoclonal antibody (mAb) requires both a highly active and also stable molecule. Whilst a range of experimental (biophysical) approaches exist to track changes in stability of proteins, routine prediction of stability remains challenging. The fluorescence red edge excitation shift (REES) phenomenon is sensitive to a range of changes in protein structure. Based on recent work, we have found that quantifying the REES effect is extremely sensitive to changes in protein conformational state and dynamics. Given the extreme sensitivity, potentially this tool could provide a ‘fingerprint’ of the structure and stability of a protein. Such a tool would be useful in the discovery and development of biopharamceuticals and so we have explored our hypothesis with a panel of therapeutic mAbs. We demonstrate that the quantified REES data show remarkable sensitivity, being able to discern between structurally identical antibodies and showing sensitivity to unfolding and aggregation. The approach works across a broad concentration range (µg–mg/ml) and is highly consistent. We show that the approach can be applied alongside traditional characterisation testing within the context of a forced degradation study (FDS). Most importantly, we demonstrate the approach is able to predict the stability of mAbs both in the short (hours), medium (days) and long-term (months). The quantified REES data will find immediate use in the biopharmaceutical industry in quality assurance, formulation and development. The approach benefits from low technical complexity, is rapid and uses instrumentation which exists in most biochemistry laboratories without modification.
Collapse
|
49
|
Pfleger C, Kusch J, Kondapuram M, Schwabe T, Sattler C, Benndorf K, Gohlke H. Allosteric signaling in C-linker and cyclic nucleotide-binding domain of HCN2 channels. Biophys J 2021; 120:950-963. [PMID: 33515603 DOI: 10.1016/j.bpj.2021.01.017] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Revised: 01/04/2021] [Accepted: 01/19/2021] [Indexed: 12/22/2022] Open
Abstract
Opening of hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels is controlled by membrane hyperpolarization and binding of cyclic nucleotides to the tetrameric cyclic nucleotide-binding domain (CNBD), attached to the C-linker (CL) disk. Confocal patch-clamp fluorometry revealed pronounced cooperativity of ligand binding among protomers. However, by which pathways allosteric signal transmission occurs remained elusive. Here, we investigate how changes in the structural dynamics of the CL-CNBD of mouse HCN2 upon cAMP binding relate to inter- and intrasubunit signal transmission. Applying a rigidity-theory-based approach, we identify two intersubunit and one intrasubunit pathways that differ in allosteric coupling strength between cAMP-binding sites or toward the CL. These predictions agree with results from electrophysiological and patch-clamp fluorometry experiments. Our results map out distinct routes within the CL-CNBD that modulate different cAMP-binding responses in HCN2 channels. They signify that functionally relevant submodules may exist within and across structurally discernable subunits in HCN channels.
Collapse
Affiliation(s)
- Christopher Pfleger
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Jana Kusch
- Institute of Physiology II, Jena University Hospital, Jena, Germany
| | | | - Tina Schwabe
- Institute of Physiology II, Jena University Hospital, Jena, Germany
| | | | - Klaus Benndorf
- Institute of Physiology II, Jena University Hospital, Jena, Germany
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany; John von Neumann Institute for Computing, Jülich Supercomputing Centre, and Institute of Biological Information Processing (IBI-7, Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany.
| |
Collapse
|
50
|
Sljoka A. Probing Allosteric Mechanism with Long-Range Rigidity Transmission Across Protein Networks. Methods Mol Biol 2021; 2253:61-75. [PMID: 33315218 DOI: 10.1007/978-1-0716-1154-8_5] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Allosteric transmission refers to regulation of protein function at a distance. "Allostery" involves regulation and/or signal transduction induced by a perturbation event. Allostery, which has been coined the "second secret of life," is a fundamental property of most dynamics proteins. Most of critical questions surrounding allostery are largely unresolved. One of the key puzzles is to describe the physical mechanism of distant coupled conformational change. Another hot research area surrounding allostery is detection of allosteric pathways or regions (residues) in the protein that are the most critical for transmission of allosteric information. Using techniques inspired by mathematical rigidity theory and mechanical linkages, we have previously proposed a mechanistic model and description of allosteric transmission and an accompanying computational method, the Rigidity Transmission Allostery (RTA) algorithm. The RTA algorithm and method are designed to predict if mechanical perturbation of rigidity, for example, due to ligand binding, at one site of the protein can transmit and propagate across a protein structure and in turn cause a change in available conformational degrees of freedom and a change in conformation at a second distant site, equivalently resulting in allosteric transmission. The RTA algorithm is computationally very fast and can rapidly scan many unknown sites for allosteric transmission, identifying potential novel allosteric sites and quantify their allosteric effect. In this chapter we will discuss the rigidity-based mechanistic model of allosteric communication. As a case illustrative study, we will demonstrate RTA analysis on a G protein coupled receptor (GPCR) human adenosine A2A receptor. Our method gives important implications and a novel prospective for general mechanistic description of allosteric communication.
Collapse
Affiliation(s)
- Adnan Sljoka
- Department of Informatics, School of Science and Technology, Kwansei Gakuin University, Sanda, Hyogo, Japan.
- CREST, Japan Science and Technology Agency (JST), Tokyo, Japan.
- RIKEN, Center for Advanced Intelligence Project, Tokyo, Japan.
| |
Collapse
|