Rational design of glycosaminoglycan binding cyclic peptides using cPEPmatch.

Cyclic peptides present a robust platform for drug design, offering high specificity and stability due to their conformationally constrained structures. In this study, we introduce an updated version of the Cyclic Peptide Matching program (cPEPmatch) tailored for the identification of cyclic peptides capable of mimicking protein-glycosaminoglycan (GAG) binding sites. We focused on engineering cyclic peptides to replicate the GAG-binding affinity of antithrombin III (ATIII), a protein that plays a crucial role in modulating anticoagulation through interaction with the GAG heparin. By integrating computational and experimental methods, we successfully identified a cyclic peptide binder with promising potential for future optimization. MD simulations and MM-GBSA calculations were used to assess binding efficacy, supplemented by umbrella sampling to approximate free energy landscapes. The binding specificity was further validated through NMR and ITC experiments. Our findings demonstrate that the computationally designed cyclic peptides effectively target GAGs, suggesting their potential as novel therapeutic agents. This study advances our understanding of peptide-GAG interactions and lays the groundwork for future development of cyclic peptide-based therapeutics.

Detecting Protein-Ligand Interactions with Nitroxide Based Paramagnetic Cosolutes.

NMR spectroscopy techniques can provide important information about protein-ligand interactions. Here we tested an NMR approach which relies on the measurement of paramagnetic relaxation enhancements (PREs) arising from analogous cationic, anionic or neutral soluble nitroxide molecules, which distribute around the protein-ligand complex depending on near-surface electrostatic potentials. We applied this approach to two protein-ligand systems, interleukin-8 interacting with highly charged glycosaminoglycans and the SH2 domain of Grb2 interacting with less charged phospho-tyrosine tripeptides. The electrostatic potential around interleukin-8 and its changes upon binding of glycosaminoglycans could be derived from the PRE data and confirmed by theoretical predictions from Poisson-Boltzmann calculations. The ligand influence on the PREs and NMR-derived electrostatic potentials of Grb2 SH2 was localized to a narrow protein region which allowed the localization of the peptide binding pocket. Our analysis suggests that experiments with nitroxide cosolutes can be useful for investigating protein-ligand electrostatic interactions and mapping ligand binding sites.

Computational modeling of the molecular basis for the calcium-dependence of the mannuronan C-5 epimerase AvAlgE6 from Azotobacter vinelandii.

The mannuronan C-5 epimerases catalyze epimerization of ß-d-mannuronic acid to α-l-guluronic acid in alginate polymers. The seven extracellular Azotobacter vinelandii epimerases (AvAlgE1-7) are calcium-dependent, and calcium is essential for the structural integrity of their carbohydrate binding R-modules. Ca2+ is also found in the crystal structures of the A-modules, where it is suggested to play a structural role. In this study, the structure of the catalytic A-module of the A. vinelandii mannuronan C-5 epimerase AvAlgE6 is used to investigate the role of this Ca2+. Molecular dynamics (MD) simulations with and without calcium reveal the possible importance of the bound Ca2+ in the hydrophobic packing of ß-sheets. In addition, a putative calcium binding site is found in the active site, indicating a potential direct role of this calcium in the catalysis. According to the literature, two of the residues coordinating calcium in this site are essential for the activity. MD simulations of the interaction with bound substrate indicate that the presence of a calcium ion in this binding site increases the binding strength. Further, explicit calculations of the substrate dissociation pathways with umbrella sampling simulations show and energetically higher dissociation barrier when calcium is present. The present study eludes to a putative catalytic role of calcium in the charge neutralizing first step of the enzymatic reaction. In addition to the importance for understanding these enzymes' molecular mechanisms, this could have implications for engineering strategies of the epimerases in industrial alginate processing.

Mannuronate C-5 epimerases and their use in alginate modification.

Alginate is a polysaccharide consisting of ß-D-mannuronate (M) and α-L-guluronate (G) produced by brown algae and some bacterial species. Alginate has a wide range of industrial and pharmaceutical applications, owing mainly to its gelling and viscosifying properties. Alginates with high G content are considered more valuable since the G residues can form hydrogels with divalent cations. Alginates are modified by lyases, acetylases, and epimerases. Alginate lyases are produced by alginate-producing organisms and by organisms that use alginate as a carbon source. Acetylation protects alginate from lyases and epimerases. Following biosynthesis, alginate C-5 epimerases convert M to G residues at the polymer level. Alginate epimerases have been found in brown algae and alginate-producing bacteria, predominantly Azotobacter and Pseudomonas species. The best characterised epimerases are the extracellular family of AlgE1-7 from Azotobacter vinelandii(Av). AlgE1-7 all consist of combinations of one or two catalytic A-modules and one to seven regulatory R-modules, but even though they are sequentially and structurally similar, they create different epimerisation patterns. This makes the AlgE enzymes promising for tailoring of alginates to have the desired properties. The present review describes the current state of knowledge regarding alginate-active enzymes with focus on epimerases, characterisation of the epimerase reaction, and how alginate epimerases can be used in alginate production.

Studying specificity in protein-glycosaminoglycan recognition with umbrella sampling.

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.

Rapamycin-induced autophagy in osteosarcoma cells is mediated via the biglycan/Wnt/ß-catenin signaling axis.

Biglycan is a class I secreted small leucine-rich proteoglycan (SLRP), which regulates signaling pathways connected to bone pathologies. Autophagy is a vital catabolic process with a dual role in cancer progression. Here, we show that biglycan inhibits autophagy in two osteosarcoma cell lines (P ≤ 0.001), while rapamycin-induced autophagy decreases biglycan expression in MG63 osteosarcoma cells and abrogates the biglycan-induced cell growth increase (P ≤ 0.001). Rapamycin also inhibits ß-catenin translocation to the nucleus, inhibiting the Wnt pathway (P ≤ 0.001) and reducing biglycan's colocalization with the Wnt coreceptor LRP6 (P ≤ 0.05). Furthermore, biglycan exhibits protective effects against the chemotherapeutic drug doxorubicin in MG63 OS cells through an autophagy-dependent manner (P ≤ 0.05). Cotreatment of these cells with rapamycin and doxorubicin enhances cells response to doxorubicin by decreasing biglycan (P ≤ 0.001) and ß-catenin (P ≤ 0.05) expression. Biglycan deficiency leads to increased caspase-3 activation (P ≤ 0.05), suggesting increased apoptosis of biglycan-deficient cells treated with doxorubicin. Computational models of LRP6 and biglycan complexes suggest that biglycan changes the receptor's ability to interact with other signaling molecules by affecting the interdomain bending angles in the receptor structure. Biglycan binding to LRP6 activates the Wnt pathway and ß-catenin nuclear translocation by disrupting ß-catenin degradation complex (P ≤ 0.01 and P ≤ 0.05). Interestingly, this mechanism is not followed in moderately differentiated, biglycan-nonexpressing U-2OS OS cells. To sum up, biglycan exhibits protective effects against the doxorubicin in MG63 OS cells by activating the Wnt signaling pathway and inhibiting autophagy.

Explicit solvent repulsive scaling replica exchange molecular dynamics (RS-REMD) in molecular modeling of protein-glycosaminoglycan complexes.

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.

Biglycan Interacts with Type I Insulin-like Receptor (IGF-IR) Signaling Pathway to Regulate Osteosarcoma Cell Growth and Response to Chemotherapy.

Osteosarcoma (OS) is a mesenchymally derived, aggressive bone cancer. OS cells produce an aberrant nonmineralized or partly mineralized extracellular matrix (ECM) whose components participate in signaling pathways connected to specific pathogenic phenotypes of this bone cancer. The expression of biglycan (BGN), a secreted small leucine-rich proteoglycan (SLRP), is correlated to aggressive OS phenotype and resistance to chemotherapy. A constitutive signaling of IGF-IR signaling input in sarcoma progression has been established. Here, we show that biglycan activates the IGF-IR signaling pathway to promote MG63 biglycan-secreting OS cell growth by forming a complex with the receptor. Computational models of IGF-IR and biglycan docking suggest that biglycan binds IGF-IR dimer via its concave surface. Our binding free energy calculations indicate the formation of a stable complex. Biglycan binding results in prolonged IGF-IR activation leading to protracted IGF-IR-dependent cell growth response of the poorly-differentiated MG63 cells. Moreover, biglycan facilitates the internalization (p ≤ 0.01, p ≤ 0.001) and sumoylation-enhanced nuclear translocation of IGF-IR (p ≤ 0.05) and its DNA binding in MG63 cells (p ≤ 0.001). The tyrosine kinase activity of the receptor mediates this mechanism. Furthermore, biglycan downregulates the expression of the tumor-suppressor gene, PTEN (p ≤ 0.01), and increases the expression of endothelial-mesenchymal transition (EMT) and aggressiveness markers vimentin (p ≤ 0.01) and fibronectin (p ≤ 0.01) in MG63 cells. Interestingly, this mechanism is not valid in moderately and well-differentiated, biglycan non-expressing U-2OS and Saos-2 OS cells. Furthermore, biglycan exhibits protective effects against the chemotherapeutic drug, doxorubicin, in MG63 OS cells (p ≤ 0.01). In conclusion, these data indicate a potential direct and adjunct therapeutical role of biglycan in osteosarcoma.

Mechanistic Basis for Understanding the Dual Activities of the Bifunctional Azotobacter vinelandii Mannuronan C-5-Epimerase and Alginate Lyase AlgE7.

The structure and functional properties of alginates are dictated by the monomer composition and molecular weight distribution. Mannuronan C-5-epimerases determine the monomer composition by catalyzing the epimerization of ß-d-mannuronic acid (M) residues into α-l-guluronic acid (G) residues. The molecular weight is affected by alginate lyases, which catalyze a ß-elimination mechanism that cleaves alginate chains. The reaction mechanisms for the epimerization and lyase reactions are similar, and some enzymes can perform both reactions. These dualistic enzymes share high sequence identity with mannuronan C-5-epimerases without lyase activity. The mechanism behind their activity and the amino acid residues responsible for it are still unknown. We investigate mechanistic determinants involved in the bifunctional epimerase and lyase activity of AlgE7 from Azotobacter vinelandii. Based on sequence analyses, a range of AlgE7 variants were constructed and subjected to activity assays and product characterization by nuclear magnetic resonance (NMR) spectroscopy. Our results show that calcium promotes lyase activity, whereas NaCl reduces the lyase activity of AlgE7. By using defined polymannuronan (polyM) and polyalternating alginate (polyMG) substrates, the preferred cleavage sites of AlgE7 were found to be M|XM and G|XM, where X can be either M or G. From the study of AlgE7 mutants, R148 was identified as an important residue for the lyase activity, and the point mutant R148G resulted in an enzyme with only epimerase activity. Based on the results obtained in the present study, we suggest a unified catalytic reaction mechanism for both epimerase and lyase activities where H154 functions as the catalytic base and Y149 functions as the catalytic acid. IMPORTANCE Postharvest valorization and upgrading of algal constituents are promising strategies in the development of a sustainable bioeconomy based on algal biomass. In this respect, alginate epimerases and lyases are valuable enzymes for tailoring the functional properties of alginate, a polysaccharide extracted from brown seaweed with numerous applications in food, medicine, and material industries. By providing a better understanding of the catalytic mechanism and of how the two enzyme actions can be altered by changes in reaction conditions, this study opens further applications of bacterial epimerases and lyases in the enzymatic tailoring of alginate polymers.

Insights into the roles of charged residues in substrate binding and mode of action of mannuronan C-5 epimerase AlgE4.

Mannuronan C-5 epimerases catalyze the epimerization of monomer residues in the polysaccharide alginate, changing the physical properties of the biopolymer. The enzymes are utilized to tailor alginate to numerous biological functions by alginate-producing organisms. The underlying molecular mechanism that control the processive movement of the epimerase along the substrate chain is still elusive. To study this, we have used an interdisciplinary approach combining molecular dynamics simulations with experimental methods from mutant studies of AlgE4, where initial epimerase activity and product formation were addressed with nuclear magnetic resonance spectroscopy, and characteristics of enzyme-substrate interactions were obtained with isothermal titration calorimetry and optical tweezers. Positive charges lining the substrate-binding groove of AlgE4 appear to control the initial binding of poly-mannuronate, and binding also seems to be mediated by both electrostatic and hydrophobic interactions. After the catalytic reaction, negatively charged enzyme residues might facilitate dissociation of alginate from the positive residues, working like electrostatic switches, allowing the substrate to translocate in the binding groove. Molecular simulations show translocation increments of two monosaccharide units before the next productive binding event resulting in mannuronate and guluronate (MG)-block formation, with the epimerase moving with its N-terminus towards the reducing end of the alginate chain. Our results indicate that the charge pair R343-D345 might be directly involved in conformational changes of a loop that can be important for binding and dissociation. The computational and experimental approaches used in this study complement each other, allowing for a better understanding of individual residues' roles in binding and movement along the alginate chains.

Identification of a Pivotal Residue for Determining the Block Structure-Forming Properties of Alginate C-5 Epimerases.

Alginate is a linear copolymer composed of 1→4 linked ß-d-mannuronic acid (M) and its epimer α-l-guluronic acid (G). The polysaccharide is first produced as homopolymeric mannuronan and subsequently, at the polymer level, C-5 epimerases convert M residues to G residues. The bacterium Azotobacter vinelandii encodes a family of seven secreted and calcium ion-dependent mannuronan C-5 epimerases (AlgE1-AlgE7). These epimerases consist of two types of structural modules: the A-modules, which contain the catalytic site, and the R-modules, which influence activity through substrate and calcium binding. In this study, we rationally designed new hybrid mannuronan C-5 epimerases constituting the A-module from AlgE6 and the R-module from AlgE4. This led to a better understanding of the molecular mechanism determining differences in MG- and GG-block-forming properties of the enzymes. A long loop with either tyrosine or phenylalanine extruding from the ß-helix of the enzyme proved essential in defining the final alginate block structure, probably by affecting substrate binding. Normal mode analysis of the A-module from AlgE6 supports the results.