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
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.