Human cytokine and coronavirus nucleocapsid protein interactivity using large-scale virtual screens

Understanding the interactions between SARS-CoV-2 and the human immune system is paramount to the characterization of novel variants as the virus co-evolves with the human host. In this study, we employed state-of-the-art molecular docking tools to conduct large-scale virtual screens, predicting the binding affinities between 64 human cytokines against 17 nucleocapsid proteins from six betacoronaviruses. Our comprehensive in silico analyses reveal specific changes in cytokine-nucleocapsid protein interactions, shedding light on potential modulators of the host immune response during infection. These findings offer valuable insights into the molecular mechanisms underlying viral pathogenesis and may guide the future development of targeted interventions. This manuscript serves as insight into the comparison of deep learning based AlphaFold2-Multimer and the semi-physicochemical based HADDOCK for protein-protein docking. We show the two methods are complementary in their predictive capabilities. We also introduce a novel algorithm for rapidly assessing the binding interface of protein-protein docks using graph edit distance: graph-based interface residue assessment function (GIRAF). The high-performance computational framework presented here will not only aid in accelerating the discovery of effective interventions against emerging viral threats, but extend to other applications of high throughput protein-protein screens.

Sulfation Effects of Chondroitin Sulfate to Bind a Chemokine in Aqueous Medium: Conformational Heterogeneity and Dynamics from Molecular Simulation

The sulfation patterns and degree of sulfation of chondroitin sulfate (CS), an important class of glycosaminoglycans (GAG), and their interactions with chemokines are accountable for various diseases. To realize the underlying mechanism of such complex biological phenomena at a molecular level and their application in rational drug design, a study on conformations and dynamics of CSs is necessary. To explore this, in this study, we performed a series of atomistic molecular dynamics (MD) simulations with different sulfated variants of octadecasaccharide CS, like CS-C, CS-E, and CS-T, in their free forms and when bound to the protein chemokine CXCL8 dimer in an aqueous medium. The calculated binding free energy of CSs with the CXCL8 dimer is favorable, and the degree of sulfation favors the complexation process further with prominent hydrophobic and hydrogen-bonded interactions. We find that the recognition is associated with the configurational entropy loss of the CS molecules as calculated from the Gaussian mixture approach, which supports that the degree of sulfation regulates the process. Cluster analysis through the k-means algorithm and end-to-end distance measurement revealed that although the free CS molecules adopted linear conformations, the nonlinear conformations during binding with protein were noted. Adaptation of nonlinear forms in the bound forms is noteworthy for the less-sulfated CS-C and CS-E. Apart from favorable 4C1 conformations, the occasional appearance of skew-boat forms from the free-energy map of ring pucker for the GlcUA unit was observed, which remains unaffected by the sulfation. We find that during recognition, the average relaxation time of intra-CS and inter-CS-CXCL8 hydrogen bonds (HBs) is about a magnitude lesser than that of CS-water HBs, most prominent on the involvement of higher sulfated CS-T analogues. The translational motion of surrounded water molecules in CSs exhibited sublinear diffusion, and the degree of sublinearity increases around the heavily sulfated molecules due to the hindrance created by them as well as the presence of the chemokine and exhibited markedly slow heterogeneous diffusion.

Pharmacology of Heparin and Related Drugs: An Update

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.

Atomic-Resolution Experimental Structural Biology and Molecular Dynamics Simulations of Hyaluronan and Its Complexes

This review summarizes the atomic-resolution structural biology of hyaluronan and its complexes available in the Protein Data Bank, as well as published studies of atomic-resolution explicit-solvent molecular dynamics simulations on these and other hyaluronan and hyaluronan-containing systems. Advances in accurate molecular mechanics force fields, simulation methods and software, and computer hardware have supported a recent flourish in such simulations, such that the simulation publications now outnumber the structural biology publications by an order of magnitude. In addition to supplementing the experimental structural biology with computed dynamic and thermodynamic information, the molecular dynamics studies provide a wealth of atomic-resolution information on hyaluronan-containing systems for which there is no atomic-resolution structural biology either available or possible. Examples of these summarized in this review include hyaluronan pairing with other hyaluronan molecules and glycosaminoglycans, with ions, with proteins and peptides, with lipids, and with drugs and drug-like molecules. Despite limitations imposed by present-day computing resources on system size and simulation timescale, atomic-resolution explicit-solvent molecular dynamics simulations have been able to contribute significant insight into hyaluronan's flexibility and capacity for intra- and intermolecular non-covalent interactions.

Protease-Triggered Release of Stabilized CXCL12 from Coated Scaffolds in an Ex Vivo Wound Model

Biomaterials are designed to improve impaired healing of injured tissue. To accomplish better cell integration, we suggest to coat biomaterial surfaces with bio-functional proteins. Here, a mussel-derived surface-binding peptide is used and coupled to CXCL12 (stromal cell-derived factor 1α), a chemokine that activates CXCR4 and consequently recruits tissue-specific stem and progenitor cells. CXCL12 variants with either non-releasable or protease-mediated-release properties were designed and compared. Whereas CXCL12 was stabilized at the N-terminus for protease resistance, a C-terminal linker was designed that allowed for specific cleavage-mediated release by matrix metalloproteinase 9 and 2, since both enzymes are frequently found in wound fluid. These surface adhesive CXCL12 derivatives were produced by expressed protein ligation. Functionality of the modified chemokines was assessed by inositol phosphate accumulation and cell migration assays. Increased migration of keratinocytes and primary mesenchymal stem cells was demonstrated. Immobilization and release were studied for bioresorbable PCL-co-LC scaffolds, and accelerated wound closure was demonstrated in an ex vivo wound healing assay on porcine skin grafts. After 24 h, a significantly improved CXCL12-specific growth stimulation of the epithelial tips was already observed. The presented data display a successful application of protein-coated biomaterials for skin regeneration.

Modeling Protein-Glycosaminoglycan Complexes: Does the Size Matter?

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.

Improvement of wound healing by the development of ECM-inspired biomaterial coatings and controlled protein release

Implant design has evolved from biochemically inert substrates, minimizing cell and protein interaction, towards sophisticated bioactive substrates, modulating the host response and supporting the regeneration of the injured tissue. Important aspects to consider are the control of cell adhesion, the discrimination of bacteria and non-local cells from the desired tissue cell type, and the stimulation of implant integration and wound healing. Here, the extracellular matrix acts as a role model providing us with inspiration for sophisticated designs. Within this scope, small bioactive peptides have proven to be miscellaneously deployable for the mediation of surface, cell and matrix interactions. Combinations of adhesion ligands, proteoglycans, and modulatory proteins should guide multiple aspects of the regeneration process and cooperativity between the different extracellular matrix components, which bears the chance to maximize the therapeutic efficiency and simultaneously lower the doses. Hence, efforts to include multiple of these factors in biomaterial design are well worth. In the following, multifunctional implant coatings based on bioactive peptides are reviewed and concepts to implement strong surface anchoring for stable cell adhesion and a dynamic delivery of modulator proteins are discussed.

New insights into the role of glycosaminoglycans in the endosteal bone microenvironment

The bone microenvironment is a complex tissue in which heterogeneous cell populations of hematopoietic and mesenchymal origin interact with environmental cues to maintain tissue integrity. Both cellular and matrix components are subject to physiologic challenges and can dynamically respond by modifying cell/matrix interactions. When either component is impaired, the physiologic balance is lost. Here, we review the current state of knowledge of how glycosaminoglycans - organic components of the bone extracellular matrix - influence the bone micromilieu. We point out how they interact with mediators of distinct signaling pathways such as the RANKL/OPG axis, BMP and WNT signaling, and affect the activity of bone remodeling cells within the endosteal niche summarizing their potential for therapeutic intervention.

Insights into structure, affinity, specificity, and function of GAG-protein interactions through the chemoenzymatic preparation of defined sulfated oligohyaluronans

High amounts of glycosaminoglycans (GAG) such as hyaluronan (HA) occur in connective tissues. There is nowadays increasing evidence that a "sulfation code" exists which mediates numerous GAG functions. High molecular weight and inhomogeneity of GAG, however, aggravated detailed studies. Thus, synthetic oligosaccharides were urgently required. We will review here chemoenzymatic and analytic strategies to provide defined sulfated and anomerically modified GAG oligosaccharides of the HA type. Representative studies of protein/GAG interactions by (bio)chemical and biophysical methods are reported yielding novel insights into GAG-protein binding. Finally, the biological conclusions and in vivo applications of defined sulfated GAG oligosaccharides will be discussed.

Tuning the network charge of biohybrid hydrogel matrices to modulate the release of SDF-1

The delivery of chemotactic signaling molecules via customized biomaterials can effectively guide the migration of cells to improve the regeneration of damaged or diseased tissues. Here, we present a novel biohybrid hydrogel system containing two different sulfated glycosaminoglycans (sGAG)/sGAG derivatives, namely either a mixture of short heparin polymers (Hep-Mal) or structurally defined nona-sulfated tetrahyaluronans (9s-HA4-SH), to precisely control the release of charged signaling molecules. The polymer networks are described in terms of their negative charge, i.e. the anionic sulfate groups on the saccharides, using two parameters, the integral density of negative charge and the local charge distribution (clustering) within the network. The modulation of both parameters was shown to govern the release characteristics of the chemotactic signaling molecule SDF-1 and allows for seamless transitions between burst and sustained release conditions as well as the precise control over the total amount of delivered protein. The obtained hydrogels with well-adjusted release profiles effectively promote MSC migration in vitro and emerge as promising candidates for new treatment modalities in the context of bone repair and wound healing.

Investigation of the structure of regulatory proteins interacting with glycosaminoglycans by combining NMR spectroscopy and molecular modeling - the beginning of a wonderful friendship

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.

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

Mannuronan C-5 epimerases catalyse 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 mechanisms 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 NMR 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 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.

NMR Characterization of the Interactions Between Glycosaminoglycans and Proteins

Glycosaminoglycans (GAGs) constitute a considerable fraction of the glycoconjugates found on cellular membranes and in the extracellular matrix of virtually all mammalian tissues. The essential role of GAG-protein interactions in the regulation of physiological processes has been recognized for decades. However, the underlying molecular basis of these interactions has only emerged since 1990s. The binding specificity of GAGs is encoded in their primary structures, but ultimately depends on how their functional groups are presented to a protein in the three-dimensional space. This review focuses on the application of NMR spectroscopy on the characterization of the GAG-protein interactions. Examples of interpretation of the complex mechanism and characterization of structural motifs involved in the GAG-protein interactions are given. Selected families of GAG-binding proteins investigated using NMR are also described.

Targeting Chemokine-Glycosaminoglycan Interactions to Inhibit Inflammation

Leukocyte migration into tissues depends on the activity of chemokines that form concentration gradients to guide leukocytes to a specific site. Interaction of chemokines with their specific G protein-coupled receptors (GPCRs) on leukocytes induces leukocyte adhesion to the endothelial cells, followed by extravasation of the leukocytes and subsequent directed migration along the chemotactic gradient. Interaction of chemokines with glycosaminoglycans (GAGs) is crucial for extravasation in vivo. Chemokines need to interact with GAGs on endothelial cells and in the extracellular matrix in tissues in order to be presented on the endothelium of blood vessels and to create a concentration gradient. Local chemokine retention establishes a chemokine gradient and prevents diffusion and degradation. During the last two decades, research aiming at reducing chemokine activity mainly focused on the identification of inhibitors of the interaction between chemokines and their cognate GPCRs. This approach only resulted in limited success. However, an alternative strategy, targeting chemokine-GAG interactions, may be a promising approach to inhibit chemokine activity and inflammation. On this line, proteins derived from viruses and parasites that bind chemokines or GAGs may have the potential to interfere with chemokine-GAG interactions. Alternatively, chemokine mimetics, including truncated chemokines and mutant chemokines, can compete with chemokines for binding to GAGs. Such truncated or mutated chemokines are characterized by a strong binding affinity for GAGs and abrogated binding to their chemokine receptors. Finally, Spiegelmers that mask the GAG-binding site on chemokines, thereby preventing chemokine-GAG interactions, were developed. In this review, the importance of GAGs for chemokine activity in vivo and strategies that could be employed to target chemokine-GAG interactions will be discussed in the context of inflammation.

Syntheses of defined sulfated oligohyaluronans reveal structural effects, diversity and thermodynamics of GAG-protein binding

High binding affinities of GAG toward extracellular regulatory proteins are governed by recognition diversity, sulfation pattern, length, and anomeric functionalization.

Binding of sulfated glycosaminoglycans (GAG) to a wide spectrum of extracellular regulatory proteins is crucial for physiological processes such as cell growth, migration, tissue homeostasis and repair. Thus, GAG derivatives exhibit great relevance in the development of innovative biomaterials for tissue regeneration therapies. We present a synthetic strategy for the preparation of libraries of defined sulfated oligohyaluronans as model GAG systematically varied in length, sulfation pattern and anomeric substitution in order to elucidate the effects of these parameters on GAG recognition by regulatory proteins. Through an experimental and computational approach using fluorescence polarization, ITC, docking and molecular dynamics simulations we investigate the binding of these functionalized GAG derivatives to ten representative regulatory proteins including IL-8, IL-10, BMP-2, sclerostin, TIMP-3, CXCL-12, TGF-β, FGF-1, FGF-2, and AT-III, and we establish structure–activity relationships for GAG recognition. Binding is mainly driven by enthalpy with only minor entropic contributions. In several cases binding is determined by GAG length, and in all cases by the position and number of sulfates. Affinities strongly depend on the anomeric modification of the GAG. Highest binding affinities are effected by anomeric functionalization with large fluorophores and by GAG dimerization. Our experimental and theoretical results suggest that the diversity of GAG binding sites and modes is responsible for the observed high affinities and other binding features. The presented new insights into GAG–protein recognition will be of relevance to guide the design of GAG derivatives with customized functions for the engineering of new biomaterials.

Structural basis, stoichiometry, and thermodynamics of binding of the chemokines KC and MIP2 to the glycosaminoglycan heparin

Keratinocyte-derived chemokine (KC or mCXCL1) and macrophage inflammatory protein 2 (MIP2 or mCXCL2) play nonredundant roles in trafficking blood neutrophils to sites of infection and injury. The functional responses of KC and MIP2 are intimately coupled to their interactions with glycosaminoglycans (GAGs). GAG interactions orchestrate chemokine concentration gradients and modulate receptor activity, which together regulate neutrophil trafficking. Here, using NMR, molecular dynamics (MD) simulations, and isothermal titration calorimetry (ITC), we characterized the molecular basis of KC and MIP2 binding to the GAG heparin. Both chemokines reversibly exist as monomers and dimers, and the NMR analysis indicates that the dimer binds heparin with higher affinity. The ITC experiments indicate a stoichiometry of two GAGs per KC or MIP2 dimer and that the enthalpic and entropic contributions vary significantly between the two chemokine-heparin complexes. NMR-based structural models of heparin-KC and heparin-MIP2 complexes reveal that different combinations of residues from the N-loop, 40s turn, β₃-strand, and C-terminal helix form a binding surface within a monomer and that both conserved residues and residues unique to a particular chemokine mediate the binding interactions. MD simulations indicate significant residue-specific differences in their contribution to binding and affinity for a given chemokine and between chemokines. On the basis of our observations that KC and MIP2 bind to GAG via distinct molecular interactions, we propose that the differences in these GAG interactions lead to differences in neutrophil recruitment and play nonoverlapping roles in resolution of inflammation.