Investigation of Carbohydrate Recognition via Computer Simulation

Screening of Poria cocos polysaccharide with immunomodulatory activity and its activation effects on TLR4/MD2/NF-κB pathway

PCP-W1, the Poria cocos polysaccharide with the strong immunomodulatory activity, was isolated through column chromatography and screened for in vitro immune activity in RAW 264.7 cells in this study. The structure analysis results revealed that the PCP-W1 were composed of galactose, glucose, fucose and mannose in a molar percentage of 35.87: 28.56: 21.77: 13.64. And it exhibited a random coil and branched conformational features with a molecular weight of 18.38 kDa. The main chain consisted of residues→3)-β-D-Glcp-(1 → 3,6)-β-D-Glcp-(1 → 3)-β-D-Glcp-(1 → 6)-β-D-Glcp-(1 → 6)-α-D-Galp-(1 → 6)-α-D-Galp-(1 → 2,6)-α-D-Galp-(1→6)-α-D-Galp-(1 → 6)-α-D-Galp-(1 → , while branching occurred at β-D-Glcp-(1→, α-D-Manp-(1→, and α-L-Fucp-(1 → 3)- α-L-Fucp-(1→. The pharmacodynamic studies demonstrated that PCP-W1 activated the release of NO, IL-6, IL-β, TNF-α, CD86, and ROS to induce polarization of RAW 264.7 murine macrophages towards M1-type through modulation of the TLR4/MD2/NF-κB pathway. The molecular docking results showed that PCP-W1 could primarily dock onto the hydrophobic binding site of TLR4/MD2 complex via its galactose chain. Furthermore, molecular dynamics simulation displayed stable modeling for TLR4-MD2-PCP-W1 complex. Overall, we screened the most immunoactive components of the polysaccharide, analyzed its structure, demonstrated its impact on TLR4/MD2/NF-kB pathway, and studied the interaction between TLR4/MD2 and the polysaccharide fragments. These results provide further support for the structure-activity relationship study of the immunomodulatory effects of Poria cocos polysaccharide.

Sulfate glycosaminoglycan from swim bladder exerts immunomodulatory potential on macrophages via toll-like receptor 4 mediated NF-κB signaling pathways

This study explored the immunomodulatory impact and potential mechanisms on macrophages RAW264.7 using a purified macromolecular sulfate glycosaminoglycan (SBSG) from the swim bladder, whose structure was similar to chondroitin sulfate A. The results showed that SBSG at 0.25-1 mg/mL increased the viability and phagocytosis of RAW264.7 cells. Meanwhile, SBSG promoted the secretion of tumor necrosis factor α (TNF-α), interleukin 10 (IL-10), and nitric oxide (NO), as well as the production of reactive oxygen species (ROS). According to the RT-PCR and Western blot data, SBSG activated TLR4-nuclear factor kappa B (NF-κB) signaling pathways, which decreased the relative mRNA and protein levels of Toll-like receptor 4 (TLR4), IκB kinase β (IKKβ), NF-κB p65, and p-NF-κB p65. The molecular docking and molecular dynamic simulation findings revealed that the main binding force between TLR4 and SBSG was conventional hydrogen bond interaction, resulting in more stable ligand receptor complexes. In summary, SBSG exhibits significant immunomodulatory potential, similar to chondroitin sulfate C. The underlying molecular mechanism involved the binding of SBSG through hydrogen bonding to TLR4 receptors, triggering the NF-κB signaling pathway to downregulate the expression of related genes and proteins. This, in turn, regulated the secretion of various cytokines that were mediated by macrophages to exert the immunity of the body.

Lectin Histochemistry: Historical Perspectives, State of the Art, and Future Directions

Lectins, discovered more than 100 years ago and defined by their ability to selectively recognize specific carbohydrate structures, are ubiquitous in living organisms. Their precise functions are as yet under-explored and incompletely understood but they are clearly involved, through recognition of their binding partners, in a myriad of biological mechanisms involved in cell identity, adhesion, signaling, and growth regulation in health and disease. Understanding the complex "sugar code" represented by the "glycome" is a major challenge and at the forefront of current biological research. Lectins have been widely employed in histochemical studies to map glycosylation in cells and tissues. Here, a brief history of the discovery of lectins and early developments in their use is presented along with a selection of some of the most interesting and significant discoveries to emerge from the use of lectin histochemistry. Further, an evaluation of the next generation of lectin-based technologies is presented, including the potential for designing recombinant lectins with more precisely defined binding characteristics, linking lectin-based studies with other technologies to answer fundamental questions in glycobiology and approaches to exploring the interactions of lectins with their binding partners in more detail.

One-Bond 13C-1H and 13C-13C Spin-Coupling Constants as Constraints in MA'AT Analysis of Saccharide Conformation

MA'AT analysis uses ensembles of redundant experimental NMR spin-coupling constants, parametrized J-coupling equations obtained from density functional theory (DFT) calculations, and circular statistics to produce probability distributions of molecular torsion angles in solution and information on librational motions about these angles (Meredith et al., J. Chem. Info. Model. 2022, 62, 3135-3141). Current DFT methods give nearly quantitative two- and three-bond J_HH, J_CH, and ¹J_CC values for use in MA'AT analysis of saccharides. In contrast, the accuracy of DFT-calculated one-bond ¹J_CH and ¹J_CC values is more difficult to determine, preventing their use in MA'AT modeling. This report describes experimental and computational studies that address this problem using two approaches (Strategies 1 and 2). Differences [¹J_CH^calc - ¹J_CH^exp] (Strategy 1) ranged from -1.2 to 2.5 Hz, giving an average difference of 0.8 ± 1.7 Hz. Percent differences ranged from -0.8% to 1.6%, giving an average % difference of 0.5 ± 1.1%. In comparison, [¹J_CH^MA'AT - ¹J_CH^exp] (Strategy 2) ranged from -1.8 to 0.2 Hz, giving an average difference of -1.2 ± 0.7 Hz. Percent differences ranged from -1.2% to 0.1%, giving an average % difference of -0.8 ± 0.5%. Strategy 1 gave an average difference of 2.1 Hz between calculated and experimental ¹J_CC values, with an average % difference of 5.1 ± 0.2%. Calculated ¹J_CC values were consistently larger than experimental values. Strategy 2 also gave calculated ¹J_CC values that were larger than the experimental values, with an average difference of 2.3 ± 0.6 Hz, and an average % difference of 5.6 ± 1.6%. The findings of both strategies are similar and indicate that ¹J_CH values in saccharides can be calculated nearly quantitatively, but ¹J_CC values appear to be consistently overestimated by ∼5% using current DFT methods.

Nonconventional NMR Spin-Coupling Constants in Oligosaccharide Conformational Modeling: Structural Dependencies Determined from Density Functional Theory Calculations

Nonconventional NMR spin-coupling constants were investigated to determine their potential as conformational constraints in MA'AT modeling of the O-glycosidic linkages of oligosaccharides. Four (¹ J _C1',H1', ¹ J _C1',C2', ² J _C1',H2', and ² J _C2',H1') and eight (¹ J _C4,H4, ¹ J _C3,C4, ¹ J _C4,C5, ² J _C3,H4, ² J _C4,H3, ² J _C5,H4, ² J _C4,H5, and ² J _C3,C5) spin-couplings in methyl β-d-galactopyranosyl-(1→4)-β-d-glucopyranoside (methyl β-lactoside) were calculated using density functional theory (DFT) to determine their dependencies on O-glycosidic linkage C-O torsion angles, ϕ and ψ, respectively. Long-range ⁴ J _H1',H4 was also examined as a potential conformational constraint of either ϕ or ψ. Secondary effects of exocyclic (hydroxyl) C-O bond rotation within or proximal to these coupling pathways were investigated. Based on the findings of methyl β-lactoside, analogous J-couplings were studied in five additional two-bond O-glycosidic linkages [βGlcNAc-(1→4)-βMan, 2-deoxy-βGlc-(1→4)-βGlc, αMan-(1→3)-βMan, αMan-(1→2)-αMan, and βGlcNAc(1→2)-αMan] to determine whether the coupling behaviors observed in methyl β-lactoside were more broadly observed. Of the 13 nonconventional J-couplings studied, 7 exhibit properties that may be useful in future MA'AT modeling of O-glycosidic linkages, none of which involve coupling pathways that include the linkage C-O bonds. The findings also provide new insights into the general effects of exocyclic C-O bond conformation on the magnitude of experimental spin-couplings in saccharides and other hydroxyl-containing molecules.

Molecular dynamic simulation: Structural insights of multi-stranded curdlan in aqueous solution

In this work, by using molecular dynamic simulation we provide microscale structure information which helps to reveal the molecular mechanisms concerning the multi-chain conformational behavior of short curdlan. Through simulations starting with different conformations of curldan dodecasaccharides, it is found that the right-handed triple helix is thermodynamically the most stable conformation in aqueous solutions, which is well maintained and stabilized by an inter-strand hydrogen bonding network of the C₂ hydroxyls. Unlike any predicted forms, the inter-strand hydrogen bonds exhibit a left-handed double helix pattern with preferred global orientations. Temperature REMD results suggest that the formation of triple helix is temperature sensitive, but the already formed triple helix is not. Investigation of curdlan with numbers of repeating units from 3 to 12 captures a critical value of 6, which in a way elucidates the relationship between the formation of triple helix and the chain length.

The barley lectin, horcolin, binds high-mannose glycans in a multivalent fashion, enabling high-affinity, specific inhibition of cellular HIV infection

N-Linked glycans are critical to the infection cycle of HIV, and most neutralizing antibodies target the high-mannose glycans found on the surface envelope glycoprotein-120 (gp120). Carbohydrate-binding proteins, particularly mannose-binding lectins, have also been shown to bind these glycans. Despite their therapeutic potency, their ability to cause lymphocyte proliferation limits their application. In this study, we report one such lectin named horcolin (Hordeum vulgare lectin), seen to lack mitogenicity owing to the divergence in the residues at its carbohydrate-binding sites, which makes it a promising candidate for exploration as an anti-HIV agent. Extensive isothermal titration calorimetry experiments reveal that the lectin was sensitive to the length and branching of mannooligosaccharides and thereby the total valency. Modeling and simulation studies demonstrate two distinct modes of binding, a monovalent binding to shorter saccharides and a bivalent mode for higher glycans, involving simultaneous interactions of multiple glycan arms with the primary carbohydrate-binding sites. This multivalent mode of binding was further strengthened by interactions of core mannosyl residues with a secondary conserved site on the protein, leading to an exponential increase in affinity. Finally, we confirmed the interaction of horcolin with recombinant gp120 and gp140 with high affinity and inhibition of HIV infection at nanomolar concentrations without mitogenicity.

Predominant Nonproductive Substrate Binding by Fungal Cellobiohydrolase I and Implications for Activity Improvement

Enzymatic conversion of the most abundant renewable source of organic compounds, cellulose to fermentable sugars is attractive for production of green fuels and chemicals. The major component of industrial enzyme systems, cellobiohydrolase I from Hypocrea jecorina (Trichoderma reesei) (HjCel7A) processively splits disaccharide units from the reducing ends of tightly packed cellulose chains. HjCel7A consists of a catalytic domain (CD) and a carbohydrate-binding module (CBM) separated by a linker peptide. A tunnel-shaped substrate-binding site in the CD includes nine subsites for β-d-glucose units, seven of which (-7 to -1) precede the catalytic center. Low catalytic activity of Cel7A is the bottleneck and the primary target for improvement. Here it is shown for the first time that, in spite of much lower apparent k_cat of HjCel7A at the hydrolysis of β-1,4-glucosidic linkages in the fluorogenic cellotetra- and -pentaose compared to the structurally related endoglucanase I (HjCel7B), the specificity constants (catalytic efficiency) k_cat /K_m for both enzymes are almost equal in these reactions. The observed activity difference appears from strong nonproductive substrate binding by HjCel7A, particularly significant for MU-β-cellotetraose (MUG₄ ). Interaction of substrates with the subsites -6 and -5 proximal to the nonconserved Gln101 residue in HjCel7A decreases K_m,ap by >1500 times. HjCel7A can be nonproductively bound onto cellulose surface with K_d ≈2-9 nM via CBM and CD that captures six terminal glucose units of cellulose chain. Decomposition of this nonproductive complex can determine the rate of cellulose conversion. MUG₄ is a promising substrate to select active cellobiohydrolase I variants with reduced nonproductive substrate binding.

Characterizing protein conformations by correlation analysis of coarse-grained contact matrices

We have developed a method to capture the essential conformational dynamics of folded biopolymers using statistical analysis of coarse-grained segment-segment contacts. Previously, the residue-residue contact analysis of simulation trajectories was successfully applied to the detection of conformational switching motions in biomolecular complexes. However, the application to large protein systems (larger than 1000 amino acid residues) is challenging using the description of residue contacts. Also, the residue-based method cannot be used to compare proteins with different sequences. To expand the scope of the method, we have tested several coarse-graining schemes that group a collection of consecutive residues into a segment. The definition of these segments may be derived from structural and sequence information, while the interaction strength of the coarse-grained segment-segment contacts is a function of the residue-residue contacts. We then perform covariance calculations on these coarse-grained contact matrices. We monitored how well the principal components of the contact matrices is preserved using various rendering functions. The new method was demonstrated to assist the reduction of the degrees of freedom for describing the conformation space, and it potentially allows for the analysis of a system that is approximately tenfold larger compared with the corresponding residue contact-based method. This method can also render a family of similar proteins into the same conformational space, and thus can be used to compare the structures of proteins with different sequences.

Magnetic beads modified with an electron-transfer carbohydrate-mimetic peptide for sensing of a galactose-dependent protein

For use in the voltammetric sensing of galactose-dependent proteins, we modified magnetic beads with a peptide that had both electroactive- and molecular recognition properties. The peptide consisted of a YXY sequence and behaved as an electron-transfer carbohydrate-mimetic peptide that would combine with proteins. With this tool, the protein could be detected via a label-free system. We synthesized several penta- and hexa-peptides with a cysteine residue on the C-terminals to examine the properties of peptides. These peptides contained amino acid residues (X) of alanine, serine, or tyrosine. The peptides were immobilized on magnetic beads via N-(8-maleimidocapryloxy) succinimide. Soybean agglutinin(SBA), the in vivo function of which has been well established in animals, was selected as a model protein. The protein was detected via the changes in electrode response due to the oxidation of tyrosine residues from the phenol group to quinone. As a result, SBA was selectively accumulated on the beads modified with YYYYC. The calibration curve of SBA was linear and ranged from 2.5 × 10-12 to 1.0 × 10-10 M. With this system, SBA was recovered in human serum at values that ranged from 98 to 103%. Furthermore, the beads with peptides were regenerated five times using a protein denaturant. Accordingly, this electrochemical system was simple and could be rapidly applied to the detection of galactose-recognition proteins.

CAMERRA: An analysis tool for the computation of conformational dynamics by evaluating residue-residue associations

A computational method which extracts the dominant motions from an ensemble of biomolecular conformations via a correlation analysis of residue-residue contacts is presented. The algorithm first renders the structural information into contact matrices, then constructs the collective modes based on the correlated dynamics of a selected set of dynamic contacts. Associated programs can bridge the results for further visualization using graphics software. The aim of this method is to provide an analysis of conformations of biopolymers from the contact viewpoint. It may assist a systematical uncovering of conformational switching mechanisms existing in proteins and biopolymer systems in general by statistical analysis of simulation snapshots. In contrast to conventional correlation analyses of Cartesian coordinates (such as distance covariance analysis and Cartesian principal component analysis), this program also provides an alternative way to locate essential collective motions in general. Herein, we detail the algorithm in a stepwise manner and comment on the importance of the method as applied to decoding allosteric mechanisms. © 2018 Wiley Periodicals, Inc.

Rational Design of an Anticalin-Type Sugar-Binding Protein Using a Genetically Encoded Boronate Side Chain

The molecular recognition of carbohydrates plays a fundamental role in many biological processes. However, the development of carbohydrate-binding reagents for biomedical research and use poses a challenge due to the generally poor affinity of proteins toward sugars in aqueous solution. Here, we describe the effective molecular recognition of pyranose monosaccharides (in particular, galactose and mannose) by a rationally designed protein receptor based on the human lipocalin scaffold (Anticalin). Complexation relies on reversible covalent cis-diol boronate diester formation with a genetically encoded l-boronophenylalanine (Bpa) residue which was incorporated as a non-natural amino acid at a sterically permissive position in the ligand pocket of the Anticalin, as confirmed by X-ray crystallography. Compared with the metal-ion and/or avidity-dependent oligovalent lectins that prevail in nature, our approach offers a novel and promising route to generate tight sugar-binding reagents both as research reagents and for biomedical applications.

Lectin Histochemistry: Historical Perspectives, State of the Art, and the Future

Lectins, discovered more than 100 years ago and defined by their ability to selectively recognize specific carbohydrate structures, are ubiquitous in living organisms. Their precise functions are as yet under-explored and incompletely understood but they are clearly involved, through recognition of their binding partners, in a myriad of biological mechanisms involved in cell identity, adhesion, signaling, growth regulation, in health and disease. Understanding the complex "sugar code" represented by the glycome is a major challenge and at the forefront of current biological research. Lectins have been widely employed in histochemical studies to map glycosylation in cells and tissues. Here, a brief history of the discovery of lectins and early developments in their use is presented along with a selection of some of the most interesting and significant discoveries to emerge from use of lectin histochemistry. Further, an evaluation of the next generation of lectin-based technologies is presented, including the potential for designing recombinant lectins with more precisely defined binding characteristics, linking lectin-based studies with other technologies to answer fundamental questions in glycobiology, and approaches to exploring the interactions of lectins with their binding partners in more detail.

Probing the role of cation-π interaction in the thermotolerance and catalytic performance of endo-polygalacturonases

Understanding the dynamics of the key pectinase, polygalacturonase, and improving its thermotolerance and catalytic efficiency are of importance for the cost-competitive bioconversion of pectic materials. By combining structure analysis and molecular dynamics (MD) simulations, eight mutagenesis sites having the potential to form cation-π interactions were identified in the widely used fungal endo-polygalacturonase PG63. In comparison to the wild-type, three single mutants H58Y, T71Y and T304Y showed improved thermostability (the apparent Tms increased by 0.6-3.9 °C) and catalytic efficiency (by up to 32-fold). Chromatogram analysis of the hydrolysis products indicated that a larger amount of shorter sugars were released from the polygalacturonic acid by these three mutants than by the wild-type. MD analysis of the enzyme-substrate complexes illustrated that the mutants with introduced cation-π interaction have modified conformations of catalytic crevice, which provide an enviable environment for the catalytic process. Moreover, the lower plasticity of T3 loop 2 at the edge of the subsite tunnel appears to recruit the reducing ends of oligogalacturonide into the active site tunnel and initiates new hydrolysis reactions. This study demonstrates the importance of cation-π interaction in protein conformation and provides a realistic strategy to enhance the thermotolerance and catalytic performance of endo-polygalacturonases.

Force Fields for Carbohydrate-Divalent Cation Interactions

We report molecular dynamics simulations to study intermolecular interactions for carbohydrate-divalent cation complexes. We observed that common force fields from literature with standard Lorentz-Berthelot combining rules are unable to reproduce the experimental stability constants for model carbohydrate monomer (α-d-Allopyranose) and alkali earth metal cation (Mg(2+), Ca(2+), Sr(2+), or Ba(2+)) complexes. A modified combining rule with rescaled effective cross-interaction radius between cations and the hydroxyl oxygens on the carbohydrates was introduced to reproduce the experimental stability constants, which the preferential carbohydrate-cation complexing structures through the ax-eq-ax sequence of O-1, O-2, and O-3 on α-d-Allopyranose were also observed. The effective radius scaling factor obtained from (α-d-Allopyranose)-Ca(2+) complexes was directly transferrable to the similar six-membered ring (α-d-Ribopyranose)-Ca(2+) complexes; however, reparameterization for the scaling factor may be necessary for the five-membered ring (α-d-Ribofuranose)-Ca(2+) complexes.

The Promiscuity of Allosteric Regulation of Nuclear Receptors by Retinoid X Receptor

The promiscuous protein retinoid X receptor (RXR) displays essential allosteric regulation of several members in the nuclear hormone receptor superfamily via heterodimerization and (anti)cooperative binding of cognate ligands. Here, the structural basis of the positive allostery of RXR and constitutive androstane receptor (CAR) is revealed. In contrast, a similar computational approach had previously revealed the mechanism for negative allostery in the complex of RXR and thyroid receptor (TR). By comparing the positive and negative allostery of RXR complexed with CAR and TR respectively, we reported the promiscuous allosteric control involving RXR. We characterize the allosteric mechanism by expressing the correlated dynamics of selected residue-residue contacts which was extracted from atomistic molecular dynamics simulation and statistical analysis. While the same set of residues in the binding pocket of RXR may initiate the residue-residue interaction network, RXR uses largely different sets of contacts (only about one-third identical) and allosteric modes to regulate TR and CAR. The promiscuity of RXR control may originate from multiple factors, including (1) the frustrated fit of cognate ligand 9c to the RXR binding pocket and (2) the different ligand-binding features of TR (loose) versus CAR (tight) to their corresponding cognate ligands.

Protein-carbohydrate interactions as part of plant defense and animal immunity

The immune system consists of a complex network of cells and molecules that interact with each other to initiate the host defense system. Many of these interactions involve specific carbohydrate structures and proteins that specifically recognize and bind them, in particular lectins. It is well established that lectin-carbohydrate interactions play a major role in the immune system, in that they mediate and regulate several interactions that are part of the immune response. Despite obvious differences between the immune system in animals and plants, there are also striking similarities. In both cases, lectins can play a role as pattern recognition receptors, recognizing the pathogens and initiating the stress response. Although plants do not possess an adaptive immune system, they are able to imprint a stress memory, a mechanism in which lectins can be involved. This review will focus on the role of lectins in the immune system of animals and plants.