The Role of Bridging Water and Hydrogen Bonding as Key Determinants of Noncovalent Protein-Carbohydrate Recognition

Conformational preferences of heparan sulfate to recognize the CXCL8 dimer in aqueous medium: degree of sulfation and hydrogen bonds

The sulfation pattern and epimerization of the long-chain sulfated polysaccharide heparan sulfate (HS) cause structural diversity and regulate various physiological and pathological processes when binding with proteins. In this work, we performed a series of molecular dynamics simulations of three variants of the octadecasaccharide HS with varying sulfation positions in aqueous medium in their free forms and in the presence of the chemokine CXCL8 dimer. The free energy of binding depicts the sulfation at the 6-O position of GlcNAc (HS6S), and both 3-O and 6-O positions of GlcNAc (HS3S6S) of HS variants are more likely to bind with the CXCL8 dimer than the triply sulfated HS2S3S6S, which is sulfated at the 2-O position of GlcUA additionally along with 3-O and 6-O positions of GlcNAc. Binding between HS and CXCL8 was driven by electrostatic and van der Waals interactions predominantly regardless of the sulfation pattern; however, unfavorable entropic contribution suppressed the interaction between HS and CXCL8. The contribution of different amino acid residues to the binding energetics suggested that basic amino acids line up the binding site of CXCL8. This study further acknowledges the role of interfacial water that is structured and bound with HS through hydrogen bonds, exhibiting differential hydrogen bond relaxation dynamics compared to when the HS molecules are free. Moreover, this study identifies that with the increase in sulfation, the HS-water hydrogen bond relaxation occurs faster with the complexation, while the reverse trend is followed in their free forms. Significant structural adaptation of the different sulfated HS molecules, as verified from the free energy landscapes generated from various reaction coordinates, root-mean-square-deviations, end-to-end distances, including ring pucker angles, dihedral flexibility, and the high conformational entropy cost arising from the glycosidic bonds, suggests that the different sulfated variants of HS undergo significant structural transformation to bind with CXCL8. The presence of a CXCL8 dimer imposes the bound forms of HS to adopt non-linear structures with skew-boat conformations. The atomistic details of the study would help in understanding the selectivity and conformational diversity, as well as the role of solvents in the recognition of CXCL8 by different sulfated variants of HS molecules.

Infrared and vibrational circular dichroism spectra of methyl β-D-glucopyranose in water: The application of the quantum cluster growth and clusters-in-a-liquid solvation models

The infrared (IR) and vibrational circular dichroism (VCD) spectra of methyl β-D-glucopyranose in water were measured. Both implicit and explicit solvation models were utilized to explain the observed spectra. The vast body of existing experimental and theoretical data suggested that about eight explicit water molecules are needed to account for the solvent effects, supported by the current Quantum Cluster Growth (QCG) analysis. Extensive manual and systematic conformational searches of the molecular target and its water clusters were carried out by using a recently developed conformational searching tool, conformer-rotamer ensemble sampling tool (CREST), and the microsolvation model in the associated QCG code. The Boltzmann averaged IR and VCD spectra of the methyl β-D-glucopyranose-(water)n (n = 8) conformers in the PCM of water provide better agreement with the experimental ones than those with n = 0, 1, and 2. The explicit solvation with eight water molecules was shown to greatly modify the conformational preference of methyl β-D-glucopyranose from its monomeric form. Further analyses show that the result is consistent with the existence of long-lived methyl β-D-glucopyranose monohydrates with the additional explicit water effects being accounted for with the quantum mechanical treatment of the other seven close-by water molecules in the PCM of water.

Accounting for Solvation Correlation Effects on the Thermodynamics of Water Networks in Protein Cavities

Macromolecular recognition and ligand binding are at the core of biological function and drug discovery efforts. Water molecules play a significant role in mediating the protein-ligand interaction, acting as more than just the surrounding medium by affecting the thermodynamics and thus the outcome of the binding process. As individual water contributions are impossible to measure experimentally, a range of computational methods have emerged to identify hydration sites in protein pockets and characterize their energetic contributions for drug discovery applications. Even though several methods model solvation effects explicitly, they focus on determining the stability of specific water sites independently and neglect solvation correlation effects upon replacement of clusters of water molecules, which typically happens in hit-to-lead optimization. In this work, we rigorously determine the conjoint effects of replacing all combinations of water molecules in protein binding pockets through the use of the RE-EDS multistate free-energy method, which combines Hamiltonian replica exchange (RE) and enveloping distribution sampling (EDS). Applications on the small bovine pancreatic trypsin inhibitor and four proteins of the bromodomain family illustrate the extent of solvation correlation effects on water thermodynamics, with the favorability of replacement of the water sites by pharmacophore probes highly dependent on the composition of the water network and the pocket environment. Given the ubiquity of water networks in biologically relevant protein targets, we believe our approach can be helpful for computer-aided drug discovery by providing a pocket-specific and a priori systematic consideration of solvation effects on ligand binding and selectivity.

Development and Evaluation of GlycanDock: A Protein-Glycoligand Docking Refinement Algorithm in Rosetta

Carbohydrate chains are ubiquitous in the complex molecular processes of life. These highly diverse chains are recognized by a variety of protein receptors, enabling glycans to regulate many biological functions. High-resolution structures of protein-glycoligand complexes reveal the atomic details necessary to understand this level of molecular recognition and inform application-focused scientific and engineering pursuits. When experimental challenges hinder high-throughput determination of quality structures, computational tools can, in principle, fill the gap. In this work, we introduce GlycanDock, a residue-centric protein-glycoligand docking refinement algorithm developed within the Rosetta macromolecular modeling and design software suite. We performed a benchmark docking assessment using a set of 109 experimentally determined protein-glycoligand complexes as well as 62 unbound protein structures. The GlycanDock algorithm can sample and discriminate among protein-glycoligand models of native-like structural accuracy with statistical reliability from starting structures of up to 7 Å root-mean-square deviation in the glycoligand ring atoms. We show that GlycanDock-refined models qualitatively replicated the known binding specificity of a bacterial carbohydrate-binding module. Finally, we present a protein-glycoligand docking pipeline for generating putative protein-glycoligand complexes when only the glycoligand sequence and unbound protein structure are known. In combination with other carbohydrate modeling tools, the GlycanDock docking refinement algorithm will accelerate research in the glycosciences.

Recent advancements in polyethyleneimine-based materials and their biomedical, biotechnology, and biomaterial applications

Precise-synthesis strategies and integration approaches of bioinspired PEI-based systems, and their biomedical, biotechnology and biomaterial applications.

The Pseudomonas aeruginosa lectin LecB binds to the exopolysaccharide Psl and stabilizes the biofilm matrix

Pseudomonas aeruginosa biofilms are composed of exopolysaccharides (EPS), exogenous DNA, and proteins that hold these communities together. P. aeruginosa produces lectins LecA and LecB, which possess affinities towards sugars found in matrix EPS and mediate adherence of P. aeruginosa to target host cells. Here, we demonstrate that LecB binds to Psl, a key matrix EPS, and this leads to increased retention of both cells and EPS in a growing biofilm. This interaction is predicted to occur between the lectin and the branched side chains present on Psl. Finally, we show that LecB coordinates Psl localization in the biofilm. This constitutes a unique function for LecB and identifies it as a matrix protein that contributes to biofilm structure through EPS interactions.