Mechanism of Binding Site Conformational Switching in the CD44–Hyaluronan Protein–Carbohydrate Binding Interaction

Atomic-Resolution Experimental Structural Biology and Molecular Dynamics Simulations of Hyaluronan and Its Complexes

This review summarizes the atomic-resolution structural biology of hyaluronan and its complexes available in the Protein Data Bank, as well as published studies of atomic-resolution explicit-solvent molecular dynamics simulations on these and other hyaluronan and hyaluronan-containing systems. Advances in accurate molecular mechanics force fields, simulation methods and software, and computer hardware have supported a recent flourish in such simulations, such that the simulation publications now outnumber the structural biology publications by an order of magnitude. In addition to supplementing the experimental structural biology with computed dynamic and thermodynamic information, the molecular dynamics studies provide a wealth of atomic-resolution information on hyaluronan-containing systems for which there is no atomic-resolution structural biology either available or possible. Examples of these summarized in this review include hyaluronan pairing with other hyaluronan molecules and glycosaminoglycans, with ions, with proteins and peptides, with lipids, and with drugs and drug-like molecules. Despite limitations imposed by present-day computing resources on system size and simulation timescale, atomic-resolution explicit-solvent molecular dynamics simulations have been able to contribute significant insight into hyaluronan's flexibility and capacity for intra- and intermolecular non-covalent interactions.

Hypoxia modulates human mast cell adhesion to hyaluronic acid

Hypoxia is an inherent factor in the inflammatory process and is important in the regulation of some immune cell functions, including the expression of mast cell pro- and anti-inflammatory mediators. Hypoxia also influences cell adhesion to the extracellular matrix (ECM). Hyaluronic acid is one of the major components of the ECM that is involved in inflammatory and tissue regeneration processes in which mast cells play a prominent role. This prompted us to investigate the effects of hypoxia on the expression of hyaluronic acid receptors in mast cells and mast cell adhesion to this ECM component. We found that human LAD2 mast cells spontaneously adhered to hyaluronic acid in a CD44-dependent manner and that reduced oxygen concentrations inhibited or even completely abolished this adhesion process. The mechanism of hypoxia downregulation of mast cell adhesion to hyaluronic acid did not involve a decrease in CD44 expression and hyaluronidase-mediated degradation of adhesion substrates but rather conformational changes in the avidity of CD44 to hyaluronic acid. Hypoxia-mediated regulation of mast cell adhesion to extracellular matrix components might be involved in the pathogenic accumulation of mast cells observed in the course of certain diseases including rheumatoid arthritis and cancer.

Comparison of Carbohydrate Force Fields in Molecular Dynamics Simulations of Protein-Carbohydrate Complexes

In this paper, we present the results of molecular dynamics simulations aimed at critical comparison of classical, biomolecular force fields (FFs) in the context of their capabilities to describe the structural and thermodynamic features of carbohydrate-protein interactions. We have considered the three main families of FFs (CHARMM, GROMOS, and GLYCAM/AMBER) by applying them to investigate the seven different carbohydrate-protein complexes. The results indicate that although the qualitative pattern of several structural descriptors (intermolecular hydrogen bonding, ligand dynamic location, etc.) is conserved among the compared FFs, there also exists a number of significant divergences (mainly the patterns of contacts between particular amino acid residues and bound carbohydrate). The carbohydrate-protein unbinding free energies also vary from one FF to another, displaying diversified trends in deviations from the experimental data. The magnitude of those deviations is not negligible and indicates the need for refinement in the currently existing combinations of carbohydrate- and protein-dedicated biomolecular force fields. In spite of the lack of explicit functional terms responsible for the corresponding intermolecular forces, all tested FFs are capable of adequately reproducing the CH-π interactions, crucial for carbohydrate-protein binding.

Oxidative damage to hyaluronan-CD44 interactions as an underlying mechanism of action of oxidative stress-inducing cancer therapy

Multiple cancer therapies nowadays rely on oxidative stress to damage cancer cells. Here we investigated the biological and molecular effect of oxidative stress on the interaction between CD44 and hyaluronan (HA), as interrupting their binding can hinder cancer progression. Our experiments demonstrated that the oxidation of HA decreased its recognition by CD44, which was further enhanced when both CD44 and HA were oxidized. The reduction of CD44-HA binding negatively affected the proliferative state of cancer cells. Our multi-level atomistic simulations revealed that the binding free energy of HA to CD44 decreased upon oxidation. The effect of HA and CD44 oxidation on CD44-HA binding was similar, but when both HA and CD44 were oxidized, the effect was much larger, in agreement with our experiments. Hence, our experiments and computations support our hypothesis on the role of oxidation in the disturbance of CD44-HA interaction, which can lead to the inhibition of proliferative signaling pathways inside the tumor cell to induce cell death.

N-Glycosylation can selectively block or foster different receptor-ligand binding modes

While DNA encodes protein structure, glycans provide a complementary layer of information to protein function. As a prime example of the significance of glycans, the ability of the cell surface receptor CD44 to bind its ligand, hyaluronan, is modulated by N-glycosylation. However, the details of this modulation remain unclear. Based on atomistic simulations and NMR, we provide evidence that CD44 has multiple distinct binding sites for hyaluronan, and that N-glycosylation modulates their respective roles. We find that non-glycosylated CD44 favors the canonical sub-micromolar binding site, while glycosylated CD44 binds hyaluronan with an entirely different micromolar binding site. Our findings show (for the first time) how glycosylation can alter receptor affinity by shielding specific regions of the host protein, thereby promoting weaker binding modes. The mechanism revealed in this work emphasizes the importance of glycosylation in protein function and poses a challenge for protein structure determination where glycosylation is usually neglected.

Synthesis of hyaluronic acid oligosaccharides with a GlcNAc-GlcA repeating pattern and their binding affinity with CD44

Hyaluronic acid (HA) is a ubiquitous glycosaminoglycan in the extracellular matrix and a ligand of CD44, a transmembrane glycoprotein that is important in cell migration. Crystal and NMR studies found a hexasaccharide of the pattern (GlcA-GlcNAc)3 as the shortest HA that could bind to CD44, but molecular dynamics simulations indicated that a tetrasaccharide of the pattern (GlcNAc-GlcA)2 is the key structure interacting with CD44. Access to oligomers with such a repeat pattern is crucial in binding studies with CD44. Here we developed a synthetic procedure to afford the HA oligosaccharides with the GlcNAc-GlcA repeating unit and measured the binding interaction between these sugars and human CD44 by isothermal titration calorimetry (ITC). During the chemical synthesis, we successfully generated the β-glycosidic bond in the absence of neighbouring group participation and overcome the issues in the oxidation step. In addition, ammonia-free dissolving metal reduction for debenzylation and azido reduction has been applied in carbohydrate synthesis for the first time. ITC analysis revealed that the HA tetrasaccharide (GlcNAc-GlcA)2 could indeed interact and bind to the human CD44.

Structural Characterization of the CD44 Stem Region for Standard and Cancer-Associated Isoforms

CD44 is widely expressed in most vertebrate cells, whereas the expression of CD44v6 is restricted to only a few tissues and has been considered to be associated with tumor progression and metastasis. Thus, CD44v6 has been recognized as a promising prognostic biomarker and therapeutic target for various cancers for more than a decade. However, despite many experimental studies, the structural dynamics and differences between CD44s and CD44v6, particularly in their stem region, still remain elusive. Here, a computational study was conducted to address these problems. We found that the stem of CD44s adopted predominantly two conformations, one featuring antiparallel β-sheets and the other featuring parallel β-sheets, whereas the stem of CD44v6 adopted mainly one conformation with relatively highly suppressed β-sheet contents. Moreover, Phe215 was found to be essential in the β-sheets of both CD44s and CD44v6. We finally found intramolecular Phe215-Trp224 hydrogen-bonding interactions and hydrophobic interactions with Phe215 that cooperatively drove conformational differences upon the addition of the v6 region to CD44. Our study elucidated the structural differences between the stem regions of CD44s and CD44v6 and thus can offer useful structural information for drug design to specifically target CD44v6 in promising clinical applications.

Molecular weight impact on the mechanical forces between hyaluronan and its receptor

Hyaluronan (HA) possesses manifold mechanical and signaling properties in the body. Most of these activities are largely regulated by its molecular weight, which often triggers opposing functions. However the molecular basis for such function distinction between HA size categories remains unclear. Using a combination of biophysical techniques, we measured the physical forces between HA ligand and its specific receptor CD44 in both normal and lateral directions, at different HA molecular weights and bound states. It was found that the impact of HA multivalency is more than just the sum of separate monovalent bindings. The HA-CD44 specific interaction enhances with HA molecular weight and the maximum binding occurs at ∼1000 kD, possibly due to the balance between multivalent HA zipping effect and conformational entropy. High friction patches, probably from CD44 protein clustering, was observed in friction force microscopy (FFM) upon HA shearing, which is also dependent on HA molecular weight. These results could help to understand the biophysical mechanism of HA in regulating CD44-induced physiological activities and thus facilitate the new design of HA-based material in fine tuning the receptor responses.

What is special about 200 kDa hyaluronan that activates hyaluronan receptor signaling?

The polydispersity of hyaluronan (HA) presents challenges for analyzing its solution properties, such as the relationship between mass and particle size. The broad mass range of natural HA (≤50-fold) makes molecular characterization difficult and ambiguous compared to molecules with known molecular weights (e.g., proteins). Biophysical studies show that large >MDa HA behaves like a random coil, whereas very small (e.g., 10 kDa) HA behaves like a rod. However, the mass range for this conformational transition is not easily determined in natural polydisperse HA. Some HA receptors (e.g., CD44 and HARE) initiate signaling responses upon binding HA in the 100-300 kDa range, but not larger MDa HA. Size-dependent responses are studied using nonnatural HA: purified narrow-size range HA [Pandey MS, Baggenstoss BA, Washburn J, Harris EN, Weigel PH. 2013. The hyaluronan receptor for endocytosis (HARE) activates NF-κB-mediated gene expression in response to 40-400 kDa, but not smaller or sarger, hyaluronans. J Biol Chem. 288:14068-14079] and very narrow size range Select-HA made chemo-enzymatically [Jing W, DeAngelis PL. 2004. Synchronized chemoenzymatic synthesis of monodisperse hyaluronan polymers. J Biol Chem. 279:42345-42349]. Here, we used size exclusion chromatography and multiangle light scattering to determine the weight-average molar mass and diameter of ~60 very narrow size preparations from 29 to 1650 kDa. The ratio of HA mass to HA diameter showed a transition in the 150-250 kDa size range (~65 nm). The HA rod-to-coil transition occurs within the size range that specifically activates cell signaling by some receptors. Thus, size-specific signaling could be due to unique external receptor•HA conformation changes that enable transmembrane-mediated activation of cytoplasmic domains. Alternatively and more likely, transition-size HA may enable multiple receptors to bind the same HA, creating new internal signal-competent cytoplasmic domain complexes.

Atomistic fingerprint of hyaluronan-CD44 binding

Hyaluronan is a polyanionic, megadalton-scale polysaccharide, which initiates cell signaling by interacting with several receptor proteins including CD44 involved in cell-cell interactions and cell adhesion. Previous studies of the CD44 hyaluronan binding domain have identified multiple widespread residues to be responsible for its recognition capacity. In contrast, the X-ray structural characterization of CD44 has revealed a single binding mode associated with interactions that involve just a fraction of these residues. In this study, we show through atomistic molecular dynamics simulations that hyaluronan can bind CD44 with three topographically different binding modes that in unison define an interaction fingerprint, thus providing a plausible explanation for the disagreement between the earlier studies. Our results confirm that the known crystallographic mode is the strongest of the three binding modes. The other two modes represent metastable configurations that are readily available in the initial stages of the binding, and they are also the most frequently observed modes in our unbiased simulations. We further discuss how CD44, fostered by the weaker binding modes, diffuses along HA when attached. This 1D diffusion combined with the constrained relative orientation of the diffusing proteins is likely to influence the aggregation kinetics of CD44. Importantly, CD44 aggregation has been suggested to be a possible mechanism in CD44-mediated signaling.

Hyaluronan is a natural sugar polymer in our bodies. Besides acting as a space-filling agent for example in multiple connective tissues, it also functions as a cellular cue in cancer and inflammation. Our tissues sense hyaluronan through receptors—proteins that sit at the surface of cells and grab the molecules they are expected to recognize. Although the knowledge associated with hyaluronan and its receptors is constantly accumulating, the molecular-level insight is largely missing or incomplete due to the lack of techniques able to probe the dynamics of protein–carbohydrate interactions with sufficiently high resolution. In this work, we characterize the binding of hyaluronan to its receptor CD44 with atomistic precision. We achieve this level of precision by employing atomistic molecular dynamics simulations. This computational technique allows one to follow the movement of atoms of a virtual system at scales beyond the resolution of any experimental technique. Our work specifically focuses on the different stages of hyaluronan–CD44 binding, and we observe the process to involve three different binding modes, making it more versatile than previously thought. Our insights, therefore, promote the understanding of the interplay between hyaluronan and HA, thereby fostering development of new drugs or inhibitors to malignancies, such as cancer metastasis.

Rigidity and flexibility in the tetrasaccharide linker of proteoglycans from atomic-resolution molecular simulation

Proteoglycans (PGs) are covalent conjugates between protein and carbohydrate (glycosaminoglycans). Certain classes of glycosaminoglycans such as chondroitin sulfate/dermatan sulfate and heparan sulfate utilize a specific tetrasaccharide linker for attachment to the protein component: GlcAβ1-3Galβ1-3Galβ1-4Xylβ1-O-Ser. Toward understanding the conformational preferences of this linker, the present work used all-atom explicit-solvent molecular dynamics (MD) simulations combined with Adaptive Biasing Force (ABF) sampling to determine high-resolution, high-precision conformational free energy maps ΔG(φ, ψ) for each glycosidic linkage between constituent disaccharides, including the variant where GlcA is substituted with IdoA. These linkages are characterized by single, predominant (> 97% occupancy), and broad (45° × 60° for ΔG(φ, ψ) < 1 kcal/mol) free-energy minima, while the Xyl-Ser linkage has two such minima similar in free-energy, and additional flexibility from the Ser sidechain dihedral. Conformational analysis of microsecond-scale standard MD on the complete tetrasaccharide-O-Ser conjugate is consistent with ABF data, suggesting (φ, ψ) probabilities are independent of the linker context, and that the tetrasaccharide acts as a relatively rigid unit whereas significant conformational heterogeneity exists with respect to rotation about bonds connecting Xyl to Ser. © 2017 Wiley Periodicals, Inc.

Ligand binding to anti-cancer target CD44 investigated by molecular simulations

CD44 is a cell-surface glycoprotein and receptor for hyaluronan, one of the major components of the tumor extracellular matrix. There is evidence that the interaction between CD44 and hyaluronan promotes breast cancer metastasis. Recently, the molecule F-19848A was shown to inhibit hyaluronan binding to receptor CD44 in a cell-based assay. In this study, we investigated the mechanism and energetics of F-19848A binding to CD44 using molecular simulation. Using the molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) method, we obtained the binding free energy and inhibition constant of the complex. The van der Waals (vdW) interaction and the extended portion of F-19848A play key roles in the binding affinity. We screened natural products from a traditional Chinese medicine database to search for CD44 inhibitors. From combining pharmaceutical requirements with docking and molecular dynamics simulations, we found ten compounds that are potentially better or equal to the F-19848A ligand at binding to CD44 receptor. Therefore, we have identified new candidates of CD44 inhibitors, based on molecular simulation, which may be effective small molecules for the therapy of breast cancer.

Native proteins trap high-energy transit conformations

During protein folding and as part of some conformational changes that regulate protein function, the polypeptide chain must traverse high-energy barriers that separate the commonly adopted low-energy conformations. How distortions in peptide geometry allow these barrier-crossing transitions is a fundamental open question. One such important transition involves the movement of a non-glycine residue between the left side of the Ramachandran plot (that is, ϕ < 0°) and the right side (that is, ϕ > 0°). We report that high-energy conformations with ϕ ~ 0°, normally expected to occur only as fleeting transition states, are stably trapped in certain highly resolved native protein structures and that an analysis of these residues provides a detailed, experimentally derived map of the bond angle distortions taking place along the transition path. This unanticipated information lays to rest any uncertainty about whether such transitions are possible and how they occur, and in doing so lays a firm foundation for theoretical studies to better understand the transitions between basins that have been little studied but are integrally involved in protein folding and function. Also, the context of one such residue shows that even a designed highly stable protein can harbor substantial unfavorable interactions.

Revealing the Mechanisms of Protein Disorder and N-Glycosylation in CD44-Hyaluronan Binding Using Molecular Simulation

The extracellular N-terminal hyaluronan binding domain (HABD) of CD44 is a small globular domain that confers hyaluronan (HA) binding functionality to this large transmembrane glycoprotein. When recombinantly expressed by itself, HABD exists as a globular water-soluble protein that retains the capacity to bind HA. This has enabled atomic-resolution structural biology experiments that have revealed the structure of HABD and its binding mode with oligomeric HA. Such experiments have also pointed to an order-to-disorder transition in HABD that is associated with HA binding. However, it had remained unclear how this structural transition was involved in binding since it occurs in a region of HABD distant from the HA-binding site. Furthermore, HABD is known to be N-glycosylated, and such glycosylation can diminish HA binding when the associated N-glycans are capped with sialic acid residues. The intrinsic flexibility of disordered proteins and of N-glycans makes it difficult to apply experimental structural biology approaches to probe the molecular mechanisms of how the order-to-disorder transition and N-glycosylation can modulate HA binding by HABD. We review recent results from molecular dynamics simulations that provide atomic-resolution mechanistic understanding of such modulation to help bridge gaps between existing experimental binding and structural biology data. Findings from these simulations include: Tyr42 may function as a molecular switch that converts the HA-binding site from a low affinity to a high affinity state; in the partially disordered form of HABD, basic amino acids in the C-terminal region can gain sufficient mobility to form direct contacts with bound HA to further stabilize binding; and terminal sialic acids on covalently attached N-glycans can form charge-paired hydrogen bonding interactions with basic amino acids that could otherwise bind to HA, thereby blocking HA binding to glycosylated CD44 HABD.

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.

Terminal sialic acids on CD44 N-glycans can block hyaluronan binding by forming competing intramolecular contacts with arginine sidechains

Specific sugar residues and their linkages form the basis of molecular recognition for interactions of glycoproteins with other biomolecules. Seemingly small changes, like the addition of a single monosaccharide in the covalently attached glycan component of glycoproteins, can greatly affect these interactions. For instance, the sialic acid capping of glycans affects protein-ligand binding involved in cell-cell and cell-matrix interactions. CD44 is a single-pass transmembrane glycoprotein whose binding with its carbohydrate ligand hyaluronan (HA), an extracellular matrix component, mediates processes such as leukocyte homing, cell adhesion, and tumor metastasis. This binding is highly regulated by glycosylation of the N-terminal extracellular hyaluronan-binding domain (HABD); specifically, sialic acid capped N-glycans of HABD inhibit ligand binding. However, the molecular mechanism behind this sialic acid mediated regulation has remained unknown. Two of the five N-glycosyation sites of HABD have been previously identified as having the greatest inhibitory effect on HA binding, but only if the glycans contain terminal sialic acid residues. These two sites, Asn25 and Asn120, were chosen for in silico glycosylation in this study. Here, from extensive standard molecular dynamics simulations and biased simulations, we propose a molecular mechanism for this behavior based on spontaneously-formed charge-paired hydrogen bonding interactions between the negatively-charged sialic acid residues and positively-charged Arg sidechains known to be critically important for binding to HA, which itself is negatively charged. Such intramolecular hydrogen bonds would preclude associations critical to hyaluronan binding. This observation suggests how CD44 and related glycoprotein binding is regulated by sialylation as cellular environments fluctuate.

CD44 receptor unfolding enhances binding by freeing basic amino acids to contact carbohydrate ligand

The extracellular carbohydrate-binding domain of the Type I transmembrane receptor CD44 is known to undergo affinity switching, where change in conformation leads to enhanced binding of its carbohydrate ligand hyaluronan. Separate x-ray crystallographic and NMR experiments have led to competing explanations, with the former supporting minor conformational changes at the binding site and the latter a major order-to-disorder unfolding transition distant from the binding site. Here, all-atom explicit-solvent molecular dynamics studies employing adaptive biasing force sampling revealed a substantial favorable free-energy change associated with contact formation between the Arg(41) side chain and hyaluronan at the binding site, independent of whether the distant site was ordered or disordered. Analogous computational experiments on Arg(41)Ala mutants showed loss of this favorable free-energy change, consistent with existing experimental data. More provocatively, the simulation data revealed the molecular mechanism by which the order-to-disorder transition enhances hyaluronan binding: in the disordered state, a number of basic residues gain sufficient conformational freedom-lacking in the ordered state-to spontaneously form side-chain contacts with hyaluronan. Mutation of these residues to Ala had been known to decrease binding affinity, but there had previously been no structural explanation, given their lack of proximity to the carbohydrate-binding site in existing structures of the complex.

Peptide backbone sampling convergence with the adaptive biasing force algorithm

Complete Boltzmann sampling of reaction coordinates in biomolecular systems continues to be a challenge for unbiased molecular dynamics simulations. A growing number of methods have been developed for applying biases to biomolecular systems to enhance sampling while enabling recovery of the unbiased (Boltzmann) distribution of states. The adaptive biasing force (ABF) algorithm is one such method and works by canceling out the average force along the desired reaction coordinate(s) using an estimate of this force progressively accumulated during the simulation. Upon completion of the simulation, the potential of mean force, and therefore Boltzmann distribution of states, is obtained by integrating this average force. In an effort to characterize the expected performance in applications such as protein loop sampling, ABF was applied to the full ranges of the Ramachandran φ/ψ backbone dihedral reaction coordinates for dipeptides of the 20 amino acids using all-atom explicit-water molecular dynamics simulations. Approximately half of the dipeptides exhibited robust and rapid convergence of the potential of mean force as a function of φ/ψ in triplicate 50 ns simulations, while the remainder exhibited varying degrees of less complete convergence. The greatest difficulties in achieving converged ABF sampling were seen in the branched-side chain amino acids threonine and valine, as well as the special case of proline. Proline dipeptide sampling was further complicated by trans-to-cis peptide bond isomerization not observed in unbiased control molecular dynamics simulations. Overall, the ABF method was found to be a robust means of sampling the entire φ/ψ reaction coordinate for the 20 amino acids, including high free-energy regions typically inaccessible in standard molecular dynamics simulations.

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.

CHARMM Additive All-Atom Force Field for Phosphate and Sulfate Linked to Carbohydrates

Presented is an extension of the CHARMM additive all-atom carbohydrate force field to enable the modeling of phosphate and sulfate linked to carbohydrates. The parameters are developed in a hierarchical fashion using model compounds containing the key atoms in the full carbohydrates. Target data for parameter optimization included full two-dimensional energy surfaces defined by the glycosidic dihedral angle pairs in the phosphate/sulfate model compound analogs of hexopyranose monosaccharide phosphates and sulfates, as determined by quantum mechanical (QM) MP2/cc-pVTZ single point energies on MP2/6-31+G(d) optimized structures. In order to achieve balanced, transferable dihedral parameters for the dihedral angles, surfaces for all possible anomeric and conformational states were included during the parametrization process. In addition, to model physiologically relevant systems both the mono- and di-anionic charged states were studied for the phosphates. This resulted in over 7000 MP2/cc-pVTZ//MP2/6-31G+(d) model compound conformational energies which, supplemented with QM geometries, were the main target data for the parametrization. Parameters were validated against crystals of relevant monosaccharide derivatives obtained from the Cambridge Structural Database (CSD) and larger systems, namely inositol-(tri/tetra/penta) phosphates non-covalently bound to the pleckstrin homology (PH) domain and oligomeric chondroitin sulfate in solution and in complex with cathepsin K protein.

CHARMM additive all-atom force field for carbohydrate derivatives and its utility in polysaccharide and carbohydrate-protein modeling

Monosaccharide derivatives such as xylose, fucose, N-acetylglucosamine (GlcNAc), N-acetylgalactosamine (GlaNAc), glucuronic acid, iduronic acid, and N-acetylneuraminic acid (Neu5Ac) are important components of eukaryotic glycans. The present work details development of force-field parameters for these monosaccharides and their covalent connections to proteins via O-linkages to serine or threonine sidechains and via N-linkages to asparagine sidechains. The force field development protocol was designed to explicitly yield parameters that are compatible with the existing CHARMM additive force field for proteins, nucleic acids, lipids, carbohydrates, and small molecules. Therefore, when combined with previously developed parameters for pyranose and furanose monosaccharides, for glycosidic linkages between monosaccharides, and for proteins, the present set of parameters enables the molecular simulation of a wide variety of biologically-important molecules such as complex carbohydrates and glycoproteins. Parametrization included fitting to quantum mechanical (QM) geometries and conformational energies of model compounds, as well as to QM pair interaction energies and distances of model compounds with water. Parameters were validated in the context of crystals of relevant monosaccharides, as well NMR and/or x-ray crystallographic data on larger systems including oligomeric hyaluronan, sialyl Lewis X, O- and N-linked glycopeptides, and a lectin:sucrose complex. As the validated parameters are an extension of the CHARMM all-atom additive biomolecular force field, they further broaden the types of heterogeneous systems accessible with a consistently-developed force-field model.