1
|
Parafioriti M, Elli S, Muñoz-García JC, Ramírez-Cárdenas J, Yates EA, Angulo J, Guerrini M. Differential Solvent DEEP-STD NMR and MD Simulations Enable the Determinants of the Molecular Recognition of Heparin Oligosaccharides by Antithrombin to Be Disentangled. Int J Mol Sci 2024; 25:4669. [PMID: 38731888 PMCID: PMC11083112 DOI: 10.3390/ijms25094669] [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: 02/28/2024] [Revised: 04/16/2024] [Accepted: 04/20/2024] [Indexed: 05/13/2024] Open
Abstract
The interaction of heparin with antithrombin (AT) involves a specific sequence corresponding to the pentasaccharide GlcNAc/NS6S-GlcA-GlcNS3S6S-IdoA2S-GlcNS6S (AGA*IA). Recent studies have revealed that two AGA*IA-containing hexasaccharides, which differ in the sulfation degree of the iduronic acid unit, exhibit similar binding to AT, albeit with different affinities. However, the lack of experimental data concerning the molecular contacts between these ligands and the amino acids within the protein-binding site prevents a detailed description of the complexes. Differential epitope mapping (DEEP)-STD NMR, in combination with MD simulations, enables the experimental observation and comparison of two heparin pentasaccharides interacting with AT, revealing slightly different bound orientations and distinct affinities of both glycans for AT. We demonstrate the effectiveness of the differential solvent DEEP-STD NMR approach in determining the presence of polar residues in the recognition sites of glycosaminoglycan-binding proteins.
Collapse
Affiliation(s)
- Michela Parafioriti
- Istituto di Ricerche Chimiche e Biochimiche “G. Ronzoni”, Via Giuseppe Colombo 81, 20133 Milano, Italy; (M.P.); (S.E.)
| | - Stefano Elli
- Istituto di Ricerche Chimiche e Biochimiche “G. Ronzoni”, Via Giuseppe Colombo 81, 20133 Milano, Italy; (M.P.); (S.E.)
| | - Juan C. Muñoz-García
- Instituto de Investigationes Químicas (IIQ)-Consejo Superior de Investigaciones Científicas (CSIC), Avenida Americo Vespucio 49, 41092 Sevilla, Spain; (J.C.M.-G.); (J.R.-C.)
| | - Jonathan Ramírez-Cárdenas
- Instituto de Investigationes Químicas (IIQ)-Consejo Superior de Investigaciones Científicas (CSIC), Avenida Americo Vespucio 49, 41092 Sevilla, Spain; (J.C.M.-G.); (J.R.-C.)
| | - Edwin A. Yates
- Department of Biochemistry and Systems Biology, Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool L69 7ZB, UK;
- Centre for Glycoscience, Keele University, Newcastle-Under-Lyme ST5 5BG, UK
| | - Jesús Angulo
- Instituto de Investigationes Químicas (IIQ)-Consejo Superior de Investigaciones Científicas (CSIC), Avenida Americo Vespucio 49, 41092 Sevilla, Spain; (J.C.M.-G.); (J.R.-C.)
| | - Marco Guerrini
- Istituto di Ricerche Chimiche e Biochimiche “G. Ronzoni”, Via Giuseppe Colombo 81, 20133 Milano, Italy; (M.P.); (S.E.)
| |
Collapse
|
2
|
Lutsyk V, Plazinski W. Exploring Ring Conformation in Uronate Monosaccharides: Insights from Ab Initio Calculations and Classical Molecular Dynamics Simulations. J Phys Chem B 2024; 128:472-491. [PMID: 38170925 DOI: 10.1021/acs.jpcb.3c06556] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
The study focuses on the conformational properties of biologically relevant monosaccharides belonging to the group of uronates: α-l-iduronate, O2-sulfated-α-l-iduronate, and O2-sulfated-α-l-guluronate, either unfunctionalized or O1-methylated. We applied the previously proposed two-step methodology, combining classical MD simulations and subsequent ab initio (QM) calculations, performed on a rationally subsampled set of molecular configurations. We found that, regardless of the number of molecular configurations considered, the level of theory, and the weighting scheme applied, none of the QM approaches is capable of predicting the correct conformational equilibrium of sulfated iduronates as long as the tight counterion binding is not considered. Multicenter, ring-shape-specific binding of either Na+ or Ca2+ ions drastically shifts the conformational equilibrium of the pyranose ring in sulfated iduronates toward 1C4 but does not significantly affect the conformation of non-sulfated compounds. A similar shift is observed upon the protonation of carboxyl groups in all iduronates. In addition, we report a set of average J-coupling constant values related to vicinal protons bound to the pyranose ring of iduronates and corresponding to each of the three main groups of ring conformers, i.e., 4C1, B/S (boat/skew boat), and 1C4. In combination with the conformational energies or with the experimental data, these values allowed the relative proportions of the ring conformers to be estimated and the Karplus-type equations linking the 3JHH-coupling constants to the torsion angles within the pyranose ring to be refined.
Collapse
Affiliation(s)
- Valery Lutsyk
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, 30-239 Krakow, Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, 30-239 Krakow, Poland
- Department of Biopharmacy, Medical University of Lublin, Chodzki 4a, 20-093 Lublin, Poland
| |
Collapse
|
3
|
Marcisz M, Anila S, Gaardløs M, Zacharias M, Samsonov SA. Studying specificity in protein-glycosaminoglycan recognition with umbrella sampling. Beilstein J Org Chem 2023; 19:1933-1946. [PMID: 38170083 PMCID: PMC10760462 DOI: 10.3762/bjoc.19.144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Accepted: 12/07/2023] [Indexed: 01/05/2024] Open
Abstract
In the past few decades, glycosaminoglycan (GAG) research has been crucial for gaining insights into various physiological, pathological, and therapeutic aspects mediated by the direct interactions between the GAG molecules and diverse proteins. The structural and functional heterogeneities of GAGs as well as their ability to bind specific proteins are determined by the sugar composition of the GAG, the size of the GAG chains, and the degree and pattern of sulfation. A deep understanding of the interactions in protein-GAG complexes is essential to explain their biological functions. In this study, the umbrella sampling (US) approach is used to pull away a GAG ligand from the binding site and then pull it back in. We analyze the binding interactions between GAGs of three types (heparin, desulfated heparan sulfate, and chondroitin sulfate) with three different proteins (basic fibroblast growth factor, acidic fibroblast growth factor, and cathepsin K). The main focus of our study was to evaluate whether the US approach is able to reproduce experimentally obtained structures, and how useful it can be for getting a deeper understanding of GAG properties, especially protein recognition specificity and multipose binding. We found that the binding free energy landscape in the proximity of the GAG native binding pose is complex and implies the co-existence of several binding poses. The sliding of a GAG chain along a protein surface could be a potential mechanism of GAG particular sequence recognition by proteins.
Collapse
Affiliation(s)
- Mateusz Marcisz
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
- Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Gdańsk, Poland
| | | | | | - Martin Zacharias
- Physics Department, Technical University of Munich, Garching, Germany
| | | |
Collapse
|
4
|
Schulze C, Danielsson A, Liwo A, Huster D, Samsonov SA, Penk A. Ligand binding of interleukin-8: a comparison of glycosaminoglycans and acidic peptides. Phys Chem Chem Phys 2023; 25:24930-24947. [PMID: 37694394 DOI: 10.1039/d3cp02457a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2023]
Abstract
Recognition and binding of regulatory proteins to glycosaminoglycans (GAGs) from the extracellular matrix is a process of high biological importance. The interaction between negatively charged sulfate or carboxyl groups of the GAGs and clusters of basic amino acids on the protein is crucial in this binding process and it is believed that electrostatics represent the key factor for this interaction. However, given the rather undirected nature of electrostatics, it is important to achieve a clear understanding of its role in protein-GAG interactions and how specificity and selectivity in these systems can be achieved, when the classical key-lock binding motif is not applicable. Here, we compare protein binding of a highly charged heparin (HP) hexasaccharide with four de novo designed decapeptides of varying negative net charge. The charge density of these peptides was comparable to typical GAGs of the extracellular matrix. We used the regulatory protein interleukin-8 (IL-8) because its interactions with GAGs are well described. All four peptide ligands bind to the same epitope of IL-8 but show much weaker binding affinity as revealed in 1H-15N HSQC NMR titration experiments. Complementary molecular docking and molecular dynamics simulations revealed further atomistic details of the interaction mode of GAG versus peptide ligands. Overall, similar contributions to the binding energy and hydrogen bond formation are determined for HP and the highly charged peptides, suggesting that the entropic loss of the peptides upon binding likely account for the remarkably different affinity of GAG versus peptide ligands to IL-8.
Collapse
Affiliation(s)
- Christian Schulze
- Institute for Medical Physics and Biophysics, University of Leipzig, Härtelstr. 16/18, 04107 Leipzig, Germany.
| | - Annemarie Danielsson
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland.
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland.
| | - Daniel Huster
- Institute for Medical Physics and Biophysics, University of Leipzig, Härtelstr. 16/18, 04107 Leipzig, Germany.
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland.
| | - Anja Penk
- Institute for Medical Physics and Biophysics, University of Leipzig, Härtelstr. 16/18, 04107 Leipzig, Germany.
| |
Collapse
|
5
|
Sankaranarayanan NV, Sistla S, Nagarajan B, Chittum JE, Lau JTY, Desai UR. Computational studies on glycosaminoglycan recognition of sialyl transferases. Glycobiology 2023; 33:579-590. [PMID: 37171590 PMCID: PMC10426320 DOI: 10.1093/glycob/cwad040] [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/29/2023] [Revised: 05/03/2023] [Accepted: 05/11/2023] [Indexed: 05/13/2023] Open
Abstract
Despite decades of research, glycosaminoglycans (GAGs) have not been known to interact with sialyl transferases (STs). Using our in-house combinatorial virtual library screening (CVLS) technology, we studied seven human isoforms, including ST6GAL1, ST6GAL2, ST3GAL1, ST3GAL3, ST3GAL4, ST3GAL5, and ST3GAL6, and predicted that GAGs, especially heparan sulfate (HS), are likely to differentially bind to STs. Exhaustive CVLS and molecular dynamics studies suggested that the common hexasaccharide sequence of HS preferentially recognized ST6GAL1 in a site overlapping the binding site of the donor substrate CMP-Sia. Interestingly, CVLS did not ascribe any special role for the rare 3-O-sulfate modification of HS in ST6GAL1 recognition. The computational predictions were tested using spectrofluorimetric studies, which confirmed preferential recognition of HS over other GAGs. A classic chain length-dependent binding of GAGs to ST6GAL1 was observed with polymeric HS displaying a tight affinity of ~65 nM. Biophysical studies also confirmed a direct competition between CMP-Sia and an HS oligosaccharide and CS polysaccharide for binding to ST6GAL1. Overall, our novel observation that GAGs bind to ST6GAL1 with high affinity and compete with the donor substrate is likely to be important because modulation of sialylation of glycan substrates on cells has considerable physiological/pathological consequences. Our work also brings forth the possibility of developing GAG-based chemical probes of ST6GAL1.
Collapse
Affiliation(s)
- Nehru Viji Sankaranarayanan
- Department of Medicinal Chemistry, Virginia Commonwealth University, Virginia 23298, Richmond, USA
- Drug Discovery and Development, Virginia Commonwealth University, Institute for Structural Biology, Virginia 23219, Richmond, USA
| | - Srinivas Sistla
- Department of Medicinal Chemistry, Virginia Commonwealth University, Virginia 23298, Richmond, USA
- Drug Discovery and Development, Virginia Commonwealth University, Institute for Structural Biology, Virginia 23219, Richmond, USA
| | - Balaji Nagarajan
- Department of Medicinal Chemistry, Virginia Commonwealth University, Virginia 23298, Richmond, USA
- Drug Discovery and Development, Virginia Commonwealth University, Institute for Structural Biology, Virginia 23219, Richmond, USA
| | - John E Chittum
- Department of Medicinal Chemistry, Virginia Commonwealth University, Virginia 23298, Richmond, USA
- Drug Discovery and Development, Virginia Commonwealth University, Institute for Structural Biology, Virginia 23219, Richmond, USA
| | - Joseph T Y Lau
- Department of Molecular and Cellular Biology, Roswell Park Comprehensive Cancer Center, Buffalo 14263, NY, USA
| | - Umesh R Desai
- Department of Medicinal Chemistry, Virginia Commonwealth University, Virginia 23298, Richmond, USA
- Drug Discovery and Development, Virginia Commonwealth University, Institute for Structural Biology, Virginia 23219, Richmond, USA
| |
Collapse
|
6
|
Bojarski KK, Samsonov SA. In silico insights into procathepsin S maturation mediated by glycosaminoglycans. J Mol Graph Model 2023; 120:108406. [PMID: 36707295 DOI: 10.1016/j.jmgm.2023.108406] [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: 10/10/2022] [Revised: 11/24/2022] [Accepted: 01/10/2023] [Indexed: 01/15/2023]
Abstract
Procathepsins, inactive precursors of cathepsins are present in the extracellular matrix (ECM) and in lysosomes. Their active forms are involved in a number of biologically relevant processes, including bone resorption, intracellular proteolysis and regulation of programmed cell death. These processes might be mediated by glycosaminoglycans (GAGs), long unbranched periodic negatively charged polysaccharides. GAGs are also present in ECM and play important role in anticoagulation, angiogenesis and tissue regeneration. GAGs not only mediate the enzymatic activity of cathepsins but can also regulate the process of procathepsin maturation, as it was shown for procathepsin B and S. In this study, we propose the molecular mechanism underlying the biological role of GAGs in procathepsin S maturation and compare our findings with computational data obtained for procathepsin B. We rigorously analyse procathepsin S-GAG complexes in terms of their dynamics, free energy and potential allosteric regulation. We conclude that the GAG binding region might have an effect on the dynamics of procathepsin S structure and so affect its maturation by two different mechanisms.
Collapse
Affiliation(s)
- Krzysztof K Bojarski
- Department of Physical Chemistry, Gdansk University of Technology, Narutowicza 11/12, Gdansk, 80-233, Poland.
| | - Sergey A Samsonov
- Department of Theoretical Chemistry, University of Gdansk, Wita Stwosza 63, Gdansk, 80-308, Poland
| |
Collapse
|
7
|
Marcisz M, Samsonov SA. Solvent Model Benchmark for Molecular Dynamics of Glycosaminoglycans. J Chem Inf Model 2023; 63:2147-2157. [PMID: 36989082 PMCID: PMC10091405 DOI: 10.1021/acs.jcim.2c01472] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
In computational studies of glycosaminoglycans (GAGs), a group of anionic, periodic linear polysaccharides, so far there has been very little discussion about the role of solvent models in the molecular dynamics simulations of these molecules. Predominantly, the TIP3P water model is commonly used as one of the most popular explicit water models in general. However, there are numerous alternative explicit and implicit water models that are neglected in the computational research of GAGs. Since solvent-mediated interactions are particularly important for GAG dynamic and structural properties, it would be of great interest for the GAG community to establish the solvent model that is suited the best in terms of the quality of theoretically obtained GAG parameters and, at the same time, would be reasonably demanding in terms of computational resources required. In this study, heparin (HP) was simulated using five implicit and six explicit solvent models with the aim to find out how different solvent models influence HP's molecular descriptors in the molecular dynamics simulations. Here, we initiate the search for the most appropriate solvent representation for GAG systems and we hope to encourage other groups to contribute to this highly relevant subject.
Collapse
Affiliation(s)
- Mateusz Marcisz
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
- Intercollegiate Faculty of Biotechnology of UG and MUG, ul. Abrahama 58, 80-307 Gdańsk, Poland
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
8
|
Hogwood J, Mulloy B, Lever R, Gray E, Page CP. Pharmacology of Heparin and Related Drugs: An Update. Pharmacol Rev 2023; 75:328-379. [PMID: 36792365 DOI: 10.1124/pharmrev.122.000684] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2022] [Revised: 11/04/2022] [Accepted: 11/08/2022] [Indexed: 02/17/2023] Open
Abstract
Heparin has been used extensively as an antithrombotic and anticoagulant for close to 100 years. This anticoagulant activity is attributed mainly to the pentasaccharide sequence, which potentiates the inhibitory action of antithrombin, a major inhibitor of the coagulation cascade. More recently it has been elucidated that heparin exhibits anti-inflammatory effect via interference of the formation of neutrophil extracellular traps and this may also contribute to heparin's antithrombotic activity. This illustrates that heparin interacts with a broad range of biomolecules, exerting both anticoagulant and nonanticoagulant actions. Since our previous review, there has been an increased interest in these nonanticoagulant effects of heparin, with the beneficial role in patients infected with SARS2-coronavirus a highly topical example. This article provides an update on our previous review with more recent developments and observations made for these novel uses of heparin and an overview of the development status of heparin-based drugs. SIGNIFICANCE STATEMENT: This state-of-the-art review covers recent developments in the use of heparin and heparin-like materials as anticoagulant, now including immunothrombosis observations, and as nonanticoagulant including a role in the treatment of SARS-coronavirus and inflammatory conditions.
Collapse
Affiliation(s)
- John Hogwood
- Sackler Institute of Pulmonary Pharmacology, Institute of Pharmaceutical Science, King's College London, London, United Kingdom (B.M., E.G., C.P.P.); National Institute for Biological Standards and Control, South Mimms, Hertfordshire, United Kingdom (J.H., E.G.) and School of Pharmacy, University College London, London, United Kingdom (R.L.)
| | - Barbara Mulloy
- Sackler Institute of Pulmonary Pharmacology, Institute of Pharmaceutical Science, King's College London, London, United Kingdom (B.M., E.G., C.P.P.); National Institute for Biological Standards and Control, South Mimms, Hertfordshire, United Kingdom (J.H., E.G.) and School of Pharmacy, University College London, London, United Kingdom (R.L.)
| | - Rebeca Lever
- Sackler Institute of Pulmonary Pharmacology, Institute of Pharmaceutical Science, King's College London, London, United Kingdom (B.M., E.G., C.P.P.); National Institute for Biological Standards and Control, South Mimms, Hertfordshire, United Kingdom (J.H., E.G.) and School of Pharmacy, University College London, London, United Kingdom (R.L.)
| | - Elaine Gray
- Sackler Institute of Pulmonary Pharmacology, Institute of Pharmaceutical Science, King's College London, London, United Kingdom (B.M., E.G., C.P.P.); National Institute for Biological Standards and Control, South Mimms, Hertfordshire, United Kingdom (J.H., E.G.) and School of Pharmacy, University College London, London, United Kingdom (R.L.)
| | - Clive P Page
- Sackler Institute of Pulmonary Pharmacology, Institute of Pharmaceutical Science, King's College London, London, United Kingdom (B.M., E.G., C.P.P.); National Institute for Biological Standards and Control, South Mimms, Hertfordshire, United Kingdom (J.H., E.G.) and School of Pharmacy, University College London, London, United Kingdom (R.L.)
| |
Collapse
|
9
|
Dhurua S, Jana M. Understanding Conformational Properties and Role of Hydrogen Bonds in Glycosaminoglycans-Interleukin8 Complexes in Aqueous Medium by Molecular Dynamics Simulation. Chemphyschem 2023; 24:e202200440. [PMID: 36239153 DOI: 10.1002/cphc.202200440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Revised: 10/13/2022] [Indexed: 11/11/2022]
Abstract
Atomistic molecular dynamics simulations were performed under ambient conditions to explore the conformational features and binding affinities of hexameric glycosaminoglycans (GAGs) with chemokine Interleukin8 (IL8) in an aqueous medium. We tried to understand the role of hydrogen bonds (HBs) involving conserved water in mediating the interactions. The Luzar-Chandler model was adopted to study the kinetics of HB breaking and formation concerning different water-mediated HBs. The conformational flexibilities of bound GAGs are due to the flexible glycosidic linkages than the occasional/rare ring pucker conformation. The free energy landscape constructed with ϕ, and ψ, depicted that different conformational minima associated with the glycosidic linkage flexibility of the GAGs in bound states are separated by energy barriers. The binding affinities of IL8 towards GAGs are favored through the electrostatic and non-polar solvation interactions. 4-different types of conserved water were explored in the solvent-mediated binding of GAGs with IL8. The average lifetime of the IL8-GAG direct HB pairs was ∼ten times less than the IL8-GAG-shared water HBs. This is due to the rapid establishment of HB breaking and reformation kinetics involving water of a shared layer. We find that despite the highly negatively charged surface of GAGs, the IL8 surface populated by non-cationic amino acids could serve as a promising binding site in addition to the cationic surface of the protein.
Collapse
Affiliation(s)
- Shakuntala Dhurua
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, 769008, Rourkela, India
| | - Madhurima Jana
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, 769008, Rourkela, India
| |
Collapse
|
10
|
Pągielska M, Samsonov SA. Molecular Dynamics-Based Comparative Analysis of Chondroitin and Dermatan Sulfates. Biomolecules 2023; 13:biom13020247. [PMID: 36830616 PMCID: PMC9953526 DOI: 10.3390/biom13020247] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Revised: 01/20/2023] [Accepted: 01/26/2023] [Indexed: 01/31/2023] Open
Abstract
Glycosaminoglycans (GAGs) are a class of linear anionic periodic polysaccharides containing disaccharide repetitive units. These molecules interact with a variety of proteins in the extracellular matrix and so participate in biochemically crucial processes such as cell signalling affecting tissue regeneration as well as the onset of cancer, Alzheimer's or Parkinson's diseases. Due to their flexibility, periodicity and chemical heterogeneity, often termed "sulfation code", GAGs are challenging molecules both for experiments and computation. One of the key questions in the GAG research is the specificity of their intermolecular interactions. In this study, we make a step forward to deciphering the "sulfation code" of chondroitin sulfates-4,6 (CS4, CS6, where the numbers correspond to the position of sulfation in NAcGal residue) and dermatan sulfate (DS), which is different from CSs by the presence of IdoA acid instead of GlcA. We rigorously investigate two sets of these GAGs in dimeric, tetrameric and hexameric forms with molecular dynamics-based descriptors. Our data clearly suggest that CS4, CS6 and DS are substantially different in terms of their structural, conformational and dynamic properties, which contributes to the understanding of how these molecules can be different when they bind proteins, which could have practical implications for the GAG-based drug design strategies in the regenerative medicine.
Collapse
|
11
|
Govind Kumar V, Polasa A, Agrawal S, Kumar TKS, Moradi M. Binding affinity estimation from restrained umbrella sampling simulations. NATURE COMPUTATIONAL SCIENCE 2023; 3:59-70. [PMID: 38177953 PMCID: PMC10766565 DOI: 10.1038/s43588-022-00389-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2022] [Accepted: 12/05/2022] [Indexed: 01/06/2024]
Abstract
The protein-ligand binding affinity quantifies the binding strength between a protein and its ligand. Computer modeling and simulations can be used to estimate the binding affinity or binding free energy using data- or physics-driven methods or a combination thereof. Here we discuss a purely physics-based sampling approach based on biased molecular dynamics simulations. Our proposed method generalizes and simplifies previously suggested stratification strategies that use umbrella sampling or other enhanced sampling simulations with additional collective-variable-based restraints. The approach presented here uses a flexible scheme that can be easily tailored for any system of interest. We estimate the binding affinity of human fibroblast growth factor 1 to heparin hexasaccharide based on the available crystal structure of the complex as the initial model and four different variations of the proposed method to compare against the experimentally determined binding affinity obtained from isothermal titration calorimetry experiments.
Collapse
Affiliation(s)
- Vivek Govind Kumar
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR, USA
| | - Adithya Polasa
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR, USA
| | - Shilpi Agrawal
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR, USA
| | | | - Mahmoud Moradi
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR, USA.
| |
Collapse
|
12
|
Kogut MM, Danielsson A, Ricard-Blum S, Samsonov SA. Impact of calcium ions on the structural and dynamic properties of heparin oligosaccharides by computational analysis. Comput Biol Chem 2022; 99:107727. [PMID: 35841830 DOI: 10.1016/j.compbiolchem.2022.107727] [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: 05/12/2022] [Revised: 06/24/2022] [Accepted: 07/07/2022] [Indexed: 11/13/2022]
Abstract
Heparin (HP) belongs to glycosaminoglycans (GAGs), anionic linear polysaccharides composed of repetitive disaccharide units. They are key players in many biological processes occurring in the extracellular matrix and at the cell surface. GAGs are challenging molecules for computational research due to their high chemical heterogeneity, flexibility, periodicity, pseudosymmetry, predominantly electrostatics-driven nature of interactions with their protein partners and potential multipose binding. The molecular mechanisms underlying GAG interactions mediated by divalent ions, which are important for GAG binding to several proteins, are not well understood. The goal of this study was to characterize the binding of Ca2+ to two HP oligosaccharides of different lengths (dp10 and dp18, dp: degree of polymerization) and their impact on HP conformational space and their dynamic behavior with the use of molecular dynamics (MD)-based approaches with two Ca2+ parameter sets. MD data suggested that the flexibility of the monosaccharides, the glycosidic linkages and ring puckering were not affected by the presence of Ca2+, in contrast to H-bond propensities and the calculated Rg for a fraction of the oligosaccharide populations in both dp10 and dp18. Moreover, the essential differences in the data obtained by using two Ca2+ parameter sets were reported.
Collapse
Affiliation(s)
- Małgorzata M Kogut
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Annemarie Danielsson
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Sylvie Ricard-Blum
- Univ Lyon, University Claude Bernard Lyon 1, CNRS, Institute of Molecular and Supramolecular Chemistry and Biochemistry, UMR 5246, Villeurbanne CEDEX F-69622, France
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, Gdańsk 80-308, Poland.
| |
Collapse
|
13
|
Marcisz M, Gaardløs M, Bojarski KK, Siebenmorgen T, Zacharias M, Samsonov SA. Explicit solvent repulsive scaling replica exchange molecular dynamics (RS-REMD) in molecular modeling of protein-glycosaminoglycan complexes. J Comput Chem 2022; 43:1633-1640. [PMID: 35796487 DOI: 10.1002/jcc.26965] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2022] [Revised: 06/20/2022] [Accepted: 06/22/2022] [Indexed: 11/10/2022]
Abstract
Glycosaminoglcyans (GAGs), linear anionic periodic polysaccharides, are crucial for many biologically relevant functions in the extracellular matrix. By interacting with proteins GAGs mediate processes such as cancer development, cell proliferation and the onset of neurodegenerative diseases. Despite this eminent importance of GAGs, they still represent a limited focus for the computational community in comparison to other classes of biomolecules. Therefore, there is a lack of modeling tools designed specifically for docking GAGs. One has to rely on existing docking software developed mostly for small drug molecules substantially differing from GAGs in their basic physico-chemical properties. In this study, we present an updated protocol for docking GAGs based on the Repulsive Scaling Replica Exchange Molecular Dynamics (RS-REMD) that includes explicit solvent description. The use of this water model improved docking performance both in terms of its accuracy and speed. This method represents a significant computational progress in GAG-related research.
Collapse
Affiliation(s)
- Mateusz Marcisz
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland.,Intercollegiate Faculty of Biotechnology, Universuty of Gdańsk and Medical University of Gdańsk, Gdańsk, Poland
| | | | - Krzysztof K Bojarski
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland.,Department of Physical Chemistry, Gdańsk University of Technology, Gdańsk, Poland
| | - Till Siebenmorgen
- Institute of Structural Biology, Helmholtz Zentrum München, Neuherberg, Germany
| | - Martin Zacharias
- Physics Department, Technical University of Munich, Garching, Germany
| | | |
Collapse
|
14
|
Molecular dynamics simulations to understand glycosaminoglycan interactions in the free- and protein-bound states. Curr Opin Struct Biol 2022; 74:102356. [DOI: 10.1016/j.sbi.2022.102356] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2021] [Revised: 02/08/2022] [Accepted: 02/13/2022] [Indexed: 11/18/2022]
|
15
|
Roy R, Jonniya NA, Kar P. Effect of Sulfation on the Conformational Dynamics of Dermatan Sulfate Glycosaminoglycan: A Gaussian Accelerated Molecular Dynamics Study. J Phys Chem B 2022; 126:3852-3866. [PMID: 35594147 DOI: 10.1021/acs.jpcb.2c01807] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Glycosaminoglycans (GAGs) are anionic biopolymers present on cell surfaces as a part of proteoglycans. The biological activities of GAGs depend on the sulfation pattern. In our study, we have considered three octadecasaccharide dermatan sulfate (DS) chains with increasing order of sulfation (dp6s, dp7s, and dp12s) to illuminate the role of sulfation on the GAG units and its chain conformation through 10 μs-long Gaussian accelerated molecular dynamics simulations. DS is composed of repeating disaccharide units of iduronic acid (IdoA) and N-acetylgalactosamine (N-GalNAc). Here, N-GalNAc is linked to IdoA via β(1-4), while IdoA is linked to N-GalNAc through α(1-3). With the increase in sulfation, the DS structure becomes more rigid and linear, as is evident from the distribution of root-mean-square deviations (RMSDs) and end-to-end distances. The tetrasaccharide linker region of the main chain shows a rigid conformation in terms of the glycosidic linkage. We have observed that upon sulfation (i.e., dp12s), the ring flip between two chair forms vanished for IdoA. The dynamic cross-correlation analysis reveals that the anticorrelation motions in dp12s are reduced significantly compared to dp6s or dp7s. An increase in sulfation generates relatively more stable hydrogen-bond networks, including water bridging with the neighboring monosaccharides. Despite the favorable linear structures of the GAG chains, our study also predicts few significant bendings related to the different puckering states, which may play a notable role in the function of the DS. The relation between the global conformation with the micro-level parameters such as puckering and water-mediated hydrogen bonds shapes the overall conformational space of GAGs. Overall, atomistic details of the DS chain provided in this study will help understand their functional and mechanical roles, besides developing new biomaterials.
Collapse
Affiliation(s)
- Rajarshi Roy
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| | - Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| |
Collapse
|
16
|
In-Depth Molecular Dynamics Study of All Possible Chondroitin Sulfate Disaccharides Reveals Key Insight into Structural Heterogeneity and Dynamism. Biomolecules 2022; 12:biom12010077. [PMID: 35053225 PMCID: PMC8773825 DOI: 10.3390/biom12010077] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2021] [Revised: 12/23/2021] [Accepted: 12/29/2021] [Indexed: 12/18/2022] Open
Abstract
GAGs exhibit a high level of conformational and configurational diversity, which remains untapped in terms of the recognition and modulation of proteins. Although GAGs are suggested to bind to more than 800 biologically important proteins, very few therapeutics have been designed or discovered so far. A key challenge is the inability to identify, understand and predict distinct topologies accessed by GAGs, which may help design novel protein-binding GAG sequences. Recent studies on chondroitin sulfate (CS), a key member of the GAG family, pinpointing its role in multiple biological functions led us to study the conformational dynamism of CS building blocks using molecular dynamics (MD). In the present study, we used the all-atom GLYCAM06 force field for the first time to explore the conformational space of all possible CS building blocks. Each of the 16 disaccharides was solvated in a TIP3P water box with an appropriate number of counter ions followed by equilibration and a production run. We analyzed the MD trajectories for torsional space, inter- and intra-molecular H-bonding, bridging water, conformational spread and energy landscapes. An in-house phi and psi probability density analysis showed that 1→3-linked sequences were more flexible than 1→4-linked sequences. More specifically, phi and psi regions for 1→4-linked sequences were held within a narrower range because of intra-molecular H-bonding between the GalNAc O5 atom and GlcA O3 atom, irrespective of sulfation pattern. In contrast, no such intra-molecular interaction arose for 1→3-linked sequences. Further, the stability of 1→4-linked sequences also arose from inter-molecular interactions involving bridged water molecules. The energy landscape for both classes of CS disaccharides demonstrated increased ruggedness as the level of sulfation increased. The results show that CS building blocks present distinct conformational dynamism that offers the high possibility of unique electrostatic surfaces for protein recognition. The fundamental results presented here will support the development of algorithms that help to design longer CS chains for protein recognition.
Collapse
|
17
|
Lutsyk V, Plazinski W. Conformational Properties of Glycosaminoglycan Disaccharides: A Molecular Dynamics Study. J Phys Chem B 2021; 125:10900-10916. [PMID: 34550710 DOI: 10.1021/acs.jpcb.1c04860] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
The structure and conformation of glycosaminoglycans (GAGs) are of central importance to understand the mechanisms behind their functions in biological systems. Due to the inherent chemical and structural heterogeneity of GAGs, focusing on longer, naturally existing GAG chains hinders drawing conclusions on the influence of the chemical functionalization on the basic conformational degree of freedom, that is, the dynamic shape of glycosidic linkage present in the particular disaccharide repeating unit. In the present study, we have considered the complete set of 106 GAG-related disaccharides, being potential building blocks for longer GAG chains (including hyaluronan, chondroitin, keratan, dermatan, and heparan). Both the unfunctionalized units and all possible combinations of either partially or fully sulfated derivatives contribute to this number. The unbiased and enhanced sampling molecular dynamics simulations provide a link to understand the influence of sulfation on the conformational properties of GAG glycosidic linkages. Residue-residue hydrogen bonding is not significant for either the glycosidic linkage conformation or its flexibility. It was found that in the majority of cases, the dominating conformation of the linkage is weakly affected by sulfation and the main role is played by the steric and stereoelectronic effects. However, there exist numerous cases where sulfation increases the contribution of alternative conformations to a nonnegligible extent and, in some rare cases (restricted to disaccharides building heparan), leads to the reorientation of the glycosidic linkage. The identified sulfation sites, being the most important in this context, are C6 and C3 at the GlcNAc residue. Finally, the full set of free energy maps relying on the glycosidic dihedral angle values for diverse GAG disaccharides are provided; they may be used for further studies, focused on longer GAG chains.
Collapse
Affiliation(s)
- Valery Lutsyk
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, Krakow 30-239, Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, Krakow 30-239, Poland
| |
Collapse
|
18
|
Marcisz M, Zacharias M, Samsonov SA. Modeling Protein-Glycosaminoglycan Complexes: Does the Size Matter? J Chem Inf Model 2021; 61:4475-4485. [PMID: 34494837 PMCID: PMC8479808 DOI: 10.1021/acs.jcim.1c00664] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
![]()
Docking glycosaminoglycans (GAGs) has been challenging because
of the complex nature of these long periodic linear and negatively
charged polysaccharides. Although standard docking tools like Autodock3
are successful when docking GAGs up to hexameric length, they experience
challenges to properly dock longer GAGs. Similar limitations concern
other docking approaches typically developed for docking ligands of
limited size to proteins. At the same time, most of more advanced
docking approaches are challenging for a user who is inexperienced
with complex in silico methodologies. In this work,
we evaluate the binding energies of complexes with different lengths
of GAGs using all-atom molecular dynamics simulations. Based on this
analysis, we propose a new docking protocol for long GAGs that consists
of conventional docking of short GAGs and further elongation with
the use of a coarse-grained representation of the GAG parts not being
in direct contact with its protein receptor. This method automated
by a simple script is straightforward to use within the Autodock3
framework but also useful in combination with other standard docking
tools. We believe that this method with some minor case-specific modifications
could also be used for docking other linear charged polymers.
Collapse
Affiliation(s)
- Mateusz Marcisz
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland.,Intercollegiate Faculty of Biotechnology of UG and MUG, ul. Abrahama 58, 80-307 Gdańsk, Poland
| | - Martin Zacharias
- Center of Functional Protein Assemblies, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
19
|
Künze G, Huster D, Samsonov SA. Investigation of the structure of regulatory proteins interacting with glycosaminoglycans by combining NMR spectroscopy and molecular modeling - the beginning of a wonderful friendship. Biol Chem 2021; 402:1337-1355. [PMID: 33882203 DOI: 10.1515/hsz-2021-0119] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2021] [Accepted: 04/09/2021] [Indexed: 11/15/2022]
Abstract
The interaction of regulatory proteins with extracellular matrix or cell surface-anchored glycosaminoglycans (GAGs) plays important roles in molecular recognition, wound healing, growth, inflammation and many other processes. In spite of their high biological relevance, protein-GAG complexes are significantly underrepresented in structural databases because standard tools for structure determination experience difficulties in studying these complexes. Co-crystallization with subsequent X-ray analysis is hampered by the high flexibility of GAGs. NMR spectroscopy experiences difficulties related to the periodic nature of the GAGs and the sparse proton network between protein and GAG with distances that typically exceed the detection limit of nuclear Overhauser enhancement spectroscopy. In contrast, computer modeling tools have advanced over the last years delivering specific protein-GAG docking approaches successfully complemented with molecular dynamics (MD)-based analysis. Especially the combination of NMR spectroscopy in solution providing sparse structural constraints with molecular docking and MD simulations represents a useful synergy of forces to describe the structure of protein-GAG complexes. Here we review recent methodological progress in this field and bring up examples where the combination of new NMR methods along with cutting-edge modeling has yielded detailed structural information on complexes of highly relevant cytokines with GAGs.
Collapse
Affiliation(s)
- Georg Künze
- Center for Structural Biology, Vanderbilt University, 465 21st Ave S, 5140 MRB3, Nashville, TN37240, USA.,Department of Chemistry, Vanderbilt University, 7330 Stevenson Center, Station B 351822, Nashville, TN37235, USA.,Institute for Drug Discovery, University of Leipzig, Brüderstr. 34, D-04103Leipzig, Germany
| | - Daniel Huster
- Institute for Medical Physics and Biophysics, University of Leipzig, Härtelstr. 16-18, D-04107Leipzig, Germany
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, Ul. Wita Stwosza 63, 80-308Gdańsk, Poland
| |
Collapse
|
20
|
Maszota-Zieleniak M, Marcisz M, Kogut MM, Siebenmorgen T, Zacharias M, Samsonov SA. Evaluation of replica exchange with repulsive scaling approach for docking glycosaminoglycans. J Comput Chem 2021; 42:1040-1053. [PMID: 33768554 DOI: 10.1002/jcc.26496] [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: 11/06/2020] [Revised: 01/19/2021] [Accepted: 01/27/2021] [Indexed: 12/14/2022]
Abstract
Glycosaminoglycans (GAGs), long linear periodic anionic polysaccharides, are key molecules in the extracellular matrix (ECM). Therefore, deciphering their role in the biologically relevant context is important for fundamental understanding of the processes ongoing in ECM and for establishing new strategies in the regenerative medicine. Although GAGs represent a number of computational challenges, molecular docking is a powerful tool for analysis of their interactions. Despite the recent development of GAG-specific docking approaches, there is plenty of room for improvement. Here, replica exchange molecular dynamics with repulsive scaling (REMD-RS) recently proved to be a successful approach for protein-protein complexes, was applied to dock GAGs. In this method, effective pairwise radii are increased in different Hamiltonian replicas. REMD-RS is shown to be an attractive alternative to classical docking approaches for GAGs. This work contributes to setting up of GAG-specific computational protocols and provides new insights into the nature of these biological systems.
Collapse
Affiliation(s)
| | | | | | - Till Siebenmorgen
- Physics Department, Technical University of Munich, Garching, Germany
| | - Martin Zacharias
- Physics Department, Technical University of Munich, Garching, Germany
| | | |
Collapse
|
21
|
Marcisz M, Huard B, Lipska AG, Samsonov SA. Further analyses of APRIL/APRIL-Receptor/Glycosaminoglycan interactions by biochemical assays linked to computational studies. Glycobiology 2021; 31:772-786. [PMID: 33682874 DOI: 10.1093/glycob/cwab016] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 02/23/2021] [Accepted: 02/23/2021] [Indexed: 12/13/2022] Open
Abstract
A proliferation-inducing ligand (APRIL) is a member of the tumor necrosis factor superfamily. APRIL is quite unique in this superfamily for at least for two reasons: i) it binds to glycosaminoglycans (GAGs) via its positively charged N-terminus; ii) one of its signaling receptor, the transmembrane activator CAML interactor (TACI) was also reported to bind GAGs. Here, as provided by biochemical evidences with the use of an APRIL deletion mutant linked to computational studies, APRIL-GAG interaction involved other regions than the APRIL N-terminus. Preferential interaction of APRIL with heparin followed by chondroitin sulfate E were confirmed by in silico analysis. Both computational and experimental approaches did not reveal heparan sulfate binding to TACI. Together, computational results corroborated experiments contributing with atomistic details to the knowledge on this biologically relevant trimolecular system. Additionally, a high-throughput rigorous analysis of the free energy calculations data was performed to critically evaluate the applied computational methodologies.
Collapse
Affiliation(s)
- Mateusz Marcisz
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland.,Intercollegiate Faculty of Biotechnology of UG and MUG, ul. Abrahama 58, 80-307 Gdańsk, Poland
| | - Bertrand Huard
- TIMC-IMAG, university Grenoble-Alpes, CNRS UMR 5525, La Tronche, France
| | - Agnieszka G Lipska
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
22
|
Yu M, Zhang T, Zhang W, Sun Q, Li H, Li JP. Elucidating the Interactions Between Heparin/Heparan Sulfate and SARS-CoV-2-Related Proteins-An Important Strategy for Developing Novel Therapeutics for the COVID-19 Pandemic. Front Mol Biosci 2021; 7:628551. [PMID: 33569392 PMCID: PMC7868326 DOI: 10.3389/fmolb.2020.628551] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Accepted: 12/16/2020] [Indexed: 12/12/2022] Open
Abstract
Owing to the high mortality and the spread rate, the infectious disease caused by SARS-CoV-2 has become a major threat to public health and social economy, leading to over 70 million infections and 1. 6 million deaths to date. Since there are currently no effective therapeutic or widely available vaccines, it is of urgent need to look for new strategies for the treatment of SARS-CoV-2 infection diseases. Binding of a viral protein onto cell surface heparan sulfate (HS) is generally the first step in a cascade of interaction that is required for viral entry and the initiation of infection. Meanwhile, interactions of selectins and cytokines (e.g., IL-6 and TNF-α) with HS expressed on endothelial cells are crucial in controlling the recruitment of immune cells during inflammation. Thus, structurally defined heparin/HS and their mimetics might serve as potential drugs by competing with cell surface HS for the prevention of viral adhesion and modulation of inflammatory reaction. In this review, we will elaborate coronavirus invasion mechanisms and summarize the latest advances in HS-protein interactions, especially proteins relevant to the process of coronavirus infection and subsequent inflammation. Experimental and computational techniques involved will be emphasized.
Collapse
Affiliation(s)
- Mingjia Yu
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, China
| | - Tianji Zhang
- Division of Chemistry and Analytical Science, National Institute of Metrology, Beijing, China
| | - Wei Zhang
- Division of Chemistry and Analytical Science, National Institute of Metrology, Beijing, China
| | - Qianyun Sun
- Division of Chemistry, Shandong Institute of Metrology, Jinan, China
| | - Hongmei Li
- Division of Chemistry and Analytical Science, National Institute of Metrology, Beijing, China
| | - Jin-ping Li
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, China
- Department of Medical Biochemistry and Microbiology, University of Uppsala, Uppsala, Sweden
| |
Collapse
|
23
|
Kogut MM, Maszota-Zieleniak M, Marcisz M, Samsonov SA. Computational insights into the role of calcium ions in protein–glycosaminoglycan systems. Phys Chem Chem Phys 2021; 23:3519-3530. [DOI: 10.1039/d0cp05438k] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The prediction power of computational methodologies for studying the role of ions in protein–glycosaminoglycan interactions was critically assessed.
Collapse
|
24
|
Bojarski KK, Samsonov SA. Role of Oligosaccharide Chain Polarity in Protein-Glycosaminoglycan Interactions. J Chem Inf Model 2020; 61:455-466. [PMID: 33375794 DOI: 10.1021/acs.jcim.0c01402] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Glycosaminoglycans (GAGs) are long unbranched anionic polysaccharides made up of repetitive disaccharide units involved in biologically relevant processes in the extracellular matrix such as cell proliferation and communication. A GAG can be bound in antiparallel energetically comparable orientations on the protein surface, and these orientations are, therefore, difficult to distinguish both experimentally and computationally. In this study, for the first time we analyzed the impact of the GAG chain polarity on the interactions with Fibroblast Growth Factors-1 and -2 (FGF-1 and FGF-2). We performed a series of 1 μs molecular dynamics simulations of the FGF-1 and FGF-2 complexes with heparin (HP), a GAG representative, of different length. We analyzed the relationship between the HP orientation, energetic, and conformational space characteristics of FGF-1-HP and FGF-2-HP complexes. We concluded that HP can be bound by these proteins in the same binding sites but in different orientations, while the orientation present in the experimental structure might be favorable. Our data presented in this study provide a novel view on the impact of GAG polarity on the specificity of protein-GAG complex formation, which is an essential aspect for the proper understanding of the intermolecular interactions in these systems.
Collapse
Affiliation(s)
- Krzysztof K Bojarski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| |
Collapse
|
25
|
Vallet SD, Clerc O, Ricard-Blum S. Glycosaminoglycan-Protein Interactions: The First Draft of the Glycosaminoglycan Interactome. J Histochem Cytochem 2020; 69:93-104. [PMID: 32757871 DOI: 10.1369/0022155420946403] [Citation(s) in RCA: 100] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The six mammalian glycosaminoglycans (GAGs), chondroitin sulfate, dermatan sulfate, heparin, heparan sulfate, hyaluronan, and keratan sulfate, are linear polysaccharides. Except for hyaluronan, they are sulfated to various extent, and covalently attached to proteins to form proteoglycans. GAGs interact with growth factors, morphogens, chemokines, extracellular matrix proteins and their bioactive fragments, receptors, lipoproteins, and pathogens. These interactions mediate their functions, from embryonic development to extracellular matrix assembly and regulation of cell signaling in various physiological and pathological contexts such as angiogenesis, cancer, neurodegenerative diseases, and infections. We give an overview of GAG-protein interactions (i.e., specificity and chemical features of GAG- and protein-binding sequences), and review the available GAG-protein interaction networks. We also provide the first comprehensive draft of the GAG interactome composed of 832 biomolecules (827 proteins and five GAGs) and 932 protein-GAG interactions. This network is a scaffold, which in the future should integrate structures of GAG-protein complexes, quantitative data of the abundance of GAGs in tissues to build tissue-specific interactomes, and GAG interactions with metal ions such as calcium, which plays a major role in the assembly of the extracellular matrix and its interactions with cells. This contextualized interactome will be useful to identify druggable GAG-protein interactions for therapeutic purpose.
Collapse
Affiliation(s)
- Sylvain D Vallet
- Univ Lyon, University Claude Bernard Lyon 1, CNRS, INSA Lyon, CPE, Institute of Molecular and Supramolecular Chemistry and Biochemistry, UMR 5246, Villeurbanne Cedex, France
| | - Olivier Clerc
- Univ Lyon, University Claude Bernard Lyon 1, CNRS, INSA Lyon, CPE, Institute of Molecular and Supramolecular Chemistry and Biochemistry, UMR 5246, Villeurbanne Cedex, France
| | - Sylvie Ricard-Blum
- Univ Lyon, University Claude Bernard Lyon 1, CNRS, INSA Lyon, CPE, Institute of Molecular and Supramolecular Chemistry and Biochemistry, UMR 5246, Villeurbanne Cedex, France
| |
Collapse
|
26
|
Bojarski KK, Karczyńska AS, Samsonov SA. Role of Glycosaminoglycans in Procathepsin B Maturation: Molecular Mechanism Elucidated by a Computational Study. J Chem Inf Model 2020; 60:2247-2256. [PMID: 32155059 PMCID: PMC7588040 DOI: 10.1021/acs.jcim.0c00023] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
![]()
Procathepsins
are an inactive, immature form of cathepsins, predominantly
cysteine proteases present in the extracellular matrix (ECM) and in
lysosomes that play a key role in various biological processes such
as bone resorption or intracellular proteolysis. The enzymatic activity
of cathepsins can be mediated by glycosaminoglycans (GAGs), long unbranched
periodic negatively charged polysaccharides found in ECM that take
part in many biological processes such as anticoagulation, angiogenesis,
and tissue regeneration. In addition to the known effects on mature
cathepsins, GAGs can mediate the maturation process of procathepsins,
in particular, procathepsin B. However, the detailed mechanism of
this mediation at the molecular level is still unknown. In this study,
for the first time, we aimed to unravel the role of GAGs in this process
using computational approaches. We rigorously analyzed procathepsin
B–GAG complexes in terms of their dynamics, energetics, and
potential allosteric regulation. We revealed that GAGs can stabilize
the conformation of the procathepsin B structure with the active site
accessible for the substrate and concluded that GAGs most probably
bind to procathepsin B once the zymogen adopts the enzymatically active
conformation. Our data provided a novel mechanistic view of the maturation
process of procathepsin B, while the approaches elaborated here might
be useful to study other procathepsins. Furthermore, our data can
serve as a rational guide for experimental work on procathepsin–GAG
systems that are not characterized in vivo and in vitro yet.
Collapse
Affiliation(s)
- Krzysztof K Bojarski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | | | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
27
|
Rat cathepsin K: Enzymatic specificity and regulation of its collagenolytic activity. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2019; 1868:140318. [PMID: 31740411 DOI: 10.1016/j.bbapap.2019.140318] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Revised: 11/08/2019] [Accepted: 11/11/2019] [Indexed: 12/25/2022]
Abstract
Human cathepsin K (hCatK), which is highly expressed in osteoclasts, has the noteworthy ability to cleave type I and II collagens in their helical domain. Its collagenase potency depends strictly on the formation of an oligomeric complex with chondroitin 4-sulfate (C4-S). Accordingly, hCatK is a pivotal protease involved in bone resorption and is an attractive target for the treatment of osteoporosis. As rat is a common animal model for the evaluation of hCatK inhibitors, we conducted a comparative analysis of rat CatK (rCatK) and hCatK, which share a high degree of identity (88%) and similarity (93%). The pH activity profile of both enzymes displayed a similar bell-shaped curve (optimal pH: 6.4). Presence of Ser134 and Val160 in the S2 pocket of rCatK instead of Ala and Leu residues, respectively, in hCatK, led to a weaker peptidase activity, as observed for mouse CatK. Also, regardless of the presence of C4-S, rCatK cleaved in the nonhelical telopeptide regions of both type I (tail) and type II (articular joint) rat collagens. Structure-based computational analyses (electrostatic potential, molecular docking, molecular dynamics, free energy calculations) sustained that the C4-S mediated collagenolytic activity of rCatK obeys distinct molecular interactions from those of hCatK. Additionally, T-kininogen (a.k.a. thiostatin), a unique rat serum acute phase molecule, acted as a tight-binding inhibitor of hCatK (Ki = 0.11 ± 0.05 nM). Taken into account the increase of T-Kininogen level in inflamed rat sera, this may raise the question of the appropriateness to evaluate pharmacological hCatK inhibitors in this peculiar animal model.
Collapse
|
28
|
Analysis of Procollagen C-Proteinase Enhancer-1/Glycosaminoglycan Binding Sites and of the Potential Role of Calcium Ions in the Interaction. Int J Mol Sci 2019; 20:ijms20205021. [PMID: 31658765 PMCID: PMC6829435 DOI: 10.3390/ijms20205021] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2019] [Revised: 10/07/2019] [Accepted: 10/09/2019] [Indexed: 12/15/2022] Open
Abstract
In this study, we characterize the interactions between the extracellular matrix protein, procollagen C-proteinase enhancer-1 (PCPE-1), and glycosaminoglycans (GAGs), which are linear anionic periodic polysaccharides. We applied molecular modeling approaches to build a structural model of full-length PCPE-1, which is not experimentally available, to predict GAG binding poses for various GAG lengths, types and sulfation patterns, and to determine the effect of calcium ions on the binding. The computational data are analyzed and discussed in the context of the experimental results previously obtained using surface plasmon resonance binding assays. We also provide experimental data on PCPE-1/GAG interactions obtained using inhibition assays with GAG oligosaccharides ranging from disaccharides to octadecasaccharides. Our results predict the localization of GAG-binding sites at the amino acid residue level onto PCPE-1 and is the first attempt to describe the effects of ions on protein-GAG binding using modeling approaches. In addition, this study allows us to get deeper insights into the in silico methodology challenges and limitations when applied to GAG-protein interactions.
Collapse
|
29
|
Samsonov SA, Lubecka EA, Bojarski KK, Ganzynkowicz R, Liwo A. Local and long range potentials for heparin‐protein systems for coarse‐grained simulations. Biopolymers 2019; 110:e23269. [DOI: 10.1002/bip.23269] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 02/21/2019] [Accepted: 02/22/2019] [Indexed: 12/31/2022]
Affiliation(s)
| | - Emilia A. Lubecka
- Faculty of ChemistryUniversity of Gdańsk Gdańsk Poland
- Faculty of Mathematics, Physics and Informatics, Institute of InformaticsUniversity of Gdańsk Gdańsk Poland
| | | | | | - Adam Liwo
- Faculty of ChemistryUniversity of Gdańsk Gdańsk Poland
| |
Collapse
|