Selectins and their ability to interact with specific ligands are a cornerstone in cell communication. Over the last three decades, a considerable wealth of experimental and molecular modeling insights into their structure and modus operandi were gathered. Nonetheless, explaining the role of individual selectin residues on a quantitative level remained elusive, despite its importance in understanding the structure-function relationship in these molecules and designing their inhibitors. This work explores essential interactions of selectin-ligand binding, employing a multiscale approach that combines molecular dynamics, quantum-chemical calculations, and residue interaction network models. Such an approach successfully reproduces most of the experimental findings. It proves to be helpful, with the potential for becoming an established tool for quantitative predictions of residue contribution to the binding of biomolecular complexes. The results empower us to quantify the importance of particular residues and functional groups in the protein-ligand interface and to pinpoint differences in molecular recognition by the three selectins. We show that mutations in the E-, L-, and P-selectins, e.g., different residues in positions 46, 85, 97, and 107, present a crucial difference in how the ligand is engaged. We assess the role of sulfation of tyrosine residues in PSGL-1 and suggest that TyrSO3- in position 51 interacting with Arg85 in P-selectin is a significant factor in the increased affinity of P-selectin to PSGL-1 compared to E- and L-selectins. We propose an original pharmacophore targeting five essential PSGL-binding sites based on the analysis of the selectin···PSGL-1 interactions.

• MeSH
• antigen sialyl-Lewis-X MeSH
• buněčná adheze MeSH
• ligandy MeSH
• P-selektin * metabolismus   MeSH
• selektiny * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• antigen sialyl-Lewis-X MeSH
• ligandy MeSH
• P-selektin * MeSH
• selektiny * MeSH

Tyrosine sulfation alters the biological activity of many proteins involved in different physiological and pathophysiological conditions, such as non-specific immune reaction, response to inflammation and ischemia, targeting of leukocytes and stem cells, or the formation of cancer metastases. Tyrosine sulfation is catalyzed by the enzymes tyrosylprotein sulfotransferases (TPST). In this study, we used QM/MM Car-Parrinello metadynamics simulations together with QM/MM potential energy calculations to investigate the catalytic mechanism of isoform TPST-1. The structural changes along the reaction coordinate are analyzed and discussed. Furthermore, both the methods supported the SN2 type of catalytic mechanism. The reaction barrier obtained from CPMD metadynamics was 12.8 kcal mol-1, and the potential energy scan led to reaction barriers of 11.6 kcal mol-1 and 13.7 kcal mol-1 with the B3LYP and OPBE functional, respectively. The comparison of the two methods (metadynamics and potential energy scan) may be helpful for future mechanistic studies. The insight into the reaction mechanism of TPST-1 might help with the rational design of transition-state TPST inhibitors.

• MeSH
• biokatalýza MeSH
• katalytická doména MeSH
• kvantová teorie MeSH
• lidé MeSH
• simulace molekulární dynamiky MeSH
• sulfotransferasy chemie   MeSH
• termodynamika MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• protein-tyrosine sulfotransferase MeSH Prohlížeč
• sulfotransferasy MeSH

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.

• Klíčová slova
• COVID-19, docking, pan-selectin inhibitors, selectins, virtual screening,
• MeSH
• antivirové látky chemie   MeSH
• biomimetické materiály chemie   MeSH
• chemické databáze * 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
• Názvy látek
• antivirové látky MeSH
• selektiny MeSH

We applied the transition path sampling (TPS) method to study the translocation step of the catalytic mechanism of galactofuranosyl transferase 2 (GlfT2). Using TPS in the field of enzymatic reactions is still relatively rare, and we show its effectiveness on this enzymatic system. We decipher an unknown mechanism of the translocation step and, thus, provide a complete understanding of the catalytic mechanism of GlfT2 at the atomistic level. The GlfT2 enzyme is involved in the formation of the mycobacterial cell wall and transfers galactofuranose (Galf) from UDP-Galf onto a growing acceptor Galf chain. The biosynthesis of the galactan chain is accomplished in a processive manner, with the growing acceptor substrate remaining bound to GlfT2. The glycosidic bond formed by GlfT2 between the two Galf residues alternates between β-(1-6) and β-(1-5) linkages. The translocation of the growing galactan between individual additions of Galf residues is crucial for the function of GlfT2. Analysis of unbiased trajectory ensembles revealed that the translocation proceeds differently depending on the glycosidic linkage between the last two Galf residues. We also showed that the protonation state of the catalytic residue Asp372 significantly influences the translocation. Approximate transition state structures and potential energy reaction barriers of the translocation process were determined. The calculated potential reaction barriers in the range of 6-14 kcal/mol show that the translocation process is not the rate-limiting step in galactan biosynthesis.

• Publikační typ
• časopisecké články 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.

• Klíčová slova
• cancer, carbohydrate determinants, glycosyltransferases, inflammation, inhibitors, selectin ligands, selectin-ligand interactions, selectins, sialyl Lewis x, transition state analog inhibitors,
• 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
• Názvy látek
• nádorové proteiny MeSH
• selektiny 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.

• MeSH
• chemické modely * MeSH
• kationické antimikrobiální peptidy chemie   MeSH
• oligosacharidy chemie   MeSH
• proteinové domény MeSH
• rostlinné lektiny chemie   MeSH
• simulace molekulární dynamiky * MeSH
• vazba proteinů MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• hevein MeSH Prohlížeč
• kationické antimikrobiální peptidy MeSH
• oligosacharidy MeSH
• rostlinné lektiny MeSH

Mycobacterium tuberculosis features a unique cell wall that protects the bacterium from the external environment. Disruption of the cell wall assembly is a promising direction for novel anti-tuberculotic drugs. A key component of the cell wall is galactan, a polysaccharide chain composed of galactofuranose (Galf) units connected by alternating β-(1-5) and β-(1-6) linkages. The majority of the galactan chain is biosynthesized by a bifunctional enzyme-galactofuranosyl transferase 2 (GlfT2). GlfT2 catalyzes two reactions: the formation of β-(1-5) and β-(1-6) linkages. It was suggested that the enzyme acts through a processive mechanism until it adds 30-35 Galf units in a single active site. We applied a QM/MM string method coupled with ab initio molecular dynamics simulations to study the two reactions catalyzed by GlfT2. We showed that both reactions proceed very similarly and feature similar transition-state structures. We also present novel information about the ring puckering behavior of the five-membered furanose ring during the glycosyltransferase reaction and a calculated transition-state structure with galactose in a furanose form that may be used as a guide for the rational design of very specific and extremely potent inhibitors, that is, transition-state analogues, for GlfT2. Due to the absence of a furanose form of galactose in humans, transition-state-analogous inhibitors represent an attractive scaffold for the development of novel antibacterial drugs.

• Klíčová slova
• enzyme catalysis, glycosyltransferases, molecular modeling, reaction mechanisms, tuberculosis,
• MeSH
• bakteriální proteiny chemie   metabolismus   MeSH
• galaktosa chemie   metabolismus   MeSH
• galaktosyltransferasy chemie   metabolismus   MeSH
• katalytická doména MeSH
• kvantová teorie * MeSH
• Mycobacterium tuberculosis enzymologie   MeSH
• simulace molekulární dynamiky * MeSH
• substrátová specifita MeSH
• termodynamika MeSH
• vazebná místa MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• bakteriální proteiny MeSH
• galaktosa MeSH
• galaktosyltransferasy MeSH

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.

• Klíčová slova
• GlcNAc-TV, N-glycosylation, glycosyltransferase, homology model, molecular dynamics,
• 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
• Názvy látek
• alpha-1,6-mannosylglycoprotein beta 1,6-N-acetylglucosaminyltransferase MeSH Prohlížeč
• glykosyltransferasy MeSH
• mannosa MeSH
• N-acetylglukosaminyltransferasy MeSH
• polysacharidy 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.

• Klíčová slova
• CH/π interactions, Carbohydrate–aromatic interactions, Hevein domain, Hydrophobic effect, Molecular dynamics simulation,
• 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
• Názvy látek
• hevein MeSH Prohlížeč
• kationické antimikrobiální peptidy MeSH
• rostlinné lektiny MeSH
• trisacharidy MeSH
• vodík 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
• Názvy látek
• glykosyltransferasy MeSH
• polysacharidy MeSH