RESUMO
Carbohydrate-binding modules (CBMs) are protein domains that typically reside near catalytic domains, increasing substrate-protein proximity by constraining the conformational space of carbohydrates. Due to the flexibility and variability of glycans, the molecular details of how these protein regions recognize their target molecules are not always fully understood. Computational methods, including molecular docking and molecular dynamics simulations, have been employed to investigate lectin-carbohydrate interactions. In this study, we introduce a novel approach that integrates multiple computational techniques to identify the critical amino acids involved in the interaction between a CBM located at the tip of bacteriophage J-1's tail and its carbohydrate counterparts. Our results highlight three amino acids that play a significant role in binding, a finding we confirmed through in vitro experiments. By presenting this approach, we offer an intriguing alternative for pinpointing amino acids that contribute to protein-sugar interactions, leading to a more thorough comprehension of the molecular determinants of protein-carbohydrate interactions.
RESUMO
Unraveling the structure of lectin-carbohydrate complexes is vital for understanding key biological recognition processes and development of glycomimetic drugs. Molecular Docking application to predict them is challenging due to their low affinity, hydrophilic nature and ligand conformational diversity. In the last decade several strategies, such as the inclusion of glycan conformation specific scoring functions or our developed solvent-site biased method, have improved carbohydrate docking performance but significant challenges remain, in particular, those related to receptor conformational diversity. In the present work we have analyzed conventional and solvent-site biased autodock4 performance concerning receptor conformational diversity as derived from different crystal structures (apo and holo), Molecular Dynamics snapshots and Homology-based models, for 14 different lectin-monosaccharide complexes. Our results show that both conventional and biased docking yield accurate lectin-monosaccharide complexes, starting from either apo or homology-based structures, even when only moderate (45%) sequence identity templates are available. An essential element for success is a proper combination of a middle-sized (10-100 structures) conformational ensemble, derived either from Molecular dynamics or multiple homology model building. Consistent with our previous works, results show that solvent-site biased methods improve overall performance, but that results are still highly system dependent. Finally, our results also show that docking can select the correct receptor structure within the ensemble, underscoring the relevance of joint evaluation of both ligand pose and receptor conformation.
Assuntos
Lectinas/química , Modelos Moleculares , Monossacarídeos/química , Cristalografia por Raios XRESUMO
Misfolded glycoprotein recognition and endoplasmic reticulum (ER) retention are mediated by the ER glycoprotein folding quality control (ERQC) checkpoint enzyme, UDP-glucose glycoprotein glucosyltransferase (UGGT). UGGT modulation is a promising strategy for broad-spectrum antivirals, rescue-of-secretion therapy in rare disease caused by responsive mutations in glycoprotein genes, and many cancers, but to date no selective UGGT inhibitors are known. The small molecule 5-[(morpholin-4-yl)methyl]quinolin-8-ol (5M-8OH-Q) binds a CtUGGTGT24 "WY" conserved surface motif conserved across UGGTs but not present in other GT24 family glycosyltransferases. 5M-8OH-Q has a 47 µM binding affinity for CtUGGTGT24in vitro as measured by ligand-enhanced fluorescence. In cellula, 5M-8OH-Q inhibits both human UGGT isoforms at concentrations higher than 750 µM. 5M-8OH-Q binding to CtUGGTGT24 appears to be mutually exclusive to M5-9 glycan binding in an in vitro competition experiment. A medicinal program based on 5M-8OH-Q will yield the next generation of UGGT inhibitors.
RESUMO
UDP-glucose:glycoprotein glucosyltransferase (UGGT) flags misfolded glycoproteins for ER retention. We report crystal structures of full-length Chaetomium thermophilum UGGT (CtUGGT), two CtUGGT double-cysteine mutants, and its TRXL2 domain truncation (CtUGGT-ΔTRXL2). CtUGGT molecular dynamics (MD) simulations capture extended conformations and reveal clamping, bending, and twisting inter-domain movements. We name "Parodi limit" the maximum distance on the same glycoprotein between a site of misfolding and an N-linked glycan that can be reglucosylated by monomeric UGGT in vitro, in response to recognition of misfold at that site. Based on the MD simulations, we estimate the Parodi limit as around 70-80 Å. Frequency distributions of distances between glycoprotein residues and their closest N-linked glycosylation sites in glycoprotein crystal structures suggests relevance of the Parodi limit to UGGT activity in vivo. Our data support a "one-size-fits-all adjustable spanner" UGGT substrate recognition model, with an essential role for the UGGT TRXL2 domain.
Assuntos
Proteínas Fúngicas/química , Glucosiltransferases/química , Simulação de Dinâmica Molecular , Domínio Catalítico , Chaetomium/enzimologia , Proteínas Fúngicas/metabolismo , Glucosiltransferases/metabolismo , Glicoproteínas/química , Glicoproteínas/metabolismo , Células HEK293 , Humanos , Dobramento de ProteínaRESUMO
Understanding enzymatic reactions with atomic resolution has proven in recent years to be of tremendous interest for biochemical research, and thus, the use of QM/MM methods for the study of reaction mechanisms is experiencing a continuous growth. Glycosyltransferases (GTs) catalyze the formation of glycosidic bonds, and are important for many biotechnological purposes, including drug targeting. Their reaction product may result with only one of the two possible stereochemical outcomes for the reacting anomeric center, and therefore, they are classified as either inverting or retaining GTs. While the inverting GT reaction mechanism has been widely studied, the retaining GT mechanism has always been controversial and several questions remain open to this day. In this work, we take advantage of our recent GPU implementation of a pure QM(DFT-PBE)/MM approach to explore the reaction and inhibition mechanism of MshA, a key retaining GT responsible for the first step of mycothiol biosynthesis, a low weight thiol compound found in pathogens like Mycobacterium tuberculosis that is essential for its survival under oxidative stress conditions. Our results show that the reaction proceeds via a front-side SNi-like concerted reaction mechanism (DNAN in IUPAC nomenclature) and has a 17.5 kcal/mol free energy barrier, which is in remarkable agreement with experimental data. Detailed analysis shows that the key reaction step is the diphosphate leaving group dissociation, leading to an oxocarbenium-ion-like transition state. In contrast, fluorinated substrate analogues increase the reaction barrier significantly, rendering the enzyme effectively inactive. Detailed analysis of the electronic structure along the reaction suggests that this particular inhibition mechanism is associated with fluorine's high electronegative nature, which hinders phosphate release and proper stabilization of the transition state.