Importance of ligand conformational energies in carbohydrate docking: Sorting the wheat from the chaff

Vina-Carb: Improving Glycosidic Angles during Carbohydrate Docking

Molecular docking programs are primarily designed to align rigid, drug-like fragments into the binding sites of macromolecules and frequently display poor performance when applied to flexible carbohydrate molecules. A critical source of flexibility within an oligosaccharide is the glycosidic linkages. Recently, Carbohydrate Intrinsic (CHI) energy functions were reported that attempt to quantify the glycosidic torsion angle preferences. In the present work, the CHI-energy functions have been incorporated into the AutoDock Vina (ADV) scoring function, subsequently termed Vina-Carb (VC). Two user-adjustable parameters have been introduced, namely, a CHI- energy weight term (chi_coeff) that affects the magnitude of the CHI-energy penalty and a CHI-cutoff term (chi_cutoff) that negates CHI-energy penalties below a specified value. A data set consisting of 101 protein-carbohydrate complexes and 29 apoprotein structures was used in the development and testing of VC, including antibodies, lectins, and carbohydrate binding modules. Accounting for the intramolecular energies of the glycosidic linkages in the oligosaccharides during docking led VC to produce acceptable structures within the top five ranked poses in 74% of the systems tested, compared to a success rate of 55% for ADV. An enzyme system was employed in order to illustrate the potential application of VC to proteins that may distort glycosidic linkages of carbohydrate ligands upon binding. VC represents a significant step toward accurately predicting the structures of protein-carbohydrate complexes. Furthermore, the described approach is conceptually applicable to any class of ligands that populate well-defined conformational states.

Nuclear Magnetic Resonance Insight into the Multiple Glycosaminoglycan Binding Modes of the Link Module from Human TSG-6

Tumor necrosis factor-stimulated gene-6 (TSG-6) is a hyaluronan (HA)-binding protein that is essential for stabilizing and remodeling the extracellular matrix (ECM) during ovulation and inflammatory disease processes such as arthritis. The Link module, one of the domains of TSG-6, is responsible for binding hyaluronan and other glycosaminoglycans found in the ECM. In this study, we used a well-defined chondroitin sulfate (CS) hexasaccharide (ΔC444S) to determine the structure of the Link module, in solution, in its chondroitin sulfate-bound state. A variety of nuclear magnetic resonance techniques were employed, including chemical shift perturbation, residual dipolar couplings (RDCs), nuclear Overhauser effects, spin relaxation measurements, and paramagnetic relaxation enhancements from a spin-labeled analogue of ΔC444S. The binding site for ΔC444S on the Link module overlapped with that of HA. Surprisingly, ΔC444S binding induced dimerization of the Link module (as confirmed by analytical ultracentrifugation), and a second weak binding site that partially overlapped with a previously identified heparin site was detected. A dimer model was generated using chemical shift perturbations and RDCs as restraints in the docking program HADDOCK. We postulate that the molecular cross-linking enhanced by the multiple binding modes of the Link module might be critical for remodeling the ECM during inflammation/ovulation and might contribute to other functions of TSG-6.

Characterization of Aspartate Kinase from Corynebacterium pekinense and the Critical Site of Arg169

Aspartate kinase (AK) is the key enzyme in the biosynthesis of aspartate-derived amino acids. Recombinant AK was efficiently purified and systematically characterized through analysis under optimal conditions combined with steady-state kinetics study. Homogeneous AK was predicted as a decamer with a molecular weight of ~48 kDa and a half-life of 4.5 h. The enzymatic activity was enhanced by ethanol and Ni²⁺. Moreover, steady-state kinetic study confirmed that AK is an allosteric enzyme, and its activity was inhibited by allosteric inhibitors, such as Lys, Met, and Thr. Theoretical results indicated the binding mode of AK and showed that Arg169 is an important residue in substrate binding, catalytic domain, and inhibitor binding. The values of the kinetic parameter V_max of R169 mutants, namely, R169Y, R169P, R169D, and R169H AK, with l-aspartate as the substrate, were 4.71-, 2.25-, 2.57-, and 2.13-fold higher, respectively, than that of the wild-type AK. Furthermore, experimental and theoretical data showed that Arg169 formed a hydrogen bond with Glu92, which functions as the entrance gate. This study provides a basis to develop new enzymes and elucidate the corresponding amino acid production.

Investigation of Carbohydrate Recognition via Computer Simulation

Carbohydrate recognition by proteins, such as lectins and other (bio)molecules, can be essential for many biological functions. Recently, interest has arisen due to potential protein and drug design and future bioengineering applications. A quantitative measurement of carbohydrate-protein interaction is thus important for the full characterization of sugar recognition. We focus on the aspect of utilizing computer simulations and biophysical models to evaluate the strength and specificity of carbohydrate recognition in this review. With increasing computational resources, better algorithms and refined modeling parameters, using state-of-the-art supercomputers to calculate the strength of the interaction between molecules has become increasingly mainstream. We review the current state of this technique and its successful applications for studying protein-sugar interactions in recent years.

The Interaction of Heparin Tetrasaccharides with Chemokine CCL5 Is Modulated by Sulfation Pattern and pH

Interactions between chemokines such as CCL5 and glycosaminoglycans (GAGs) are essential for creating haptotactic gradients to guide the migration of leukocytes into inflammatory sites, and the GAGs that interact with CCL5 with the highest affinity are heparan sulfates/heparin. The interaction between CCL5 and its receptor on monocytes, CCR1, is mediated through residues Arg-17 and -47 in CCL5, which overlap with the GAG-binding (44)RKNR(47) "BBXB" motifs. Here we report that heparin and tetrasaccharide fragments of heparin are able to inhibit CCL5-CCR1 binding, with IC50 values showing strong dependence on the pattern and extent of sulfation. Modeling of the CCL5-tetrasaccharide complexes suggested that interactions between specific sulfate and carboxylate groups of heparin and residues Arg-17 and -47 of the protein are essential for strong inhibition; tetrasaccharides lacking the specific sulfation pattern were found to preferentially bind CCL5 in positions less favorable for inhibition of the interaction with CCR1. Simulations of a 12-mer heparin fragment bound to CCL5 indicated that the oligosaccharide preferred to interact simultaneously with both (44)RKNR(47) motifs in the CCL5 homodimer and engaged residues Arg-47 and -17 from both chains. Direct engagement of these residues by the longer heparin oligosaccharide provides a rationalization for its effectiveness as an inhibitor of CCL5-CCR1 interaction. In this mode, histidine (His-23) may contribute to CCL5-GAG interactions when the pH drops just below neutral, as occurs during inflammation. Additionally, an examination of the contribution of pH to modulating CCL5-heparin interactions suggested a need for careful interpretation of experimental results when experiments are performed under non-physiological conditions.

Impact of sulfation pattern on the conformation and dynamics of sulfated fucan oligosaccharides as revealed by NMR and MD

Sulfated fucans from sea urchin egg jelly express well-defined chemical structures that vary with species. This species specificity regulates the sperm acrosome reaction, a critical step to assure intra-specific fertilization. In addition, these polysaccharides are involved in other biological activities such as anticoagulation. Although sulfation patterns are relevant to the levels of response in both activities, conformation and dynamics of these glycans are also contributing factors. However, data about these features of sulfated fucans are very rare. To address this, we have employed nuclear magnetic resonance experiments combined with molecular dynamics on structurally defined oligosaccharides derived from two sulfated fucans. The results have indicated that the oligosaccharides are flexible in solution. Ring conformation of their composing units displays just the (1)C4 chair configuration. In a particular octasaccharide, composed of two tetrasaccharide sequences, inter-residual hydrogen bonds play a role to decrease dynamics in these repeating units. Conversely, the linking disaccharide [-3)-α-L-Fucp-2(OSO3(-))-(1-3)-α-L-Fucp-4(OCO3(-))-(1-] located right between the two tetrasaccharide units has amplified motions suggested to be promoted by electrostatic repulsion of sulfates on opposite sides of the central glycosidic bond. This conjunction of information about conformation and dynamics of sulfated fucan oligosaccharides provides new insights to explain how these glycans behave free in solution and influenced by sulfation patterns. It may also serve for future studies concerning structure-function relationship of sulfated fucans, especially those involving sea urchin fertilization and anticoagulation.

Predicting the Origins of Anti-Blood Group Antibody Specificity: A Case Study of the ABO A- and B-Antigens

The ABO blood group system is the most important blood type system in human transfusion medicine. Here, we explore the specificity of antibody recognition toward ABO blood group antigens using computational modeling and biolayer interferometry. Automated docking and molecular dynamics simulations were used to explore the origin of the specificity of an anti-blood group A antibody variable fragment (Fv AC1001). The analysis predicts a number of Fv-antigen interactions that contribute to affinity, including a hydrogen bond between a His^L49 and the carbonyl moiety of the GalNAc in antigen A. This interaction was consistent with the dependence of affinity on pH, as measured experimentally; at lower pH there is an increase in binding affinity. Binding energy calculations provide unique insight into the origin of interaction energies at a per-residue level in both the scFv and the trisaccharide antigen. The calculations indicate that while the antibody can accommodate both blood group A and B antigens in its combining site, the A antigen is preferred by 4 kcal/mol, consistent with the lack of binding observed for the B antigen.

Recent advances in employing molecular modelling to determine the specificity of glycan-binding proteins

Impressive improvements in docking performance can be achieved by applying energy bonuses to poses in which glycan hydroxyl groups occupy positions otherwise preferred by bound waters. In addition, inclusion of glycosidic conformational energies allows unlikely glycan conformations to be appropriately penalized. A method for predicting the binding specificity of glycan-binding proteins has been developed, which is based on grafting glycan branches onto a minimal binding determinant in the binding site. Grafting can be used either to screen virtual libraries of glycans, such as the known glycome, or to identify docked poses of minimal binding determinants that are consistent with specificity data. The reviewed advances allow accurate modelling of carbohydrate-protein 3D co-complexes, but challenges remain in ranking the affinity of congeners.