Coronavirus disease 2019 (COVID-19) has spread rapidly throughout the globe. The spectrum of disease is broad but among hospitalized patients with COVID-19, respiratory failure from acute respiratory distress syndrome is the leading cause of mortality. There is an urgent need for an effective treatment. The current focus has been developing novel therapeutics, including antivirals, protease inhibitors, vaccines and targeting the overactive cytokine response with anti-cytokine therapy. The overproduction of early response proinflammatory cytokines results in what has been described as a "cytokine storm" is leading eventually to death when the cells fail to terminate the inflammatory response. Accumulating evidence shows that inflammatory cytokines induce selectin ligands that play a crucial role in the pathogenesis of inflammatory diseases by mediating leukocyte migration from the blood into the tissue. Thus, the selectins and selectin ligands represent a promising therapeutic target for the treatment of COVID-19. In this paper, potential pan-selectin inhibitors were identified employing a virtual screening using a docking procedure. For this purpose, the Asinex and ZINC databases of ligands, including approved drugs, biogenic compounds and glycomimetics, altogether 923,602 compounds, were screened against the P-, L- and E-selectin. At first, the experimentally confirmed inhibitors were docked into all three selectins' carbohydrate recognition domains to assess the suitability of the screening procedure. Finally, based on the evaluation of ligands binding, we propose 10 purchasable pan-selectin inhibitors to develop COVID-19 therapeutics.
- MeSH
- antivirové látky chemie MeSH
- biomimetické materiály chemie MeSH
- chemické databáze * MeSH
- COVID-19 MeSH
- farmakoterapie COVID-19 MeSH
- lidé MeSH
- počítačová simulace * MeSH
- preklinické hodnocení léčiv MeSH
- SARS-CoV-2 chemie metabolismus MeSH
- selektiny chemie MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Selectins belong to a group of adhesion molecules that fulfill an essential role in immune and inflammatory responses and tissue healing. Selectins are glycoproteins that decode the information carried by glycan structures, and non-covalent interactions of selectins with these glycan structures mediate biological processes. The sialylated and fucosylated tetrasaccharide sLex is an essential glycan recognized by selectins. Several glycosyltransferases are responsible for the biosynthesis of the sLex tetrasaccharide. Selectins are involved in a sequence of interactions of circulated leukocytes with endothelial cells in the blood called the adhesion cascade. Recently, it has become evident that cancer cells utilize a similar adhesion cascade to promote metastases. However, like Dr. Jekyll and Mr. Hyde's two faces, selectins also contribute to tissue destruction during some infections and inflammatory diseases. The most prominent function of selectins is associated with the initial stage of the leukocyte adhesion cascade, in which selectin binding enables tethering and rolling. The first adhesive event occurs through specific non-covalent interactions between selectins and their ligands, with glycans functioning as an interface between leukocytes or cancer cells and the endothelium. Targeting these interactions remains a principal strategy aimed at developing new therapies for the treatment of immune and inflammatory disorders and cancer. In this review, we will survey the significant contributions to and the current status of the understanding of the structure of selectins and the role of selectins in various biological processes. The potential of selectins and their ligands as therapeutic targets in chronic and acute inflammatory diseases and cancer will also be discussed. We will emphasize the structural characteristic of selectins and the catalytic mechanisms of glycosyltransferases involved in the biosynthesis of glycan recognition determinants. Furthermore, recent achievements in the synthesis of selectin inhibitors will be reviewed with a focus on the various strategies used for the development of glycosyltransferase inhibitors, including substrate analog inhibitors and transition state analog inhibitors, which are based on knowledge of the catalytic mechanism.
- MeSH
- buněčná adheze * MeSH
- leukocyty metabolismus patologie MeSH
- lidé MeSH
- nádorové proteiny metabolismus MeSH
- nádory metabolismus patologie MeSH
- rolling leukocytů * MeSH
- selektiny metabolismus MeSH
- zánět metabolismus patologie MeSH
- zvířata MeSH
- Check Tag
- lidé MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- přehledy MeSH
Interactions between proteins and their small molecule ligands are of great importance for the process of drug design. Here we report an unbiased molecular dynamics simulation of systems containing hevein domain (HEV32) with N-acetylglucosamine mono-, di- or trisaccharide. Carbohydrate molecules were placed outside the binding site. Three of six simulations (6 × 2 μs) led to binding of a carbohydrate ligand into the binding mode in agreement with the experimentally determined structure. Unbinding was observed in one simulation (monosaccharide). There were no remarkable intermediates of binding for mono and disaccharide. Trisaccharide binding was initiated by formation of carbohydrate-aromatic CH/π interactions. Our results indicate that binding of ligands followed the model of conformational selection because the conformation of the protein ready for ligand binding was observed before the binding. This study extends the concept of docking by dynamics on carbohydrate-protein interactions.
The enzyme UDP-N-acetylglucosamine: α-d-mannoside β-1-6 N-acetylglucosaminyltransferase V (GnT-V) catalyzes the transfer of GlcNAc from the UDP-GlcNAc donor to the α-1-6-linked mannose of the trimannosyl core structure of glycoproteins to produce the β-1-6-linked branching of N-linked oligosaccharides. β-1-6-GlcNAc-branched N-glycans are associated with cancer growth and metastasis. Therefore, the inhibition of GnT-V represents a key target for anti-cancer drug development. However, the development of potent and specific inhibitors of GnT-V is hampered by the lack of information on the three-dimensional structure of the enzyme and on the binding characteristics of its substrates. Here we present the first 3D structure of GnT-V as a result of homology modeling. Various alignment methods, docking the donor and acceptor substrates, and molecular dynamics simulation were used to construct seven homology models of GnT-V and characterize the binding of its substrates. The best homology model is consistent with available experimental data. The three-dimensional model, the structure of the enzyme catalytic site and binding information obtained for the donor and acceptor can be useful in studies of the catalytic mechanism and design of inhibitors of GnT-V.
- MeSH
- glykosyltransferasy chemie MeSH
- lidé MeSH
- mannosa chemie MeSH
- molekulární konformace * MeSH
- molekulární modely MeSH
- N-acetylglukosaminyltransferasy chemie MeSH
- polysacharidy chemie genetika MeSH
- simulace molekulární dynamiky MeSH
- substrátová specifita MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
The glycosylation of cell surface proteins plays a crucial role in a multitude of biological processes, such as cell adhesion and recognition. To understand the process of protein glycosylation, the reaction mechanisms of the participating enzymes need to be known. However, the reaction mechanism of retaining glycosyltransferases has not yet been sufficiently explained. Here we investigated the catalytic mechanism of human isoform 2 of the retaining glycosyltransferase polypeptide UDP-GalNAc transferase by coupling two different QM/MM-based approaches, namely a potential energy surface scan in two distance difference dimensions and a minimum energy reaction path optimisation using the Nudged Elastic Band method. Potential energy scan studies often suffer from inadequate sampling of reactive processes due to a predefined scan coordinate system. At the same time, path optimisation methods enable the sampling of a virtually unlimited number of dimensions, but their results cannot be unambiguously interpreted without knowledge of the potential energy surface. By combining these methods, we have been able to eliminate the most significant sources of potential errors inherent to each of these approaches. The structural model is based on the crystal structure of human isoform 2. In the QM/MM method, the QM region consists of 275 atoms, the remaining 5776 atoms were in the MM region. We found that ppGalNAcT2 catalyzes a same-face nucleophilic substitution with internal return (SNi). The optimized transition state for the reaction is 13.8 kcal/mol higher in energy than the reactant while the energy of the product complex is 6.7 kcal/mol lower. During the process of nucleophilic attack, a proton is synchronously transferred to the leaving phosphate. The presence of a short-lived metastable oxocarbenium intermediate is likely, as indicated by the reaction energy profiles obtained using high-level density functionals.
- MeSH
- algoritmy MeSH
- chemické modely * MeSH
- glykosylace MeSH
- glykosyltransferasy chemie ultrastruktura MeSH
- katalýza MeSH
- kinetika MeSH
- konformace proteinů MeSH
- molekulární modely * MeSH
- počítačová simulace MeSH
- polysacharidy chemie ultrastruktura MeSH
- vazba proteinů MeSH
- vazebná místa MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Carbohydrate-protein complexes are often characterized by interactions via aromatic amino acid residues. Several mechanisms have been proposed to explain these stacking-like interactions between pyranose sugars and aromatic moieties. The physical basis of these interactions is being explained as either dispersion CH/π or hydrophobic. In order to elucidate the nature of these interactions, we performed a series of molecular dynamics simulation of hevein domain (HEV32) in complex with (β-D-GlcNAc)3. Selected OH/O and CH/π hydrogen bonds involved in carbohydrate recognition were artificially weakened in 100 ns molecular dynamics simulations. Separate weakening of either OH/O or CH/π hydrogen bonds was not sufficient to destabilize the complex. This indicates that other effects, not solely CH/π dispersion interactions, contribute significantly to the stability of the complex. Significant destabilization of complexes was reached only by simultaneous weakening of OH/O and CH/π hydrogen bonds. This also shows that classical hydrogen bonds and CH/π interactions are working in concert to stabilize this carbohydrate-protein test case.
- MeSH
- kationické antimikrobiální peptidy chemie metabolismus MeSH
- konformace sacharidů MeSH
- molekulární modely MeSH
- rostlinné lektiny chemie metabolismus MeSH
- sacharidové sekvence MeSH
- simulace molekulární dynamiky MeSH
- trisacharidy chemie MeSH
- vazebná místa MeSH
- vodík chemie MeSH
- vodíková vazba MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Protein-ligand affinities can be significantly influenced not only by the interaction itself but also by conformational equilibrium of both binding partners, free ligand and free protein. Identification of important conformational families of a ligand and prediction of their thermodynamics is important for efficient ligand design. Here we report conformational free energy modeling of nine small-molecule drugs in explicitly modeled water by metadynamics with a bias potential applied in the space of weighted holistic invariant molecular (WHIM) descriptors. Application of metadynamics enhances conformational sampling compared to unbiased molecular dynamics simulation and allows to predict relative free energies of key conformations. Selected free energy minima and one example of transition state were tested by a series of unbiased molecular dynamics simulation. Comparison of free energy surfaces of free and target-bound Imatinib provides an estimate of free energy penalty of conformational change induced by its binding to the target.
Modelling of carbohydrate conformations is a challenging task for force field developers. Three carbohydrate force fields, namely GLYCAM06, GROMOS 45a4 and OPLS were evaluated. Free energies of different ring conformations of beta-D-glucopyranose were calculated using metadynamics in vacuum as well as in explicitly modelled water. All three force fields model the (4)C(1) conformation as the most stable by at least 6kJ/mol, as compared to other conformations. Interconversion from the (4)C(1) to any other conformation is associated with a barrier of no lower than 26kJ/mol. The free energy surface calculated in the GLYCAM06 force field is in remarkably good agreement with the recent Car-Parrinello metadynamics study. The effect of a water environment is relatively low and analogous in all tested force fields. Namely, the presence of water stabilizes the upper-left ((3,O)B) versus bottom-right (B(3,O)) area of Stoddard's plot, relative to the situation in vacuum. Comparison of free and potential surfaces is also provided for vacuum calculations.
Conformational behavior of the [(2S,3R,4R,5S)-3,4,5-trihydroxy-2-(phenylsulfanyl)tetrahydrofuran-2-yl]methyl sulfate anion (2), which is the potential transition state (TS) analogue of the inverting glycosyltransferases, was studied by means of two-dimensional potential-energy maps, using a density functional theory method at the B3LYP/6-31+G* level. The maps revealed the presence of eight low-energy domains which were refined at the B3LYP/6-311++G** level and led to six conformers in vacuum. In aqueous solution, two conformers dominate at equilibrium. The preferred conformers superimpose well with the transition state structure, as determined previously for glycosyltransferase GnT-I. The conformations of 2 in the active site of glycosyltransferase GnT-I were obtained by docking methods. It was found that one of the two best docking poses mimics the binding mode of TS. These results suggest that the proposed TS mimics 2 have the potential to be used as a scaffold for the design of TS analogue inhibitors.
