Involvement of water in carbohydrate-protein binding

Halogen Bonding on Water─A Drop in the Ocean?

Halogen bonding is a valuable interaction in drug design, offering an unconventional way to influence affinity and selectivity by leveraging the halogen atoms' ability to form directional bonds. The present study evaluates halogen-water interactions within protein binding sites, demonstrating that targeting a water molecule via halogen bonding can in specific cases contribute beneficially to ligand binding. In solving and examining the crystal structure of 2-cyclopentyl-7-iodo-1H-indole-3-carbonitrile bound to DYRK1a kinase, we identified a notable iodine-water interaction, where water accepts a halogen bond with good geometric and energetic features. This starting point triggered further investigations into the prevalence of such interactions across various halogen-bearing ligands (chlorine, bromine, iodine) in the PDB. Using QM calculations (MP2/TZVPP), we highlight the versatility and potential benefits of such halogen-water interactions, particularly when the water molecule is a stable part of the binding site's structured environment. While the interaction energies with water are lower compared to other typical halogen bond acceptors, we deem this different binding strength essential for reducing desolvation costs. We suggest that "interstitial" water molecules, as stable parts of the binding site engaging in multiple strong interactions, could be prime targets for halogen bonding. Further systematic studies, combining high-resolution crystal structures and quantum chemistry, are required to scrutinize whether halogen bonding on water is more than a "drop in the ocean".

Evaluation of phyto-gallic acid as a potential inhibitor of Staphylococcus aureus efflux pump mediated tetracycline resistance: an in-vitro and in-silico study

Plants contain many bioactive compounds with potent antibacterial and efflux pump inhibitory activity (EPI). In this study, gallic acid extracted from pomegranate molasses by analytical HPLC holds promise as an EPI drug for Staphylococcus aureus mediated tetracycline resistance, it lowered the bacterial resistance and reversed the mechanism via tet family efflux pump, using molecular technique and in-silico molecular docking analysis. Extracted gallic acid combined with tetracycline demonstrated a significant decrease in the minimal inhibitory concentration MIC compared to its single activity. Similarly, little growth and lower fluorescence of S. aureus were observed on ethidium bromide (2.5 mg/mL) agar plates, indicating a reversible efflux pump mechanism and a potent EPI activity. Molecular docking demonstrated a promising affinity binding energy between gallic acid and tet efflux genes, opening a new baseline in bacterial infection treatment. PCR for tetK and Qac A/B genes failed to show any relation between tet genes and gallic acid.

Outstanding Properties of the Hydration Shell around β-d-Glucose: A Computational Study

Ab initio molecular dynamics (AIMD) simulations have been performed on aqueous solutions of four simple sugars, α-d-glucose, β-d-glucose, α-d-mannose, and α-d-galactose. Hydrogen-bonding (HB) properties, such as the number of donor- and acceptor-type HB-s, and the lengths and strengths of hydrogen bonds between sugar and water molecules, have been determined. Related electronic properties, such as the dipole moments of water molecules and partial charges of the sugar O atoms, have also been calculated. The hydrophilic and hydrophobic shells were characterized by means of spatial distribution functions. β-d-Glucose was found to form the highest number of hydrophilic and the smallest number of hydrophobic connections to neighboring water molecules. The average sugar-water H-bond length was the shortest for β-d-glucose, which suggests that these are the strongest such H-bonds. Furthermore, β-d-glucose appears to stand out in terms of the symmetry properties of both its hydrophilic and hydrophobic hydration shells. In summary, in all aspects considered here, there seems to be a correlation between the distinct characteristics of β-d-glucose reported here and its outstanding solubility in water. Admittedly, our findings represent only some of the important factors that influence the solubility.

HydraProt: A New Deep Learning Tool for Fast and Accurate Prediction of Water Molecule Positions for Protein Structures

Water molecules are integral to the structural stability of proteins and vital for facilitating molecular interactions. However, accurately predicting their precise position around protein structures remains a significant challenge, making it a vibrant research area. In this paper, we introduce HydraProt (deep Hydration of Proteins), a novel methodology for predicting precise positions of water molecule oxygen atoms around protein structures, leveraging two interconnected deep learning architectures: a 3D U-net and a Multi-Layer Perceptron (MLP). Our approach starts by introducing a coarse voxel-based representation of the protein, which allows for rapid sampling of candidate water positions via the 3D U-net. These water positions are then assessed by embedding the water-protein relationship in the Euclidean space by means of an MLP. Finally, a postprocessing step is applied to further refine the MLP predictions. HydraProt surpasses existing state-of-the-art approaches in terms of precision and recall and has been validated on large data sets of protein structures. Notably, our method offers rapid inference runtime and should constitute the method of choice for protein structure studies and drug discovery applications. Our pretrained models, data, and the source code required to reproduce these results are accessible at https://github.com/azamanos/HydraProt.

Kinetics and thermodynamics of thermal inactivation for recombinant Escherichia coli cellulases, cel12B, cel8C, and polygalacturonase, peh28; biocatalysts for biofuel precursor production

Lignocellulosic biomass conversion using cellulases/polygalacturonases is a process that can be progressively influenced by several determinants involved in cellulose microfibril degradation. This article focuses on the kinetics and thermodynamics of thermal inactivation of recombinant Escherichia coli cellulases, cel12B, cel8C and a polygalacturonase, peh 28, derived from Pectobacterium carotovorum sub sp. carotovorum. Several consensus motifs conferring the enzymes' thermal stability in both cel12B and peh28 model structures have been detailed earlier, which were confirmed for the three enzymes through the current study of their thermal inactivation profiles over the 20-80°C range using the respective activities on carboxymethylcellulose and polygalacturonic acid. Kinetic constants and half-lives of thermal inactivation, inactivation energy, plus inactivation entropies, enthalpies and Gibbs free energies, revealed high stability, less conformational change and protein unfolding for cel12B and peh28 due to thermal denaturation compared to cel8C. The apparent thermal stability of peh28 and cel12B, along with their hydrolytic efficiency on a lignocellulosic biomass conversion as reported previously, makes these enzymes candidates for various industrial applications. Analysis of the Gibbs free energy values suggests that the thermal stabilities of cel12B and peh28 are entropy-controlled over the tested temperature range.

Effect of hydroxylysine-O-glycosylation on the structure of type I collagen molecule: A computational study

Collagen undergoes many types of post-translational modifications (PTMs), including intracellular modifications and extracellular modifications. Among these PTMs, glycosylation of hydroxylysine (Hyl) is the most complicated. Experimental studies demonstrated that this PTM ceases once the collagen triple helix is formed and that Hyl-O-glycosylation modulates collagen fibrillogenesis. However, the underlying atomic-level mechanisms of these phenomena remain unclear. In this study, we first adapted the force field parameters for O-linkages between Hyl and carbohydrates and then investigated the influence of Hyl-O-glycosylation on the structure of type I collagen molecule, by performing comprehensive molecular dynamic simulations in explicit solvent of collagen molecule segment with and without the glycosylation of Hyl. Data analysis demonstrated that (i) collagen triple helices remain in a triple-helical structure upon glycosylation of Hyl; (ii) glycosylation of Hyl modulates the peptide backbone conformation and their solvation environment in the vicinity and (iii) the attached sugars are arranged such that their hydrophilic faces are well exposed to the solvent, while their hydrophobic faces point towards the hydrophobic portions of collagen. The adapted force field parameters for O-linkages between Hyl and carbohydrates will aid future computational studies on proteins with Hyl-O-glycosylation. In addition, this work, for the first time, presents the detailed effect of Hyl-O-glycosylation on the structure of human type I collagen at the atomic level, which may provide insights into the design and manufacture of collagenous biomaterials and the development of biomedical therapies for collagen-related diseases.

Targeting prolyl tripeptidyl peptidase from Porphyromonas gingivalis with the bioactive compounds from Rosmarinus officinalis

Background

Complications in periodontitis and other systemic infections related to Porphyromonas gingivalis poses a serious impediment in the treatment process. This leads to the search of novel target proteins to develop newer drugs against P. gingivalis. Prolyl tripeptidyl peptidase (ptp-A) seem to be a vital protein in P. gingivalis virulence and can be a good target for the novel natural bioactive compounds.

Objectives

To explore the inhibitory potential of Rosmarinus officinalis biocompounds against the ptp-A of P. gingivalis.

Methods

Three-dimensional structure of ptp-A was retrieved from the Protein Data Bank with further optimization of both the protein and ligands. In silico inhibitory potential of the selected ligands against ptp-A was done by AutoDock 2.0 and was visualized with Biovia discovery studio visualizing tool with the assessment of the molecular properties of the ligands against ptp-A by molinspiration calculations and drug likeliness.

Results

High ptp-A inhibitory effect was observed using rosmarinic acid and luteolin with a bonding energy of −9.81 kcal/mol with 10 hydrogen bond interactions and −9.99 kcal/mol with 7 hydrogen bond interactions, respectively. Carnosic acid and p-coumaric acid showed a binding energy of −7.14 kcal/mol and −6.34 kcal/mol, respectively, with 5 hydrogen bond interactions. Molinspiration assessments showed R. officinalis compounds as the best drug candidates with the topological polar surface area scores <140 Å toward the best oral bioavailability.

Conclusion

The carnosic acid, rosmarinic acid, p-coumaric acid, and luteolin from R. officinalis seem to possess a promising inhibitory effect against ptp-A of Candida albicans suggesting ptp-A as the best target to combat P. gingivalis with further in vivo validation.

The role of hydration effects in 5-fluorouridine binding to SOD1: insight from a new 3D-RISM-KH based protocol for including structural water in docking simulations

Misfolded Cu/Zn superoxide dismutase enzyme (SOD1) shows prion-like propagation in neuronal cells leading to neurotoxic aggregates that are implicated in amyotrophic lateral sclerosis (ALS). Tryptophan-32 (W32) in SOD1 is part of a potential site for templated conversion of wild type SOD1. This W32 binding site is located on a convex, solvent exposed surface of the SOD1 suggesting that hydration effects can play an important role in ligand recognition and binding. A recent X-ray crystal structure has revealed that 5-Fluorouridine (5-FUrd) binds at the W32 binding site and can act as a pharmacophore scaffold for the development of anti-ALS drugs. In this study, a new protocol is developed to account for structural (non-displaceable) water molecules in docking simulations and successfully applied to predict the correct docked conformation binding modes of 5-FUrd at the W32 binding site. The docked configuration is within 0.58 Å (RMSD) of the observed configuration. The docking protocol involved calculating a hydration structure around SOD1 using molecular theory of solvation (3D-RISM-KH, 3D-Reference Interaction Site Model-Kovalenko-Hirata) whereby, non-displaceable water molecules are identified for docking simulations. This protocol was also used to analyze the hydrated structure of the W32 binding site and to explain the role of solvation in ligand recognition and binding to SOD1. Structural water molecules mediate hydrogen bonds between 5-FUrd and the receptor, and create an environment favoring optimal placement of 5-FUrd in the W32 binding site.

A computational study targeting the mutated L321F of ERG11 gene in C. albicans, associated with fluconazole resistance with bioactive compounds from Acacianilotica

BACKGROUND

Emergence of fluconazole resistance mediated by L321F mutation in the ERG11 gene poses a serious impediment in the candidal treatment process. This leads to the search of novel target proteins to develop newer drugs against fluconazole resistant C. albicans. The present investigation is thus aimed to explore the inhibitory potential of bioactive compounds from A. nilotica against the wild and mutated ERG11 gene of C. albicans.

MATERIALS AND METHODS

Homology modelling of the wild and mutated ERG11 target gene was done by Modeller 9.2, with SDM I-mutant form of ERG11 together with SAVES server and Ramachandran plot validation. 2D and 3D structures of the bioactive compounds from A. nilotica were optimized by ACD chemsketch. Molinspirational assessments for the molecular properties of the ligands and their drug likeliness were estimated. In silico inhibitory potential of the selected ligands against wild and mutated ERG11 was done by AutoDock 2.0 and was visualized with Accelrys discovery studio tool.

RESULTS

Apigenin proved to be the best candidate to target mutant ERG11 with a binding energy of -8.33 Kcal/mol followed by catechin with six hydrogen bonds with more drug likeliness. Molinspiration assessments showed zero violations for all the bioactive compounds from A. nilotica and the TPSA scores of the ligands showed the values<140Å towards the best oral bio-availability.

CONCLUSION

The findings of the study emphasize that kaempferol, apigenin and catechin from A. nilotica seem to possess a promising inhibitory effect against the wild and mutated ERG11 of C. albicans suggesting ERG11 as the best target to combat fluconazole resistant C. albicans with further in vivo validation targeting the same.

Distribution of toxicity values across different species and modes of action of pesticides from PESTIMEP and PPDB databases

The continuous use of compounds contained in commodities such as processed food, medicines, and pesticides, demands safety measures, in particular, for those in direct contact with humans and the environment. Safety measures have evolved and regulations are now in place around the globe. In the case of pesticides, attempts have been made to use toxicological data to inform of potentially harmful compounds either across species, on different routes of exposure, or entirely new chemicals. The generation of models, based on statistical and molecular modeling studies, allows for such predictions. However, the use of these models is framed by the available data, the experimental errors, the complexity of the measurement, and the available computational algorithms, among other factors. In this work, we present the methodologies used for extrapolation across different species and routes of administration and show the appropriateness of developing predictive models of pesticides based on their type and mode of action. The analyses include comparisons based on structural characteristics and physicochemical properties. Whenever possible, the scope and limitations of the methodologies are discussed. We expect that this work will serve as a useful introductory guide of the tools employed in the toxicity assessment of agrochemical compounds.

Exploring interfacial water trapping in protein-ligand complexes with multithermal titration calorimetry

In this work, we examine the hypothesis about how trapped water molecules at the interface between triosephosphate isomerase (TIM) and either of two phosphorylated inhibitors, 2-phosphoglycolate (2PG) or phosphoglycolohydroxamate (PGH), can explain the anomalous highly negative binding heat capacities (ΔC_p,b) of both complexes, TIM-2PG and TIM-PGH. We performed fluorimetric titrations of the enzyme with PGH inhibitor under osmotic stress conditions, using various concentrations of either osmolyte: sucrose, ethylene glycol or glycine betaine. We also analyze the binding processes under various stressor concentrations using a novel calorimetric methodology that allows ΔC_p,b determinations in single experiments: Multithermal Titration Calorimetry. The binding constant of the TIM-PGH complex decreased gradually with the concentration of all osmolytes, but at diverse extents depending on the osmolyte nature. According to the osmotic stress theory, this decrease indicates that the number of water molecules associated with the enzyme increases with inhibitor binding, i.e. some solvent molecules became trapped. Additionally, the binding heat capacities became less negative at higher osmolyte concentrations, their final values depending on the osmolyte. These effects were also observed in the TIM-2PG complex using sucrose as stressor. Our results strongly suggest that some water molecules became immobilized when the TIM-inhibitor complexes were formed. A computational analysis of the hydration state of the binding site of TIM in both its free state and its complexed form with 2PG or PGH, based on molecular dynamics (MD) simulations in explicit solvent, showed that the binding site effectively immobilized additional water molecules after binding these inhibitors.

Effective lead optimization targeting the displacement of bridging receptor–ligand water molecules

Enhancing the binding affinities of ligands by means of lead modifications that displace bridging water molecules at protein–ligand interfaces is an important and widely studied lead optimization strategy.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Testing the CP-correction procedure with different DFT methods on H-bonding complexes of κ-carrabiose with water molecules

Interaction of water molecules with κ-carrabiose disaccharide, within three H-bonding complexes, was investigated. Particular interest was focused on the way with which the BSSE correction has to be performed. Two strategies were used, either performing BSSE correction during or after optimization. For this aim, several DFT-functionals (hybrid GGA and hybrid meta-GGA) and 6-31 + G* basis set were considered. The results demonstrated the uselessness of including of BSSE-CP correction during optimization for all complexes. From a structural point of view, a proper H-bonding description was obtained using the PBE0 functional for all complexes. The basis set effect on the BSSE using B3LYP functional was also investigated. The reliability of B3LYP/6-31 + G** and B3LYP/6-31++G** models for the complexes involving one or two water molecules was reported while the use of B3LYP/6-311 + G** or B3LYP/6-311++G** levels was shown to be more appropriate for larger complexes equivalent to that involving three water molecules. CP-corrected interaction energies were demonstrated to be closer to CBS-4 M interaction energies than the uncorrected ones. Graphical abstract Functional and basis set effects on BSSE.

Calculations of the absolute binding free energies for Ralstonia solanacearum lectins bound with methyl-α-l-fucoside at molecular mechanical and quantum mechanical/molecular mechanical levels

In this work, both a molecular mechanical (MM) method and a hybrid quantum mechanical/molecular mechanical (QM/MM) method have been applied in the study of the binding affinities of methyl-α-l-fucoside to Ralstonia solanacearum lectins.

Photoelectron spectroscopy and density functional theory studies of (fructose + (H2O)n)− (n = 1–5) anionic clusters

(Fructose + (H₂O)_n)⁻ (n = 1–5) cluster anions mainly exist as open chain structures. Some cyclic structures of (fructose + (H₂O)_n)⁻ (n = 3, 4) are present too.

An insight into the thermodynamic characteristics of human thrombopoietin complexation with TN1 antibody

Human thrombopoietin (hTPO) primarily stimulates megakaryocytopoiesis and platelet production and is neutralized by the mouse TN1 antibody. The thermodynamic characteristics of TN1 antibody-hTPO complexation were analyzed by isothermal titration calorimetry (ITC) using an antigen-binding fragment (Fab) derived from the TN1 antibody (TN1-Fab). To clarify the mechanism by which hTPO is recognized by TN1-Fab the conformation of free TN1-Fab was determined to a resolution of 2.0 Å using X-ray crystallography and compared with the hTPO-bound form of TN1-Fab determined by a previous study. This structural comparison revealed that the conformation of TN1-Fab does not substantially change after hTPO binding and a set of 15 water molecules is released from the antigen-binding site (paratope) of TN1-Fab upon hTPO complexation. Interestingly, the heat capacity change (ΔCp) measured by ITC (-1.52 ± 0.05 kJ mol(-1)  K(-1) ) differed significantly from calculations based upon the X-ray structure data of the hTPO-bound and unbound forms of TN1-Fab (-1.02 ∼ 0.25 kJ mol(-1)  K(-1) ) suggesting that hTPO undergoes an induced-fit conformational change combined with significant desolvation upon TN1-Fab binding. The results shed light on the structural biology associated with neutralizing antibody recognition.

How Does Solvation Affect the Binding of Hydrophilic Amino Saccharides to Cucurbit[7]uril with Exceptional Anomeric Selectivity?

Cucurbit[7]uril (CB[7]) is known to bind strongly to hydrophilic amino saccharide guests with exceptional α-anomer selectivities under aqueous conditions. Single-crystal X-ray crystallography and computational methods were used to elucidate the reason behind this interesting phenomenon. The crystal structures of protonated galactosamine (GalN) and glucosamine (GluN) complexes confirm the inclusion of α anomers inside CB[7] and disclose the details of the host-guest binding. Whereas computed gas-phase structures agree with these crystal structures, gas-phase binding free energies show preferences for the β-anomer complexes over their α counterparts, in striking contrast to the experimental results under aqueous conditions. However, when the solvation effect is considered, the binding structures drastically change and the preference for the α anomers is recovered. The α anomers also tend to bind more tightly and leave less space in the CB[7] cavity toward inclusion of only one water molecule, whereas loosely bound β anomers leave more space toward accommodating two water molecules, with markedly different hydrogen-bonding natures. Surprisingly, entropy seems to contribute significantly to both anomeric discrimination and binding. This suggests that of all the driving factors for the strong complexation of the hydrophilic amino saccharide guests, water mediation plays a crucial role in the anomer discrimination.

Dowser++, a new method of hydrating protein structures

A new method of hydrating protein structures, which we call Dowser++, is presented. The method is based on a semi-empirical modification of a popular program for protein hydration Dowser, and the usage of protocols AutoDock Vina, and WaterDock. The positions of water molecules predicted by Dowser++ were compared with experimental data for a set of 14 high-resolution crystal structures of oligopeptide-binding protein (OppA) containing a large number of resolved internal water molecules, as well as for the D- and K-channels of cytochrome c oxidase, and the recent data on PSII. Comparison is also made with the predictions of the original Dowser, and its improved version, Dowser+, described in our previous publication. We also present a model for quantitative estimation of the quality of water molecules placement made by a program, which includes an assumption of possible false negative data from the crystallographic analysis. The comparison of predictions made by Dowser++, Dowser and Dowser+ demonstrates significant improvement of predictive power of the new method. Proteins 2016; 84:1347-1357. © 2016 Wiley Periodicals, Inc.

Glucose induced variation of water structure from temperature dependent near infrared spectra

The different effects of glucose on water species provide evidence to explain the bioprotective function of carbohydrates in aqueous solutions.

Beyond Membrane Protein Structure: Drug Discovery, Dynamics and Difficulties

Most of the previous content of this book has focused on obtaining the structures of membrane proteins. In this chapter we explore how those structures can be further used in two key ways. The first is their use in structure based drug design (SBDD) and the second is how they can be used to extend our understanding of their functional activity via the use of molecular dynamics. Both aspects now heavily rely on computations. This area is vast, and alas, too large to consider in depth in a single book chapter. Thus where appropriate we have referred the reader to recent reviews for deeper assessment of the field. We discuss progress via the use of examples from two main drug target areas; G-protein coupled receptors (GPCRs) and ion channels. We end with a discussion of some of the main challenges in the area.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Acα(2-3)Gal: a study by in silico mutations and molecular dynamics simulations

Wheat germ agglutinin (WGA) is a plant lectin, which specifically recognizes the sugars NeuNAc and GlcNAc. Mutated WGA with enhanced binding specificity can be used as biomarkers for cancer. In silico mutations are performed at the active site of WGA to enhance the binding specificity towards sialylglycans, and molecular dynamics simulations of 20 ns are carried out for wild type and mutated WGAs (WGA1, WGA2, and WGA3) in complex with sialylgalactose to examine the change in binding specificity. MD simulations reveal the change in binding specificity of wild type and mutated WGAs towards sialylgalactose and bound conformational flexibility of sialylgalactose. The mutated polar amino acid residues Asn114 (S114N), Lys118 (G118K), and Arg118 (G118R) make direct and water mediated hydrogen bonds and hydrophobic interactions with sialylgalactose. An analysis of possible hydrogen bonds, hydrophobic interactions, total pair wise interaction energy between active site residues and sialylgalactose and MM-PBSA free energy calculation reveals the plausible binding modes and the role of water in stabilizing different binding modes. An interesting observation is that the binding specificity of mutated WGAs (cyborg lectin) towards sialylgalactose is found to be higher in double point mutation (WGA3). One of the substituted residues Arg118 plays a crucial role in sugar binding. Based on the interactions and energy calculations, it is concluded that the order of binding specificity of WGAs towards sialylgalactose is WGA3 > WGA1 > WGA2 > WGA. On comparing with the wild type, double point mutated WGA (WGA3) exhibits increased specificity towards sialylgalactose, and thus, it can be effectively used in targeted drug delivery and as biological cell marker in cancer therapeutics.

Calculating binding free energies for protein-carbohydrate complexes

A variety of computational techniques may be applied to compute theoretical binding free energies for protein-carbohydrate complexes. Elucidation of the intermolecular interactions, as well as the thermodynamic effects, that contribute to the relative strength of receptor binding can shed light on biomolecular recognition, and the resulting initiation or inhibition of a biological process. Three types of free energy methods are discussed here, including MM-PB/GBSA, thermodynamic integration, and a non-equilibrium alternative utilizing SMD. Throughout this chapter, the well-known concanavalin A lectin is employed as a model system to demonstrate the application of these methods to the special case of carbohydrate binding.

Analysis of factors influencing hydration site prediction based on molecular dynamics simulations

Water contributes significantly to the binding of small molecules to proteins in biochemical systems. Molecular dynamics (MD) simulation based programs such as WaterMap and WATsite have been used to probe the locations and thermodynamic properties of hydration sites at the surface or in the binding site of proteins generating important information for structure-based drug design. However, questions associated with the influence of the simulation protocol on hydration site analysis remain. In this study, we use WATsite to investigate the influence of factors such as simulation length and variations in initial protein conformations on hydration site prediction. We find that 4 ns MD simulation is appropriate to obtain a reliable prediction of the locations and thermodynamic properties of hydration sites. In addition, hydration site prediction can be largely affected by the initial protein conformations used for MD simulations. Here, we provide a first quantification of this effect and further indicate that similar conformations of binding site residues (RMSD < 0.5 Å) are required to obtain consistent hydration site predictions.

Dipole Moment and Binding Energy of Water in Proteins from Crystallographic Analysis

The energetics of water molecules in proteins is studied using the water placement software Dowser. We compared the water position predictions for 14 high-resolution crystal structures of oligopeptide-binding protein (OppA) containing a large number of resolved internal water molecules. From the analysis of the outputs of Dowser with variable parameters and comparison with experimental X-ray data, we derived an estimate of the average dipole moment of water molecules located in the internal cavities of the protein and their binding energies. The water parameters thus obtained from the experimental data are then analyzed within the framework of charge-scaling theory developed recently by this group; the parameters are shown to be in good agreement with the predictions that the theory makes for the dipole moment in a protein environment. The water dipole in the protein environment is found to be much different from that in the bulk and in such models as SPC or TIPnP. The role of charge scaling due to electronic polarizability of the protein is discussed.

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.

Glycopolymer probes of signal transduction

Glycans are key participants in biological processes ranging from reproduction to cellular communication to infection. Revealing glycan roles and the underlying molecular mechanisms by which glycans manifest their function requires access to glycan derivatives that vary systematically. To this end, glycopolymers (polymers bearing pendant carbohydrates) have emerged as valuable glycan analogs. Because glycopolymers can readily be synthesized, their overall shape can be varied, and they can be altered systematically to dissect the structural features that underpin their activities. This review provides examples in which glycopolymers have been used to effect carbohydrate-mediated signal transduction. Our objective is to illustrate how these powerful tools can reveal the molecular mechanisms that underlie carbohydrate-mediated signal transduction.

Restricted dynamics of water around a protein-carbohydrate complex: computer simulation studies

Water-mediated protein-carbohydrate interaction is a complex phenomenon responsible for different biological processes in cellular environment. One of the unexplored but important issues in this area is the role played by water during the recognition process and also in controlling the microscopic properties of the complex. In this study, we have carried out atomistic molecular dynamics simulations of a protein-carbohydrate complex formed between the hyaluronan binding domain of the murine Cd44 protein and the octasaccharide hyaluronan in explicit water. Efforts have been made to explore the heterogeneous influence of the complex on the dynamic properties of water present in different regions around it. It is revealed from our analyses that the heterogeneous dynamics of water around the complex are coupled with differential time scales of formation and breaking of hydrogen bonds at the interface. Presence of a highly rigid thin layer of motionally restricted water molecules bridging the protein and the carbohydrate in the common region of the complex has been identified. Such water molecules are expected to play a crucial role in controlling properties of the complex. Importantly, it is demonstrated that the formation of the protein-carbohydrate complex affects the transverse and longitudinal degrees of freedom of the interfacial water molecules in a heterogeneous manner.

Quantifying the role of water in protein-carbohydrate interactions

Water-mediated interactions play a key role in carbohydrate-lectin binding, where the interactions involve a conserved water that is separated from the bulk solvent and present a bridge between the side chains of the protein and the carbohydrate ligand. To apply quantum mechanical methods to examine the role of conserved waters, we present an analysis in which the relevant carbohydrate atoms are modeled by methanol, and in which the protein is replaced by a limited number of amino acid side chains. Clusters containing a conserved water and a representative amino acid fragment were also examined to determine the influence of amino acid side chains on interaction energies. To quantify the differential binding energies of methanol versus water, quantum mechanical calculations were performed at the B3LYP/6-311++G(3df,3pd)//B3LYP/6-31+G(d) level in which either a methanol molecule was bound to the conserved water (liganded state) or in which a water molecule replaces the methanol (unliganded state). Not surprisingly, the binding of a water to clusters containing charged amino acid side chains was more favorable by 1.55 to 7.23 kcal/mol than that for the binding of a water to the corresponding pure water clusters. In contrast, the binding energy of water to clusters containing polar-uncharged amino acid side chains ranged from 4.35 kcal/mol less favorable to 4.72 kcal/mol more favorable than for binding to the analogous pure water clusters. The overall trend for the binding of methanol versus water, in any of the clusters, favored methanol by an average value of 1.05 kcal/mol. To extend these studies to a complex between a protein (Concanavalin A) and its carbohydrate ligand, a cluster was examined that contained the side chains of three key amino acids, namely asparagine, aspartate, and arginine, as well as a key water molecule, arranged as in the X-ray diffraction structure of Con A. Again, using methanol as a model for the endogenous carbohydrate ligand, energies of -5.94 kcal/mol and -5.70 kcal/mol were obtained for the binding of methanol and water, respectively, to the Con A-water cluster. The extent to which cooperativity enhanced the binding energies has been quantified in terms of nonadditive three-body contributions. In general, the binding of water or methanol to neutral dimers formed cooperative clusters; in contrast, the cooperativity in charged clusters depended on the overall geometry as well as the charge.

Synthesis of 4" manipulated Lewis X trisaccharide analogues

Three analogues of the Le^x trisaccharide antigen (β-D-Galp(1→4)[α-L-Fucp(1→3)]-D-GlcNAcp) in which the galactosyl residue is modified at O-4 as a methyloxy, deoxychloro or deoxyfluoro, were synthesized. We first report the preparation of the modified 4-OMe, 4-Cl and 4-F trichloroacetimidate galactosyl donors and then report their use in the glycosylation of an N-acetylglucosamine glycosyl acceptor. Thus, we observed that the reactivity of these donors towards the BF₃·OEt₂-promoted glycosylation at O-4 of the N-acetylglucosamine glycosyl acceptors followed the ranking 4-F > 4-OAc ≈ 4-OMe > 4-Cl. The resulting disaccharides were deprotected at O-3 of the glucosamine residue and fucosylated, giving access to the desired protected Le^x analogues. One-step global deprotection (Na/NH₃) of the protected 4”-methoxy analogue, and two-step deprotections (removal of a p-methoxybenzyl with DDQ, then Zemplén deacylation) of the 4”-deoxychloro and 4”-deoxyfluoro protected Le^x analogues gave the desired compounds in good yields.

Conformational flexibility of a protein-carbohydrate complex and the structure and ordering of surrounding water

Protein-carbohydrate non-covalent interactions are important to understand various biological processes in living organisms. One of the important issues in protein-carbohydrate binding is how the protein identifies the target carbohydrate and recognizes its conformational features. Surrounding water molecules are expected to play a critical role not only in mediating the recognition process but also in maintaining the structure of the complex. We carried out atomistic molecular dynamics (MD) simulations of an aqueous solution of the protein-carbohydrate complex formed between the hyaluronan binding domain (HABD) of the murine Cd44 protein and the octasaccharide hyaluronan (HA(8)). The conformational flexibilities of the protein and the carbohydrate, and the microscopic structure and ordering of water molecules around them in the complexed form have been explored. It is revealed that the formation of the complex is associated with significant immobilization of the monosaccharide units of the carbohydrate moiety that are involved in binding. Further, reduction in water densities around the binding residues of the two molecules in the complex with respect to their free forms clearly demonstrated that the recognition between the protein and the carbohydrate is facilitated by removal of a fraction of water molecules from regions around the binding domains.