GLYCAM06: a generalizable biomolecular force field. Carbohydrates

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.

Improved physicochemical properties of bigels produced with ethyl cellulose-based oleogel and moderately deacetylated konjac glucomannan hydrogel

The ideal physicochemical properties of bigels are important for food applications. Therefore, a new bigel was prepared based on mixed beef tallow and soybean oil oleogel and deacetylated konjac glucomannan (KGM) hydrogel. The effect of the deacetylation degree of KGM on the physicochemical properties and microstructure of bigels was studied. The bigel containing moderate deacetylation degree of KGM had better rheological properties and hardness (319.84 g) than that with low and high deacetylation degrees of KGM. The interactions among the bigel components were analyzed by Fourier transform infrared spectroscopy and molecular dynamics simulation, indicating that the formation of the bigels was dominated by electrostatic interactions. Overall, the bigels containing moderate deacetylation degree of KGM had better physical properties, which may provide a theoretical foundation to develop bigels with low cholesterol, trans and saturated fats levels to replace traditional solid fats in food industry.

Antigenic determinants driving serogroup-specific antibody response to Neisseria meningitidis C, W, and Y capsular polysaccharides: Insights for rational vaccine design

Meningococcal glycoconjugate vaccines sourced from capsular polysaccharides (CPSs) of pathogenic Neisseria meningitidis strains are well-established measures to prevent meningococcal disease. However, the exact structural factors responsible for antibody recognition are not known. CPSs of Neisseria meningitidis serogroups Y and W differ by a single stereochemical center, yet they evoke specific immune responses. Herein, we developed specific monoclonal antibodies (mAbs) targeting serogroups C, Y, and W and evaluated their ability to kill bacteria. We then used these mAbs to dissect structural elements responsible for carbohydrate-protein interactions. First, Men oligosaccharides were screened against the mAbs using ELISA to select putative lengths representing the minimal antigenic determinant. Next, molecular interaction features between the mAbs and serogroup-specific sugar fragments were elucidated using STD-NMR. Moreover, X-ray diffraction data with the anti-MenW CPS mAb enabled the elucidation of the sugar-antibody binding mode. Our findings revealed common traits in the epitopes of all three sialylated serogroups. The minimal binding epitopes typically comprise five to six repeating units. Moreover, the O-acetylation of the neuraminic acid moieties was fundamental for mAb binding. These insights hold promise for the rational design of optimized meningococcal oligosaccharides, opening new avenues for novel production methods, including chemical or enzymatic approaches.

Chemical Resolution of an Epitope Recognized by Blood Group Antibodies Capable of Binding Both A and B Red Blood Cells

The high specificity of human antibodies to blood group A and B antigens is impressive, especially when considering the structural difference between these antigens (tetrasaccharides) is a NHAc versus a hydroxyl group on the terminal monosaccharide residue. It is well established that in addition to anti-A and anti-B there is a third antibody, anti-A,B capable of recognizing both A and B antigens. To analyze this AB specificity, we synthesized a tetrasaccharide, where the NHAc of the A antigen was replaced with an NH2. This NH2 group was then used to attach the glycan to an affinity resin, creating an AB epitope (ABep) adsorbent where the critical site for recognition by A and B antibodies was not accessible, while the rest of the (conformationally compact) tetrasaccharide remained accessible. Anti-ABep antibodies were then isolated from blood group O donors and found to have expected A,B specificity against immobilized and red cell bound synthetic antigens, including ABep, and were able to agglutinate both A and B red cells. The amount of these anti-ABep (anti-A,B) antibodies found in the blood of group O donors was comparable to levels of anti-A and anti-B found in group B and A individuals. Using STD-NMR the location for the AB epitope on the tetrasaccharide was found.

In silico approaches for better understanding cysteine cathepsin-glycosaminoglycan interactions

Cysteine cathepsins constitute the largest cathepsin family, with 11 proteases in human that are present primarily within acidic endosomal and lysosomal compartments. They are involved in the turnover of intracellular and extracellular proteins. They are synthesized as inactive procathepsins that are converted to mature active forms. Cathepsins play important roles in physiological and pathological processes and, therefore, receive increasing attention as potential therapeutic targets. Their maturation and activity can be regulated by glycosaminoglycans (GAGs), long linear negatively charged polysaccharides composed of recurring dimeric units. In this review, we summarize recent computational progress in the field of (pro)cathepsin-GAG complexes analyses.

Computational toolbox for the analysis of protein-glycan interactions

Protein-glycan interactions play pivotal roles in numerous biological processes, ranging from cellular recognition to immune response modulation. Understanding the intricate details of these interactions is crucial for deciphering the molecular mechanisms underlying various physiological and pathological conditions. Computational techniques have emerged as powerful tools that can help in drawing, building and visualising complex biomolecules and provide insights into their dynamic behaviour at atomic and molecular levels. This review provides an overview of the main computational tools useful for studying biomolecular systems, particularly glycans, both in free state and in complex with proteins, also with reference to the principles, methodologies, and applications of all-atom molecular dynamics simulations. Herein, we focused on the programs that are generally employed for preparing protein and glycan input files to execute molecular dynamics simulations and analyse the corresponding results. The presented computational toolbox represents a valuable resource for researchers studying protein-glycan interactions and incorporates advanced computational methods for building, visualising and predicting protein/glycan structures, modelling protein-ligand complexes, and analyse MD outcomes. Moreover, selected case studies have been reported to highlight the importance of computational tools in studying protein-glycan systems, revealing the capability of these tools to provide valuable insights into the binding kinetics, energetics, and structural determinants that govern specific molecular interactions.

Small lectin ligands as a basis for applications in glycoscience and glycomedicine

Glycan recognition by lectins mediates important biological events. This Tutorial Review aims to introduce lectin-ligand interactions and show how these molecular recognition events inspire innovations such as: (i) glycomimetic ligands; (ii) multivalent ligand agonists/antagonists; (iii) ligands for precision delivery of therapies to cells, where therapies include vaccines, siRNA and LYTACs (iv) development of diagnostics. A small number of case studies are selected to demonstrate principles for development of new ligands for applications inspired by knowledge of natural glycan ligand structure and function.

Molecular Dynamics Study on Adsorption and Desorption of the Model Oligosaccharide above Polymer Antifouling Membranes

Polysaccharide foulants play a key role in the adhesion of many fouling organisms, which may cause severe marine biofouling. However, the detailed interaction mechanism between polysaccharides and antifouling membranes is still indistinct compared with that between the fouling protein and antifouling surfaces. In this paper, a model oligosaccharide built based on the monosaccharide composition found in diatom extracellular polymeric substances (EPS) was used as a model foulant to investigate its adsorption and desorption above three T4 antifouling membranes. It was found that the anionic poly(3-(methacryloyloxy)propane-1-sulfonate) (T4-SP) antifouling membrane had excellent antifouling ability with respect to the model oligosaccharide, while the oligosaccharide can be easily adsorbed on the poly(2-(dimethylamino)ethyl methacrylate) (T4-DM) membrane with vdW attraction and on the zwitterionic poly(sulfobetaine methacrylate) (T4-SB) membrane with electrostatic attraction. As little is known about the details of polysaccharides' adsorption above antifouling membranes at the molecular level, we hope this work will serve as a theoretical basis for finding more effective materials to prevent or control marine biofouling.

Scorpion α-toxin LqhαIT specifically interacts with a glycan at the pore domain of voltage-gated sodium channels

Voltage-gated sodium (Nav) channels sense membrane potential and drive cellular electrical activity. The deathstalker scorpion α-toxin LqhαIT exerts a strong action potential prolonging effect on Nav channels. To elucidate the mechanism of action of LqhαIT, we determined a 3.9 Å cryoelectron microscopy (cryo-EM) structure of LqhαIT in complex with the Nav channel from Periplaneta americana (NavPas). We found that LqhαIT binds to voltage sensor domain 4 and traps it in an "S4 down" conformation. The functionally essential C-terminal epitope of LqhαIT forms an extensive interface with the glycan scaffold linked to Asn330 of NavPas that augments a small protein-protein interface between NavPas and LqhαIT. A combination of molecular dynamics simulations, structural comparisons, and prior mutagenesis experiments demonstrates the functional importance of this toxin-glycan interaction. These findings establish a structural basis for the specificity achieved by scorpion α-toxins and reveal the conserved glycan as an essential component of the toxin-binding epitope.

Alteration of Neuropilin-1 and Heparan Sulfate Interaction Impairs Murine B16 Tumor Growth

Neuropilin-1 acts as a coreceptor with vascular endothelial growth factor receptors to facilitate binding of its ligand, vascular endothelial growth factor. Neuropilin-1 also binds to heparan sulfate, but the functional significance of this interaction has not been established. A combinatorial library screening using heparin oligosaccharides followed by molecular dynamics simulations of a heparin tetradecasaccharide suggested a highly conserved binding site composed of amino acid residues extending across the b1 and b2 domains of murine neuropilin-1. Mutagenesis studies established the importance of arginine513 and lysine514 for binding of heparin to a recombinant form of Nrp1 composed of the a1, a2, b1, and b2 domains. Recombinant Nrp1 protein bearing R513A,K514A mutations showed a significant loss of heparin-binding, heparin-induced dimerization, and heparin-dependent thermal stabilization. Isothermal calorimetry experiments suggested a 1:2 complex of heparin tetradecasaccharide:Nrp1. To study the impact of altered heparin binding in vivo, a mutant allele of Nrp1 bearing the R513A,K514A mutations was created in mice (Nrp1D) and crossbred to Nrp1+/- mice to examine the impact of altered heparan sulfate binding. Analysis of tumor formation showed variable effects on tumor growth in Nrp1D/D mice, resulting in a frank reduction in tumor growth in Nrp1D/- mice. Expression of mutant Nrp1D protein was normal in tissues, suggesting that the reduction in tumor growth was due to the altered binding of heparin/heparan sulfate to neuropilin-1. These findings suggest that the interaction of neuropilin-1 with heparan sulfate modulates its stability and its role in tumor formation and growth.

Gut microbiota biotransformation of drug glucuronides leading to gastrointestinal toxicity: Therapeutic potential of bacterial β-glucuronidase inhibition in mycophenolate-induced enteropathy

AIMS

Drug-induced enteropathy is often associated with the therapeutic use of certain glucuronidated drugs. One such drug is mycophenolic acid (MPA), a well-established immunosuppressant of which gastrointestinal adverse effects are a major concern. The role of bacterial β-glucuronidase (β-G) from the gut microbiota in MPA-induced enteropathy has recently been discovered. Bacterial β-G hydrolyzes MPAG, the glucuronide metabolite of MPA excreted in the bile, leading to the digestive accumulation of MPA that would favor in turn these adverse events. We therefore hypothesized that taming bacterial β-G activity might reduce MPA digestive exposure and prevent its toxicity.

MAIN METHODS

By using a multiscale approach, we evaluated the effect of increasing concentrations of MPA on intestinal epithelial cells (Caco-2 cell line) viability, proliferation, and migration. Then, we investigated the inhibitory properties of amoxapine, a previously described bacterial β-G inhibitor, by using molecular dynamics simulations, and evaluated its efficiency in blocking MPAG hydrolysis in an Escherichia coli-based β-G activity assay. The pharmacological effect of amoxapine was evaluated in a mouse model.

KEY FINDINGS

We observed that MPA impairs intestinal epithelial cell homeostasis. Amoxapine efficiently blocks the hydrolysis of MPAG to MPA and significantly reduces digestive exposure to MPA in mice. As a result, administration of amoxapine in MPA-treated mice significantly attenuated gastrointestinal lesions.

SIGNIFICANCE

Collectively, these results suggest that the digestive accumulation of MPA is involved in the pathophysiology of MPA-gastrointestinal adverse effects. This study provides a proof-of-concept of the therapeutic potential of bacterial β-G inhibitors in glucuronidated drug-induced enteropathy.

Conformational preferences of heparan sulfate to recognize the CXCL8 dimer in aqueous medium: degree of sulfation and hydrogen bonds

The sulfation pattern and epimerization of the long-chain sulfated polysaccharide heparan sulfate (HS) cause structural diversity and regulate various physiological and pathological processes when binding with proteins. In this work, we performed a series of molecular dynamics simulations of three variants of the octadecasaccharide HS with varying sulfation positions in aqueous medium in their free forms and in the presence of the chemokine CXCL8 dimer. The free energy of binding depicts the sulfation at the 6-O position of GlcNAc (HS6S), and both 3-O and 6-O positions of GlcNAc (HS3S6S) of HS variants are more likely to bind with the CXCL8 dimer than the triply sulfated HS2S3S6S, which is sulfated at the 2-O position of GlcUA additionally along with 3-O and 6-O positions of GlcNAc. Binding between HS and CXCL8 was driven by electrostatic and van der Waals interactions predominantly regardless of the sulfation pattern; however, unfavorable entropic contribution suppressed the interaction between HS and CXCL8. The contribution of different amino acid residues to the binding energetics suggested that basic amino acids line up the binding site of CXCL8. This study further acknowledges the role of interfacial water that is structured and bound with HS through hydrogen bonds, exhibiting differential hydrogen bond relaxation dynamics compared to when the HS molecules are free. Moreover, this study identifies that with the increase in sulfation, the HS-water hydrogen bond relaxation occurs faster with the complexation, while the reverse trend is followed in their free forms. Significant structural adaptation of the different sulfated HS molecules, as verified from the free energy landscapes generated from various reaction coordinates, root-mean-square-deviations, end-to-end distances, including ring pucker angles, dihedral flexibility, and the high conformational entropy cost arising from the glycosidic bonds, suggests that the different sulfated variants of HS undergo significant structural transformation to bind with CXCL8. The presence of a CXCL8 dimer imposes the bound forms of HS to adopt non-linear structures with skew-boat conformations. The atomistic details of the study would help in understanding the selectivity and conformational diversity, as well as the role of solvents in the recognition of CXCL8 by different sulfated variants of HS molecules.

Machine-learned molecular mechanics force fields from large-scale quantum chemical data

The development of reliable and extensible molecular mechanics (MM) force fields-fast, empirical models characterizing the potential energy surface of molecular systems-is indispensable for biomolecular simulation and computer-aided drug design. Here, we introduce a generalized and extensible machine-learned MM force field, espaloma-0.3, and an end-to-end differentiable framework using graph neural networks to overcome the limitations of traditional rule-based methods. Trained in a single GPU-day to fit a large and diverse quantum chemical dataset of over 1.1 M energy and force calculations, espaloma-0.3 reproduces quantum chemical energetic properties of chemical domains highly relevant to drug discovery, including small molecules, peptides, and nucleic acids. Moreover, this force field maintains the quantum chemical energy-minimized geometries of small molecules and preserves the condensed phase properties of peptides and folded proteins, self-consistently parametrizing proteins and ligands to produce stable simulations leading to highly accurate predictions of binding free energies. This methodology demonstrates significant promise as a path forward for systematically building more accurate force fields that are easily extensible to new chemical domains of interest.

Computational insights into the mechanisms underlying structural destabilization and recovery in trafficking-deficient hERG mutants

Cardiovascular diseases are a major global health concern, responsible for a significant number of deaths each year, often linked to cardiac arrhythmias resulting from dysfunction in ion channels. Hereditary Long QT Syndrome (LQTS) is a condition characterized by a prolonged QT interval on ECG, increasing the risk of sudden cardiac death. The most common type of LQTS, LQT2, is caused by mutations in the hERG gene, affecting a potassium ion channel. The majority of these mutations disrupt the channel's trafficking to the cell membrane, leading to intracellular retention. Specific high-affinity hERG blockers (e.g., E-4031) can rescue this mutant phenotype, but the exact mechanism is unknown. This study used accelerated molecular dynamics simulations to investigate how these mutations affect the hERG channel's structure, folding, endoplasmic reticulum (ER) retention, and trafficking. We reveal that these mutations induce structural changes in the channel, narrowing its central pore and altering the conformation of the intracellular domains. These changes expose internalization signals that contribute to ER retention and degradation of the mutant hERG channels. Moreover, the study found that the trafficking rescue drug E-4031 can inhibit these structural changes, potentially rescuing the mutant channels. This research offers valuable insights into the structural issues responsible for the degradation of rescuable transmembrane trafficking mutants. Understanding the defective trafficking structure of the hERG channel could help identify binding sites for small molecules capable of restoring proper folding and facilitating channel trafficking. This knowledge has the potential to lead to mechanism-based therapies that address the condition at the cellular level, which may prove more effective than treating clinical symptoms, ultimately offering hope for individuals with hereditary Long QT Syndrome.

Dependence of lactose adsorption on the exposed crystal facets of metals: a comparative study of gold, silver and copper

In this theoretical work, we investigated the adsorption of a lactose molecule on metal-based surfaces, with a focus on the influence of the nature of the metal and of the type of exposed crystal facet on the adsorbed structures and energetics. More precisely, we considered three flat crystallographic facets of three face-centered cubic metals (gold, silver, and copper). For the global exploration of the energy landscape, we employed a multi-stage procedure where high-throughput searches, using a stochastic method that performs global optimization by iterating local searches, are followed by a refinement of the most probable adsorption conformations of the molecule at the ab initio level. We predicted the optimal conformation of lactose on each of the nine metal-surface combinations, classified the many low-energy minima into possible adsorption modes, and analyzed the structural, electronic and energetic aspects of the lactose molecule on the surface, as well as their dependence on the type of metal and exposed crystal facet. We observed structural similarities between the various minimum-energy conformations of lactose in vacuum and on the surface, a rough correlation between adsorption and interaction energies of the molecule, and a small charge transfer between molecule and surface whose direction is metal-dependent. During adsorption, an electronic reorganization occurs at the metal-molecule interface only, without affecting the vacuum-pointing atoms of the lactose molecule. For all types of surfaces, lactose exhibits the weakest adsorption on silver substrates, while for each coinage metal the adsorption is strongest on the (110) crystal facet. This study demonstrates that the control of exposed facets can allow to modulate the interaction between metals and small saccharides.

Improved moisture barrier and mechanical properties of rice protein/sodium alginate films for banana and oil preservation: Effect of the type and addition form of fatty acid

Attenuating the moisture sensitivity of hydrophilic protein/polysaccharide-based films without impairing other properties remains a challenge. Fatty acid dispersed in Pickering emulsion was proposed to overcome such issue. An increase in fatty acid chain length slightly reduced the water vapor permeability (WVP) of emulsion films. As the number of fatty acid double bonds increased from 0 to 1, the WVP of emulsion films was significantly decreased by 14.02% while mechanical properties were significantly enhanced. More hydrogen bonds and stronger electrostatic interactions in the presence of fatty acids were observed by molecular dynamics simulation. The weight loss of bananas coated with oleic acid-incorporated film-forming emulsion was 6.81% lower than that of uncoated group after 4 days, and the corresponding film was more effective to delay oil oxidation than the commercial polypropylene film, indicating that the film is a promising alternative to food coating and packaging material.

Unravelling biochemical and structural features of Bacillus licheniformis GH5 mannanase using site-directed mutagenesis and high-resolution protein crystallography studies

Glycoside hydrolase family 5 (GH5) encompasses enzymes with several different activities, including endo-1,4-β-mannosidases. These enzymes are involved in mannan degradation, and have a number of biotechnological applications, such as mannooligosaccharide prebiotics production, stain removal and dyes decolorization, to name a few. Despite the importance of GH5 enzymes, only a few members of subfamily 7 were structurally characterized. In the present work, biochemical and structural characterization of Bacillus licheniformis GH5 mannanase, BlMan5_7 were performed and the enzyme cleavage pattern was analyzed, showing that BlMan5_7 requires at least 5 occupied subsites to perform efficient hydrolysis. Additionally, crystallographic structure at 1.3 Å resolution was determined and mannoheptaose (M7) was docked into the active site to investigate the interactions between substrate and enzyme through molecular dynamic (MD) simulations, revealing the existence of a - 4 subsite, which might explain the generation of mannotetraose (M4) as an enzyme product. Biotechnological application of the enzyme in stain removal was investigated, demonstrating that BlMan5_7 addition to washing solution greatly improves mannan-based stain elimination.

Cellulose nanocrystal-stabilized Pickering emulsion gels as vehicles for follicular delivery of minoxidil

Minoxidil (MXD) is the only topical over-the-counter medicine approved by the United States Food and Drug Administration for the treatment of androgenetic alopecia. For the purpose of targeting the delivery of MXD to dermal papilla in the hair follicle, MXD Pickering emulsion gels were fabricated based on the designability of deep eutectic solvent (DES) and the versatility of cellulose nanocrystal (CNC) and sodium carboxymethyl cellulose (CMC-Na). Structural studies and theoretical calculations results suggest that CNC can stabilize the interface between the MXD-DES and water, leading to the formation of Pickering emulsions. The rheological properties and stabilities of MXD Pickering emulsions were enhanced through gelation using CMC-Na, which highlights the good compatibility and effectiveness of natural polysaccharides in emulsion gels. Due to the particle size of emulsion droplets (679 nm) and the rheological properties of emulsion gel, the fabricated MXD formulations show in vivo hair regrowth promotion and hair follicle targeting capabilities. Interestingly, the MXD Pickering emulsion-based formulations exert therapeutic effects by upregulating the expression of hair growth factors. The proposed nanodrug strategy based on supramolecular strategies of CNC and CMC-Na provides an interesting avenue for androgenetic alopecia treatment.

Soluble Human Lectins at the Host-Microbe Interface

Human lectins are integral to maintaining microbial homeostasis on the skin, in the blood, and at mucosal barriers. These proteins can recognize microbial glycans and inform the host about its microbial status. In accordance with their roles, their production can vary with tissue type. They also can have unique structural and biochemical properties, and they can influence microbial colonization at sites proximal and distal to their tissue of origin. In line with their classification as innate immune proteins, soluble lectins have long been studied in the context of acute infectious disease, but only recently have we begun to appreciate their roles in maintaining commensal microbial communities (i.e., the human microbiota). This review provides an overview of soluble lectins that operate at host-microbe interfaces, their glycan recognition properties, and their roles in physiological and pathological mechanisms.

The effect of host size on binding in host-guest complexes of cyclodextrins and polyoxometalates

Harnessing flexible host cavities opens opportunities for the design of novel supramolecular architectures that accommodate nanosized guests. This research examines unprecedented gas-phase structures of Keggin-type polyoxometalate PW12O40 3- (WPOM) and cyclodextrins (X-CD, X = α, β, γ, δ, ε, ζ) including previously unexplored large, flexible CDs. Using ion mobility spectrometry coupled to mass spectrometry (IM-MS) in conjunction with molecular dynamics (MD) simulations, we provide first insights into the binding modes between WPOM and larger CD hosts as isolated structures. Notably, γ-CD forms two distinct structures with WPOM through binding to its primary and secondary faces. We also demonstrate that ε-CD forms a deep inclusion complex, which encapsulates WPOM within its annular inner cavity. In contrast, ζ-CD adopts a saddle-like conformation in its complex with WPOM, which resembles its free form in solution. More intriguingly, the gas-phase CD-WPOM structures are highly correlated with their counterparts in solution as characterized by nuclear magnetic resonance (NMR) spectroscopy. The strong correlation between the gas- and solution phase structures of CD-WPOM complexes highlight the power of gas-phase IM-MS for the structural characterization of supramolecular complexes with nanosized guests, which may be difficult to examine using conventional approaches.

In vitro and in silico studies of the inclusion complexation of 8-bromobaicalein with β-cyclodextrins

Baicalein, a flavone derived from Scutellaria baicalensis Georgi, exhibits potent anti-inflammatory, antiviral, and anticancer properties. Its derivative, known as 8-bromobaicalein (BB), has been found to have strong cytotoxic effect on MCF-7 human breast cancer cells. However, its limited solubility in water has hindered its potential for wider applications. To address this issue, we investigated the use of cyclodextrins specifically βCD, 2,6-di-O-methyl-β-cyclodextrin (DMβCD), and hydroxypropyl-β-cyclodextrin (HPβCD) to improve the solubility of BB through inclusion complexation. During 250 ns molecular dynamics simulations, it was found that BB can form inclusion complexes with all βCDs. These complexes exhibit two distinct orientations: chromone group insertion (C-form) and phenyl group insertion (P-form). The formation of these complexes is primarily driven by van der Waals interactions. DMβCD has the highest number of atom contacts with BB and the lowest solvent accessibility in the hydrophobic cavity. These results coincide with the highest binding affinity from the MM/GBSA-based free energy calculation method. Experimental phase solubility diagrams revealed a 1:1 stoichiometric ratio (AL type) between BB and βCDs, in which BB/DMβCD showed the highest stability. The formation of inclusion complexes was confirmed by differential scanning calorimetry and scanning electron microscope methods. Additionally, the BB/DMβCD inclusion complex demonstrated significantly higher anticancer activity against MCF-7 human breast cancer cells compared to BB alone. These findings underscore the potential of DMβCD for formulating BB in pharmaceutical and medical applications.

O-glycosylation of SARS-CoV-2 spike protein by host O-glycosyltransferase strengthens its trimeric structure

Protein O-glycosylation, also known as mucin-type O-glycosylation, is one of the most abundant glycosylation in mammalian cells. It is initially catalyzed by a family of polypeptide GalNAc transferases (ppGalNAc-Ts). The trimeric spike protein (S) of SARS-CoV-2 is highly glycosylated and facilitates the virus's entry into host cells and membrane fusion of the virus. However, the functions and relationship between host ppGalNAc-Ts and O-glycosylation on the S protein remain unclear. Herein, we identify 15 O-glycosites and 10 distinct O-glycan structures on the S protein using an HCD-product-dependent triggered ETD mass spectrometric analysis. We observe that the isoenzyme T6 of ppGalNAc-Ts (ppGalNAc-T6) exhibits high O-glycosylation activity for the S protein, as demonstrated by an on-chip catalytic assay. Overexpression of ppGalNAc-T6 in HEK293 cells significantly enhances the O-glycosylation level of the S protein, not only by adding new O-glycosites but also by increasing O-glycan heterogeneity. Molecular dynamics simulations reveal that O-glycosylation on the protomer-interface regions, modified by ppGalNAc-T6, potentially stabilizes the trimeric S protein structure by establishing hydrogen bonds and non-polar interactions between adjacent protomers. Furthermore, mutation frequency analysis indicates that most O-glycosites of the S protein are conserved during the evolution of SARS-CoV-2 variants. Taken together, our finding demonstrate that host O-glycosyltransferases dynamically regulate the O-glycosylation of the S protein, which may influence the trimeric structural stability of the protein. This work provides structural insights into the functional role of specific host O-glycosyltransferases in regulating the O-glycosylation of viral envelope proteins.

Molecular Mechanism of Double-Displacement Retaining β-Kdo Glycosyltransferase WbbB

Glycosyltransferases (GTs) are pivotal enzymes involved in glycosidic bond synthesis, which can lead to either retention or inversion of the glycosyl moiety's anomeric configuration. However, the catalytic mechanism for retaining GTs remains a subject of controversy. In this study, we employ MD and QM/MM metadynamics to investigate the double-displacement catalytic mechanism of the retaining β-Kdo transferase WbbB. Our findings demonstrate that the nucleophile Asp232 initiates the reaction by attacking the sugar ring containing a carboxylate at the anomeric position, forming a covalent adduct. Subsequently, the adduct undergoes a rotational rearrangement, ensuring proper orientation of the anomeric carbon for the acceptor substrate. In the second step, Glu158 acts as the catalytic base to abstract the proton of the acceptor substrate to complete the transglycosylation reaction. Notably, His265 does not function as the anticipated catalytic acid; instead, it stabilizes the phosphate group through H-bonding interactions. Our simulations support the double-displacement mechanism implicated from the crystallographic studies of WbbB. This mechanism deviates from the common SNi-type and retaining glycoside hydrolase mechanisms, thereby expanding our understanding of GT catalytic mechanisms.

Deciphering the Role of Trehalose in Chroococcidiopsis sp. 029's High-Desiccation Resistance: Sequence Determination, Structural Modelling and Simulative Analysis of the 30S Ribosomal Subunit

Desert strains of the genus Chroococcidiopsis are among the most desiccation-resistant cyanobacteria capable of anhydrobiosis. The accumulation of two sugars, sucrose and trehalose, facilitates the entrance of anhydrobiotes into a reversible state of dormancy by stabilizing cellular components upon water removal. This study aimed to evaluate, at the atomistic level, the role of trehalose in desiccation resistance by using as a model system the 30S ribosomal subunit of the desert cyanobacterium Chroococcidiopsis sp. 029. Molecular dynamic simulations provided atomistic evidence regarding its protective role on the 30S molecular structure. Trehalose forms an enveloping shell around the ribosomal subunit and stabilizes the structures through a network of direct interactions. The simulation confirmed that trehalose actively interacts with the 30S ribosomal subunit and that, by replacing water molecules, it ensures ribosomal structural integrity during desiccation, thus enabling protein synthesis to be carried out upon rehydration.

The Conformation of Glycosidic Linkages According to Various Force Fields: Monte Carlo Modeling of Polysaccharides Based on Extrapolation of Short-Chain Properties

The conformational features of the glycosidic linkage are the most important variable to consider when studying di-, oligo-, and polysaccharide molecules using molecular dynamics (MD) simulations. The accuracy of the theoretical model describing this degree of freedom influences the quality of the results obtained from MD calculations based on this model. This article focuses on the following two issues related to the conformation of the glycosidic linkage. First, we describe the results of a comparative analysis of the predictions of three carbohydrate-dedicated classical force fields for MD simulations, namely, CHARMM, GLYCAM, and GROMOS, in the context of different parameters of structural and energetic nature related to the conformation of selected types of glycosidic linkages, α(1 → 4), β(1 → 3), and β(1 → 4), connecting glucopyranose units. This analysis revealed several differences, mainly concerning the energy levels of the secondary and tertiary conformers and the linkage flexibility within the dominant exo-syn conformation for α(1 → 4) and β(1 → 3) linkages. Some aspects of the comparative analysis also included the newly developed, carbohydrate-dedicated Martini 3 coarse-grained force field. Second, to overcome the time-scale problem associated with sampling slow degrees of freedom in polysaccharide chains during MD simulations, we developed a coarse-grained (CG) model based on the data from MD simulations and designed for Monte Carlo modeling. This model (CG MC) is based on information from simulations of short saccharide chains, effectively sampled in atomistic MD simulations, and is capable of extrapolating local conformational properties to the case of polysaccharides of arbitrary length. The CG MC model has the potential to estimate the conformations of very long polysaccharide chains, taking into account the influence of secondary and tertiary conformations of glycosidic linkages. With respect to the comparative analysis of force fields, the application of CG MC modeling showed that relatively small differences in the predictions of individual force fields with respect to a single glycosidic linkage accumulate when considering their effect on the structure of longer chains, leading to drastically different predictions with respect to parameters describing the polymer conformation, such as the persistence length.

Accumulation of Hydrogen Bonds and van der Waals Interactions Determines Force Response between Two Parallel Cellulose Chains: Steered Molecular Dynamics Simulations

We investigated the response forces between two parallel cellulose chains during the shearing and tearing processes by using steered molecular dynamics simulations. It was found that there are two logarithmic dependencies between response force and pulling speed in shearing processes but only one in tearing, according to Bell's equation by fitting the f-ln v curve. The mechanism is that there are 2-fold interactions determining the force response between two parallel cellulose chains resisting chain separation during a shearing process. Our results indicate that hydrogen bonds dominate the interchain interactions in the fast pull mode (FPM) for shearing, while van der Waals interactions dominate in the slow pull mode (SPM). For tearing, the one-by-one breaking of hydrogen bonds and van der Waals interactions plays a main role.

Three-Step Enzymatic Remodeling of Chitin into Bioactive Chitooligomers

Here we describe a complex enzymatic approach to the efficient transformation of abundant waste chitin, a byproduct of the food industry, into valuable chitooligomers with a degree of polymerization (DP) ranging from 6 to 11. This method involves a three-step process: initial hydrolysis of chitin using engineered variants of a novel fungal chitinase from Talaromyces flavus to generate low-DP chitooligomers, followed by an extension to the desired DP using the high-yielding Y445N variant of β-N-acetylhexosaminidase from Aspergillus oryzae, achieving yields of up to 57%. Subsequently, enzymatic deacetylation of chitooligomers with DP 6 and 7 was accomplished using peptidoglycan deacetylase from Bacillus subtilis BsPdaC. The innovative enzymatic procedure demonstrates a sustainable and feasible route for converting waste chitin into unavailable bioactive chitooligomers potentially applicable as natural pesticides in ecological and sustainable agriculture.

Descriptive molecular pharmacology of the δ opioid receptor (DOR): A computational study with structural approach

This work focuses on the δ receptor (DOR), a G protein-coupled receptor (GPCR) belonging to the opioid receptor group. DOR is expressed in numerous tissues, particularly within the nervous system. Our study explores computationally the receptor's interactions with various ligands, including opiates and opioid peptides. It elucidates how these interactions influence the δ receptor response, relevant in a wide range of health and pathological processes. Thus, our investigation aims to explore the significance of DOR as an incoming drug target for pain relief and neurodegenerative diseases and as a source for novel opioid non-narcotic analgesic alternatives. We analyze the receptor's structural properties and interactions using Molecular Dynamics (MD) simulations and Gaussian-accelerated MD across different functional states. To thoroughly assess the primary differences in the structural and conformational ensembles across our different simulated systems, we initiated our study with 1 μs of conventional Molecular Dynamics. The strategy was chosen to encompass the full activation cycle of GPCRs, as activation processes typically occur within this microsecond range. Following the cMD, we extended our study with an additional 100 ns of Gaussian accelerated Molecular Dynamics (GaMD) to enhance the sampling of conformational states. This simulation approach allowed us to capture a comprehensive range of dynamic interactions and conformational changes that are crucial for GPCR activation as influenced by different ligands. Our study includes comparing agonist and antagonist complexes to uncover the collective patterns of their functional states, regarding activation, blocking, and inactivation of DOR, starting from experimental data. In addition, we also explored interactions between agonist and antagonist molecules from opiate and opioid classifications to establish robust structure-activity relationships. These interactions have been systematically quantified using a Quantitative Structure-Activity Relationships (QSAR) model. This research significantly contributes to our understanding of this significant pharmacological target, which is emerging as an attractive subject for drug development.

O-glycans Expand Lubricin and Attenuate Its Viscosity and Shear Thinning

Lubricin, an intrinsically disordered glycoprotein, plays a pivotal role in facilitating smooth movement and ensuring the enduring functionality of synovial joints. The central domain of this protein serves as a source of this excellent lubrication and is characterized by its highly glycosylated, negatively charged, and disordered structure. However, the influence of O-glycans on the viscosity of lubricin remains unclear. In this study, we employ molecular dynamics simulations in the absence and presence of shear, along with continuum simulations, to elucidate the intricate interplay between O-glycans and lubricin and the impact of O-glycans on lubricin's conformational properties and viscosity. We found the presence of O-glycans to induce a more extended conformation in fragments of the disordered region of lubricin. These O-glycans contribute to a reduction in solution viscosity but at the same time weaken shear thinning at high shear rates, compared to nonglycosylated systems with the same density. This effect is attributed to the steric and electrostatic repulsion between the fragments, which prevents their conglomeration and structuring. Our computational study yields a mechanistic mechanism underlying previous experimental observations of lubricin and paves the way to a more rational understanding of its function in the synovial fluid.

Interactions of a Novel Anthracycline with Oligonucleotide DNA and Cyclodextrins in an Aqueous Environment

Berubicin, a chemotherapy medication belonging to the class of anthracyclines, is simulated in double-stranded DNA sequences and cyclodextrins in an aqueous environment via full-atom molecular dynamics simulations on the time scale of microseconds. The drug is studied in both the neutral and protonated states so as to better comprehend the role of its charge in the formed complexes. The noncovalent berubicin-DNA and berubicin-cyclodextrin complexes are investigated in detail, paying special attention to their thermodynamic description by employing the double decoupling method, the solvent balance method, the weighted solvent accessible surface model, and the linear interaction energy method. A novel approach for extracting the desolvation thermodynamics of the binding process is also presented. Both the binding and desolvation Gibbs energies are decomposed into entropic and enthalpic contributions so as to elucidate the nature of complexation and its driving forces. Selected structural and geometrical properties of all the complexes, which are all stable, are analyzed. Both cyclodextrins under consideration are widely utilized for drug delivery purposes, and a comparative investigation between their bound states with berubicin is carried out.

Methyl α-D-galactopyranosyl-(1→3)-β-D-galactopyranoside and methyl β-D-galactopyranosyl-(1→3)-β-D-galactopyranoside: Glycosidic linkage conformation determined from MA'AT analysis

MA'AT analysis has been applied to two biologically-important O-glycosidic linkages in two disaccharides, α-D-Galp-(1→3)-β-D-GalpOMe (3) and β-D-Galp-(1→3)-β-D-GalpOMe (4). Using density functional theory (DFT) to obtain parameterized equations relating a group of trans-O-glycosidic NMR spin-couplings to either phi (ϕ') or psi (ψ'), and experimental 3JCOCH, 2JCOC, and 3JCOCC spin-couplings measured in aqueous solution in 13C-labeled isotopomers, probability distributions of ϕ' and ψ' in each linkage were determined and compared to those determined by aqueous 1-μs molecular dynamics (MD) simulation. Good agreement was found between the MA'AT and single-state MD conformational models of these linkages for the most part, with modest (approximately <15°) differences in the mean values of ϕ' and ψ', although the envelope of allowed angles (encoded in circular standard deviations or CSDs) is consistently larger for ϕ' determined from MA'AT analysis than from MD for both linkages. The MA'AT model of the α-Galp-(1→3)-β-Galp linkage agrees well with those determined previously using conventional NMR methods (3JCOCH values and/or 1H-1H NOEs), but some discrepancy was observed for the β-Galp-(1→3)-β-Galp linkage, which may arise from errors in the conventions used to describe the linkage torsion angles. Statistical analyses of X-ray crystal structures show ranges of ϕ' and ψ' for both linkages that include the mean angles determined from MA'AT analyses, although both angles adopt a wide range of values in the crystalline state, with ϕ' in β-Galp-(1→3)-β-Galp linkages showing greater-than-expected conformational variability.

Predicting glycan structure from tandem mass spectrometry via deep learning

Glycans constitute the most complicated post-translational modification, modulating protein activity in health and disease. However, structural annotation from tandem mass spectrometry (MS/MS) data is a bottleneck in glycomics, preventing high-throughput endeavors and relegating glycomics to a few experts. Trained on a newly curated set of 500,000 annotated MS/MS spectra, here we present CandyCrunch, a dilated residual neural network predicting glycan structure from raw liquid chromatography-MS/MS data in seconds (top-1 accuracy: 90.3%). We developed an open-access Python-based workflow of raw data conversion and prediction, followed by automated curation and fragment annotation, with predictions recapitulating and extending expert annotation. We demonstrate that this can be used for de novo annotation, diagnostic fragment identification and high-throughput glycomics. For maximum impact, this entire pipeline is tightly interlaced with our glycowork platform and can be easily tested at https://colab.research.google.com/github/BojarLab/CandyCrunch/blob/main/CandyCrunch.ipynb . We envision CandyCrunch to democratize structural glycomics and the elucidation of biological roles of glycans.

Improved physicochemical properties and in vitro digestion of walnut oil microcapsules with soy protein isolate and highly oxidized konjac glucomannan as wall materials

This study investigated the effect of the oxidation degrees of oxidized konjac glucomannan (OKGM) on the encapsulation efficiency (EE), physicochemical and in vitro digestive properties of soy protein isolate (SPI)-based microcapsules walnut oil using experimental and computational approaches. Microcapsules had the highest EE when the ratio of OKGM and SPI to oil was 2.5:1. With increasing the oxidation degree of OKGM, the EE of microcapsules was increased and the hygroscopicity was decreased. Molecular dynamics simulation results showed that SPI/oil/highly OKGM had relatively low binding energy (-4.03 × 106 kJ/mol) and strong electrostatic interactions, which may contribute to a higher EE and lower hygroscopicity of microcapsules, respectively. The oxidative stability of the oil was markedly improved by SPI and OKGM, and microcapsules prepared with SPI and highly OKGM had the highest in vitro digestion. This study provided theoretical support for broadening the application of microcapsules prepared with SPI and OKGM.

Fast Quantitative Validation of 3D Models of Low-Affinity Protein-Ligand Complexes by STD NMR Spectroscopy

Low-affinity protein-ligand interactions are important for many biological processes, including cell communication, signal transduction, and immune responses. Structural characterization of these complexes is also critical for the development of new drugs through fragment-based drug discovery (FBDD), but it is challenging due to the low affinity of fragments for the binding site. Saturation transfer difference (STD) NMR spectroscopy has revolutionized the study of low-affinity receptor-ligand interactions enabling binding detection and structural characterization. Comparison of relaxation and exchange matrix calculations with 1H STD NMR experimental data is essential for the validation of 3D structures of protein-ligand complexes. In this work, we present a new approach based on the calculation of a reduced relaxation matrix, in combination with funnel metadynamics MD simulations, that allows a very fast generation of experimentally STD-NMR-validated 3D structures of low-affinity protein-ligand complexes.

Carbohydrate Structure Database: current state and recent developments

Carbohydrate Structure Database (CSDB) is a curated glycan data collection and a glycoinformatic platform. In this report, its database, analytical, and other components that have appeared for the recent years are reviewed. The major improvements were achieving close-to-full coverage on glycans from microorganisms, launching modules for glycosyltransferases and saccharide conformations, online glycan builder and 3D modeler, NMR simulator, NMR-based structure predictor, and other tools.

Deciphering the conformational dynamics of Myelin Oligodendrocyte glycoprotein in the myelin sheath

Myelin Oligodendrocyte Glycoprotein (MOG) is a transmembrane protein in the myelin sheath. It acts as an auto-antigen under certain unknown conditions causing demyelination, thus resulting in Myelin Oligodendrocyte Glycoprotein Antibody-associated Disease (MOGAD). The significance of glycosylation in the conformational dynamics of the extracellular region (EC1) of the MOG were evident from the previous computational studies. Here, in this study, we performed the molecular dynamics simulation of the entire human MOG in the myelin sheath for 100 ns using the NAMD program. The results indicated that the EC1 and cytoplasmic region (CP) dominate the conformational rigidity of the protein, and enhance its interaction with lipids. This in turn helps in maintaining the myelin integrity in the presence of glycan. The transmembrane regions have reduced interaction with lipids in the glycosylated system. Moreover, the C-terminal extracellular region 2 (EC2) behaves exactly opposite to that of EC1 in the glycan presence. This may be attributed to the glycosylation site in the EC1 region. Hence, not only the region EC1 (having 3 crucial epitopes) but even the CP region were important for understanding the proper function of MOG in the glycan presence.Communicated by Ramaswamy H. Sarma.

Computational Study on the Wealth Conformations of Hepatitis C Virus Envelope Glycoprotein 2 and Role of Its Glycan Coat

Hepatitis C virus (HCV) is a major cause of chronic liver disease and hepatocellular carcinoma. Antibody development efforts mainly revolve around HCV envelope glycoprotein 2 (E2), which mediates host cell entry by interacting with several cell surface receptors, including CD81. We still have limited knowledge about the structural ensembles and the dynamic behavior of both the CD81 binding sites and the glycans on E2. Here, multiple microsecond-long, all-atom molecular dynamics (MD) simulations, as well as a Markov state model (MSM), were performed to provide an atomistic perspective on the dynamic nature of E2 and its glycans. End-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of E2, which may be exploited in therapeutic efforts. Additionally, the Markov state model built from the simulation maps four metastable states for AS412 and three metastable states for the front layer in CD81 binding sites, while binding with HEPC3 would induce a conformation selection for both of them. Overall, this work presents hitherto unseen functional and structural insights into E2 and its glycan coat, providing a new theoretical foundation to control the conformational plasticity of E2 that could be harnessed for vaccine development.

Structural insights into the role and targeting of EGFRvIII

The epidermal growth factor receptor (EGFR) is a well-known oncogenic driver in lung and other cancers. In glioblastoma multiforme (GBM), the EGFR deletion variant III (EGFRvIII) is frequently found alongside EGFR amplification. Agents targeting the EGFR axis have shown limited clinical benefits in GBM and the role of EGFRvIII in GBM is poorly understood. To shed light on the role of EGFRvIII and its potential as a therapeutic target, we determined X-ray crystal structures of a monomeric EGFRvIII extracellular region (ECR). The EGFRvIII ECR resembles the unliganded conformation of EGFR, including the orientation of the C-terminal region of domain II. Domain II is mostly disordered, but the ECR structure is compact. We selected a nanobody with preferential binding to EGFRvIII relative to EGFR and structurally defined an epitope on domain IV that is occluded in the unliganded intact EGFR. These findings suggest new avenues for EGFRvIII targeting in GBM.

Modeling Chromatography Binding through Molecular Dynamics Simulations with Resin Fragments

Accurate atomistic modeling of the interactions of a chromatography resin with a solute can inform the selection of purification conditions for a product, an important problem in the biotech and pharmaceutical industries. We present a molecular dynamics simulation-based approach for the qualitative prediction of interaction sites (specificity) and retention times (affinity) of a protein for a given chromatography resin. We mimicked the resin with an unrestrained ligand composed of the resin headgroup coupled with successively larger fragments of the agarose backbone. The interactions of the ligand with the protein are simulated in an explicit solvent using the Replica Exchange Molecular Dynamics enhanced sampling approach in conjunction with Hydrogen Mass Repartitioning (REMD-HMR). We computed the ligand interaction surface from the simulation trajectories and correlated the features of the interaction surface with experimentally determined retention times. The simulation and analysis protocol were first applied to a series of ubiquitin mutants for which retention times on Capto MMC resin are available. The ubiquitin simulations helped identify the optimal ligand that was used in subsequent simulations on six proteins for which Capto MMC elution times are available. For each of the six proteins, we computed the interaction surface and characterized it in terms of a range of simulation-averaged residue-level physicochemical descriptors. Modeling of the salt concentrations required for elution with respect to the descriptors resulted in a linear fit in terms of aromaphilicity and Kyte-Doolittle hydrophobicity that was robust to outliers, showed high correlation, and correctly ranked the protein elution order. The physics-based model building approach described here does not require a large experimental data set and can be readily applied to different resins and diverse biomolecules.

Targeting host O-linked glycan biosynthesis affects Ebola virus replication efficiency and reveals differential GalNAc-T acceptor site preferences on the Ebola virus glycoprotein

Ebola virus glycoprotein (EBOV GP) is one of the most heavily O-glycosylated viral glycoproteins, yet we still lack a fundamental understanding of the structure of its large O-glycosylated mucin-like domain and to what degree the host O-glycosylation capacity influences EBOV replication. Using tandem mass spectrometry, we identified 47 O-glycosites on EBOV GP and found similar glycosylation signatures on virus-like particle- and cell lysate-derived GP. Furthermore, we performed quantitative differential O-glycoproteomics on proteins produced in wild-type HEK293 cells and cell lines ablated for the three key initiators of O-linked glycosylation, GalNAc-T1, -T2, and -T3. The data show that 12 out of the 47 O-glycosylated sites were regulated, predominantly by GalNAc-T1. Using the glycoengineered cell lines for authentic EBOV propagation, we demonstrate the importance of O-linked glycan initiation and elongation for the production of viral particles and the titers of progeny virus. The mapped O-glycan positions and structures allowed to generate molecular dynamics simulations probing the largely unknown spatial arrangements of the mucin-like domain. The data highlight targeting GALNT1 or C1GALT1C1 as a possible way to modulate O-glycan density on EBOV GP for novel vaccine designs and tailored intervention approaches.IMPORTANCEEbola virus glycoprotein acquires its extensive glycan shield in the host cell, where it is decorated with N-linked glycans and mucin-type O-linked glycans. The latter is initiated by a family of polypeptide GalNAc-transferases that have different preferences for optimal peptide substrates resulting in a spectrum of both very selective and redundant substrates for each isoform. In this work, we map the exact locations of O-glycans on Ebola virus glycoprotein and identify subsets of sites preferentially initiated by one of the three key isoforms of GalNAc-Ts, demonstrating that each enzyme contributes to the glycan shield integrity. We further show that altering host O-glycosylation capacity has detrimental effects on Ebola virus replication, with both isoform-specific initiation and elongation playing a role. The combined structural and functional data highlight glycoengineered cell lines as useful tools for investigating molecular mechanisms imposed by specific glycans and for steering the immune responses in future vaccine designs.

Role of N343 glycosylation on the SARS-CoV-2 S RBD structure and co-receptor binding across variants of concern

Glycosylation of the SARS-CoV-2 spike (S) protein represents a key target for viral evolution because it affects both viral evasion and fitness. Successful variations in the glycan shield are difficult to achieve though, as protein glycosylation is also critical to folding and structural stability. Within this framework, the identification of glycosylation sites that are structurally dispensable can provide insight into the evolutionary mechanisms of the shield and inform immune surveillance. In this work, we show through over 45 μs of cumulative sampling from conventional and enhanced molecular dynamics (MD) simulations, how the structure of the immunodominant S receptor binding domain (RBD) is regulated by N-glycosylation at N343 and how this glycan's structural role changes from WHu-1, alpha (B.1.1.7), and beta (B.1.351), to the delta (B.1.617.2), and omicron (BA.1 and BA.2.86) variants. More specifically, we find that the amphipathic nature of the N-glycan is instrumental to preserve the structural integrity of the RBD hydrophobic core and that loss of glycosylation at N343 triggers a specific and consistent conformational change. We show how this change allosterically regulates the conformation of the receptor binding motif (RBM) in the WHu-1, alpha, and beta RBDs, but not in the delta and omicron variants, due to mutations that reinforce the RBD architecture. In support of these findings, we show that the binding of the RBD to monosialylated ganglioside co-receptors is highly dependent on N343 glycosylation in the WHu-1, but not in the delta RBD, and that affinity changes significantly across VoCs. Ultimately, the molecular and functional insight we provide in this work reinforces our understanding of the role of glycosylation in protein structure and function and it also allows us to identify the structural constraints within which the glycosylation site at N343 can become a hotspot for mutations in the SARS-CoV-2 S glycan shield.

Amorphous solid dispersion of a binary formulation with felodipine and HPMC for 3D printed floating tablets

This study focuses on the combination of three-dimensional printing (3DP) and amorphous solid dispersion (ASD) technologies for the manufacturing of gastroretentive floating tablets. Employing hot melt extrusion (HME) and fused deposition modeling (FDM), the study investigates the development of drug-loaded filaments and 3D printed (3DP) tablets containing felodipine as model drug and hydroxypropyl methylcellulose (HPMC) as the polymeric carrier. Prior to fabrication, solubility parameter estimation and molecular dynamics simulations were applied to predict drug-polymer interactions, which are crucial for ASD formation. Physical bulk and surface characterization complemented the quality control of both drug-loaded filaments and 3DP tablets. The analysis confirmed a successful amorphous dispersion of felodipine within the polymeric matrix. Furthermore, the low infill percentage and enclosed design of the 3DP tablet allowed for obtaining low-density systems. This structure resulted in buoyancy during the entire drug release process until a complete dissolution of the 3DP tablets (more than 8 h) was attained. The particular design made it possible for a single polymer to achieve a zero-order controlled release of the drug, which is considered the ideal kinetics for a gastroretentive system. Accordingly, this study can be seen as an advancement in ASD formulation for 3DP technology within pharmaceutics.

The Legionella collagen-like protein employs a distinct binding mechanism for the recognition of host glycosaminoglycans

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual trimer arrangement with a positively charged external surface and negatively charged solvent exposed internal cavity. Through molecular dynamics simulations, we show how the glycosaminoglycan chondroitin-4-sulphate associates with the Lcl-CTD surface via distinct binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate-binding mechanism.

Structural and genetic basis for the binding of a mouse monoclonal antibody to Flavobacterium psychrophilum lipopolysaccharide

A mouse monoclonal antibody (mAb FL100A) previously prepared against Flavobacterium psychrophilum (Fp) CSF259-93 has now been examined for binding to lipopolysaccharides (LPS) of this strain and Fp 950106-1/1. The corresponding O-polysaccharides (O-PS) of these strains are formed by identical trisaccharide repeats composed of l-Rhamnose (l-Rha), 2-acetamido-2-deoxy-l-fucose (l-FucNAc) and 2-acetamido-4-R1-2,4-dideoxy-d-quinovose (d-Qui2NAc4NR1) where R1 represents a dihydroxyhexanamido moiety. The O-PS loci of these strains are also identical except for the gene (wzy1 or wzy2) that encodes the polysaccharide polymerase. Accordingly, adjacent O-PS repeats are joined through d-Qui2NAc4NR1 and l-Rha by wzy2-dependent α(1-2) linkages in Fp CSF259-93 versus wzy1-dependent β(1-3) linkages in Fp 950106-1/1. mAb FL100A reacted strongly with Fp CSF259-93 O-PS and LPS but weakly or not at all with Fp 950106-1/1 LPS and O-PS. Importantly, it also labelled cell surface blebs on the former but not the latter strain. Additionally, mAb binding was approximately 5-times stronger to homologous Fp CSF259-93 LPS than to LPS from a strain with a different R-group gene. A conformational epitope for mAb FL100A binding was suggested from molecular dynamic simulations of each O-PS. Thus, Fp CSF259-93 O-PS formed a stable well-defined compact helix in which the R1 groups were displayed in a regular pattern on the helix exterior while unreactive Fp 950106-1/1 O-PS adopted a flexible extended linear conformation. Taken together, the findings establish the specificity of mAb FL100A for Wzy2-linked F. psychrophilum O-PS and LPS.

Evaluating solubility, stability, and inclusion complexation of oxyresveratrol with various β-cyclodextrin derivatives using advanced computational techniques and experimental validation

Oxyresveratrol (OXY), a natural stilbenoid in mulberry fruits, is known for its diverse pharmacological properties. However, its clinical use is hindered by low water solubility and limited bioavailability. In the present study, the inclusion complexes of OXY with β-cyclodextrin (βCD) and its three analogs, dimethyl-β-cyclodextrin (DMβCD), hydroxypropyl-β-cyclodextrin (HPβCD) and sulfobutylether-β-cyclodextrin (SBEβCD), were investigated using in silico and in vitro studies. Molecular docking revealed two binding orientations of OXY, namely, 4',6'-dihydroxyphenyl (A-form) and 5,7-benzenediol ring (B-form). Molecular Dynamics simulations suggested the formation of inclusion complexes with βCDs through two distinct orientations, with OXY/SBEβCD exhibiting maximum atom contacts and the lowest solvent-exposed area in the hydrophobic cavity. These results corresponded well with the highest binding affinity observed in OXY/SBEβCD when assessed using the MM/GBSA method. Beyond traditional simulation methods, Ligand-binding Parallel Cascade Selection Molecular Dynamics method was employed to investigate how the drug enters and accommodates within the hydrophobic cavity. The in silico results aligned with stability constants: SBEβCD (2060 M-1), HPβCD (1860 M-1), DMβCD (1700 M-1), and βCD (1420 M-1). All complexes exhibited a 1:1 binding mode (AL type), with SBEβCD enhancing OXY solubility (25-fold). SEM micrographs, DSC thermograms, FT-IR and 1H NMR spectra confirm the inclusion complex formation, revealing novel surface morphologies, distinctive thermal behaviors, and new peaks. Notably, the inhibitory impact on the proliferation of breast cancer cell lines, MCF-7, exhibited by inclusion complexes particularly OXY/DMβCD, OXY/HPβCD, and OXY/SBEβCD were markedly superior compared to that of OXY alone.

Preparation and characterization of cellulose fluorescent material: Experiment and simulation

The extensive use of fossil based materials has caused serious pollution problems, the full utilization of biomass resources to prepare high value-added new materials is of great significance for the environmental protection and sustainable social development. For this purpose, this study explored the preparation process and molecular dynamics simulation of cellulose fluorescent materials. Firstly, bacterial cellulose was dissolved in a solution of NaOH and urea at low temperature, followed by a solution blending and hot pressing with hyperbranched polyamide. It was found that the addition of hyperbranched polyamide could effectively filled in the internal pores of cellulose hydrogel, thereby enhancing the fluorescence effects and tensile properties, especially the elongation at break of cellulose materials. The optimal amount of hyperbranched polyamide added was 5 wt%. Molecular dynamics simulation showed that the hydrogen bonds and interaction with cellulose increased as the concentration of hyperbranched polyamide increased.

Growing Glycans in Rosetta: Accurate de novo glycan modeling, density fitting, and rational sequon design

Carbohydrates and glycoproteins modulate key biological functions. However, experimental structure determination of sugar polymers is notoriously difficult. Computational approaches can aid in carbohydrate structure prediction, structure determination, and design. In this work, we developed a glycan-modeling algorithm, GlycanTreeModeler, that computationally builds glycans layer-by-layer, using adaptive kernel density estimates (KDE) of common glycan conformations derived from data in the Protein Data Bank (PDB) and from quantum mechanics (QM) calculations. GlycanTreeModeler was benchmarked on a test set of glycan structures of varying lengths, or "trees". Structures predicted by GlycanTreeModeler agreed with native structures at high accuracy for both de novo modeling and experimental density-guided building. We employed these tools to design de novo glycan trees into a protein nanoparticle vaccine to shield regions of the scaffold from antibody recognition, and experimentally verified shielding. This work will inform glycoprotein model prediction, glycan masking, and further aid computational methods in experimental structure determination and refinement.

Rotational Dynamics of Water near Osmolytes by Molecular Dynamics Simulations

The behavior of water molecules around organic molecules has attracted considerable attention as a crucial factor influencing the properties and functions of soft matter and biomolecules. Recently, it has been suggested that the change in protein stability upon the addition of small organic molecules (osmolytes) is dominated by the change in the water dynamics caused by the osmolyte, where the dynamics of not only the directly interacting water molecules but also the long-range hydration layer affect the protein stability. However, the relation between the long-range structure of hydration water in various solutions and the water dynamics remains unclear at the molecular level. We performed density-functional tight-binding molecular dynamics simulations to elucidate the varying rotational dynamics of water molecules in 15 osmolyte solutions. A positive correlation was observed between the rotational relaxation time and our proposed normalized parameter obtained by dividing the number of hydrogen bonds between water molecules by the number of nearest-neighbor water molecules. For the 15 osmolyte solutions, an increase or a decrease in the value of the normalized parameter for the second hydration shell tended to result in water molecules with slow and fast rotational dynamics, respectively, thus illustrating the importance of the second hydration shell for the rotational dynamics of water molecules. Our simulation results are anticipated to advance the current understanding of water dynamics around organic molecules and the long-range structure of water molecules.

Enhancing Protein Solubility via Glycosylation: From Chemical Synthesis to Machine Learning Predictions

Glycosylation is a valuable tool for modulating protein solubility; however, the lack of reliable research strategies has impeded efficient progress in understanding and applying this modification. This study aimed to bridge this gap by investigating the solubility of a model glycoprotein molecule, the carbohydrate-binding module (CBM), through a two-stage process. In the first stage, an approach involving chemical synthesis, comparative analysis, and molecular dynamics simulations of a library of glycoforms was employed to elucidate the effect of different glycosylation patterns on solubility and the key factors responsible for the effect. In the second stage, a predictive mathematical formula, innovatively harnessing machine learning algorithms, was derived to relate solubility to the identified key factors and accurately predict the solubility of the newly designed glycoforms. Demonstrating feasibility and effectiveness, this two-stage approach offers a valuable strategy for advancing glycosylation research, especially for the discovery of glycoforms with increased solubility.

Interfacial Glucose to Regulate Hydrated Lipid Bilayer Properties: Influence of Concentrations

A series of atomistic molecular dynamics (MD) simulations were carried out with a hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer with the variation of glucose concentrations from 0 to 30 wt % in the presence of 0.3 M NaCl. The study suggested that although the thickness of the lipid bilayer dropped significantly with the increase in glucose concentration, it expanded laterally at high glucose levels due to the intercalation of glucose between the headgroups of adjacent lipids. We adopted the surface assessment via the grid evaluation method to compute the deviation of the bilayer's key structural features for the different amounts of glucose present. This suggested that the accumulation of glucose molecules near the headgroups influences the local lipid bilayer undulation and crimping of the lipid tails. We find that the area compressibility modulus increases with the glucose level, causing enhanced bilayer rigidity arising from the slow lateral diffusion of lipids. The restricted lipid motion at high glucose concentrations controls the sustainability of the curved bilayer surface. Calculations revealed that certain orientations of CO → of interfacial glucose with the PN → of lipid headgroups are preferred, which helps the glucose to form direct hydrogen bonds (HBs) with the lipid headgroups. Such lipid-glucose (LG) HBs relax slowly at low glucose concentrations and exhibit a higher lifetime, whereas fast structural relaxation of LG HBs with a shorter lifetime was noticed at a higher glucose level. In contrast, lipid-water (LW) HBs exhibited a higher lifetime at a higher glucose level, which gradually decreased with the glucose level lowering. The study interprets that the glucose concentration-driven LW and LG interactions are mutually inclusive. Our detailed analysis will exemplify small saccharide concentration-driven membrane stabilizing efficiency, which is, in general, helpful for drug delivery study.