Conformational determinants of the activity of antiproliferative factor glycopeptide

The effect of a methyl group on structure and function: Serine vs. threonine glycosylation and phosphorylation

A variety of glycan structures cover the surface of all cells and are involved in myriad biological processes, including but not limited to, cell adhesion and communication, protein quality control, signal transduction and metabolism, while also being intimately involved in innate and adaptive immune functions. Immune surveillance and responses to foreign carbohydrate antigens, such as capsular polysaccharides on bacteria and surface protein glycosylation of viruses, are the basis of microbial clearance, and most antimicrobial vaccines target these structures. In addition, aberrant glycans on tumors called Tumor-Associated Carbohydrate Antigens (TACAs) elicit immune responses to cancer, and TACAs have been used in the design of many antitumor vaccine constructs. A majority of mammalian TACAs are derived from what are referred to as mucin-type O-linked glycans on cell-surface proteins and are linked to the protein backbone through the hydroxyl group of either serine or threonine residues. A small group of structural studies that have compared mono- and oligosaccharides attached to each of these residues have shown that there are distinct differences in conformational preferences assumed by glycans attached to either "unmethylated" serine or ß-methylated threonine. This suggests that the linkage point of antigenic glycans will affect their presentation to the immune system as well as to various carbohydrate binding molecules (e.g., lectins). This short review, followed by our hypothesis, will examine this possibility and extend the concept to the presentation of glycans on surfaces and in assay systems where recognition of glycans by proteins and other binding partners can be defined by different attachment points that allow for a range of conformational presentations.

Glycoamino Acid Analogues of the Thomsen-Friedenreich Tumor-Associated Carbohydrate Antigen: Synthesis and Evaluation of Novel Antiproliferative Factor Glycopeptides

Glycoamino acid analogues of the Thomsen-Friedenreich antigen disaccharide, where the 4' and 4″ hydroxyl groups were substituted with fluorine or hydrogen, were synthesized and incorporated into the asialylated antiproliferative factor (as-APF), a biologically active form of APF, a glycopeptide found in the urine of patients with interstitial cystitis. Various strategies were employed to incorporate the fluorine atom at the 4-positions of either the galactose or N-acetylgalactosamine unit of the disaccharide antigen, based on stereochemistry and reactivity. These glycopeptides were evaluated in antiproliferative assays on both primary normal bladder epithelial cells and T24 bladder carcinoma cells. Unlike many previously published substitutions to APF, mono-4'-fluorination of the GalNAc residue did not affect the activity, whereas fluoro-derivatives of the galactose 4″-position or both 4' and 4″ hydroxyls showed a reduced potency relative to the monosubstituted GalNAc derivative. A fourth compound where the 4″ position of galactose was deoxygenated showed a lower potency than the parent and monosubstituted compounds. These results suggest that specific substitutions in the sugar moieties in the APF can be tolerated, and the glycomimetic design of APF analogues can include fluorine in the GalNAc sugar of the disaccharide.

Antiproliferative factor (APF) binds specifically to sites within the cytoskeleton-associated protein 4 (CKAP4) extracellular domain

Background

Antiproliferative factor (APF) is a sialoglycopeptide elevated in the urine of patients with interstitial cystitis—a chronic, painful bladder disease. APF inhibits the proliferation of normal bladder epithelial cells and cancer cells in vitro, presumably by binding to its cellular receptor, cytoskeleton associated-protein 4 (CKAP4); however, the biophysical interaction of APF with CKAP4 has not been characterized previously. In this study, we used surface plasmon resonance (SPR) to explore the binding kinetics of the interaction of APF and as-APF (a desialylated APF analogue with full activity) to CKAP4.

Results

We immobilized non-glycosylated APF (TVPAAVVVA) to the Fc1 channel as the control and as-APF to Fc2 channel as the ligand in order to measure the binding of CKAP4 recombinant proteins encompassing only the extracellular domain (Aa 127–602) or the extracellular domain plus the transmembrane domain (Aa 106–602). Positive binding was detected to both CKAP4_126–602 and CKAP4_106–602, suggesting that as-APF can bind specifically to CKAP4 and that the potential binding site(s) are located within the extracellular domain. To identify the primary APF binding site(s) within the CKAP4 extracellular domain, deletion mutants were designed according to structural predictions, and the purified recombinant proteins were immobilized on a CM5 chip through amine-coupling to measure as-APF binding activity. Importantly, both CKAP4_127–360 and CKAP4_361–524 exhibited a fast association rate (k_on) and a slow dissociation rate (k_off), thus generating high binding affinity and suggesting that both regions contribute relatively equally to overall as-APF binding. Therefore, two or more as-APF binding sites may exist within the Aa 127–524 region of the CKAP4 extracellular domain.

Conclusions

We determined that the CKAP4_127–360 and CKAP4_361–524 mutants exhibit improved binding activity to as-APF as compared to the full-length extracellular domain, making it possible to detect low concentrations of as-APF in urine, thereby establishing a foundation for a non-invasive diagnostic assay for IC. Further, these data have revealed novel APF binding site(s) suggesting that targeting this region of CKAP4 to inhibit APF binding may be a useful strategy for treating IC-related bladder pathology.

Delineating the conformational flexibility of trisaccharides from NMR spectroscopy experiments and computer simulations

The conformation of saccharides in solution is challenging to characterize in the context of a single well-defined three-dimensional structure. Instead, they are better represented by an ensemble of conformations associated with their structural diversity and flexibility. In this study, we delineate the conformational heterogeneity of five trisaccharides via a combination of experimental and computational techniques. Experimental NMR measurements target conformationally sensitive parameters, including J couplings and effective distances around the glycosidic linkages, while the computational simulations apply the well-calibrated additive CHARMM carbohydrate force field in combination with efficient enhanced sampling molecular dynamics simulation methods. Analysis of conformational heterogeneity is performed based on sampling of discreet states as defined by dihedral angles, on root-mean-square differences of Cartesian coordinates and on the extent of volume sampled. Conformational clustering, based on the glycosidic linkage dihedral angles, shows that accounting for the full range of sampled conformations is required to reproduce the experimental data, emphasizing the utility of the molecular simulations in obtaining an atomic detailed description of the conformational properties of the saccharides. Results show the presence of differential conformational preferences as a function of primary sequence and glycosidic linkage types. Significant differences in conformational ensembles associated with the anomeric configuration of a single glycosidic linkage reinforce the impact of such changes on the conformational properties of carbohydrates. The present structural insights of the studied trisaccharides represent a foundation for understanding the range of conformations adopted in larger oligosaccharides and how these molecules encode their conformational heterogeneity into the monosaccharide sequence.

Enhanced conformational sampling using replica exchange with concurrent solute scaling and hamiltonian biasing realized in one dimension

Replica exchange (REX) is a powerful computational tool for overcoming the quasi-ergodic sampling problem of complex molecular systems. Recently, several multidimensional extensions of this method have been developed to realize exchanges in both temperature and biasing potential space or the use of multiple biasing potentials to improve sampling efficiency. However, increased computational cost due to the multidimensionality of exchanges becomes challenging for use on complex systems under explicit solvent conditions. In this study, we develop a one-dimensional (1D) REX algorithm to concurrently combine the advantages of overall enhanced sampling from Hamiltonian solute scaling and the specific enhancement of collective variables using Hamiltonian biasing potentials. In the present Hamiltonian replica exchange method, termed HREST-BP, Hamiltonian solute scaling is applied to the solute subsystem, and its interactions with the environment to enhance overall conformational transitions and biasing potentials are added along selected collective variables associated with specific conformational transitions, thereby balancing the sampling of different hierarchical degrees of freedom. The two enhanced sampling approaches are implemented concurrently allowing for the use of a small number of replicas (e.g., 6 to 8) in 1D, thus greatly reducing the computational cost in complex system simulations. The present method is applied to conformational sampling of two nitrogen-linked glycans (N-glycans) found on the HIV gp120 envelope protein. Considering the general importance of the conformational sampling problem, HREST-BP represents an efficient procedure for the study of complex saccharides, and, more generally, the method is anticipated to be of general utility for the conformational sampling in a wide range of macromolecular systems.

Conformational sampling of oligosaccharides using Hamiltonian replica exchange with two-dimensional dihedral biasing potentials and the weighted histogram analysis method (WHAM)

Oligosaccharides and polysaccharides exert numerous functional roles in biology through their structural diversity and conformational properties. To investigate their conformational properties using computational methods, Hamiltonian replica exchange (H-REX) combined with two-dimensional grid-based correction maps as biasing potentials (bpCMAP) significantly improves the sampling efficiency about glycosidic linkages. In the current study, we extend the application of H-REX with bpCMAP to complex saccharides and establish systematic procedures for bpCMAP construction, determination of replica distribution, and data analysis. Our main findings are that (1) the bpCMAP for each type of glycosidic linkage can be constructed from the corresponding disaccharide using gas-phase umbrella sampling simulations, (2) the replica distribution can be conveniently determined following the exact definition of the average acceptance ratio based on the assigned distribution of biasing potentials, and (3) the extracted free energy surface (or potential of mean force (PMF)) can be improved using the weighted histogram analysis method (WHAM) allowing for the inclusion of data from the excited state replicas in the calculated probability distribution. The method is applied to a branched N-glycan found on the HIV gp120 protein, and a linear N-glycan. Considering the general importance of N-glycans and the wide appreciation of the sampling problem, the present method represents an efficient procedure for the conformational sampling of complex oligo- and polysaccharides under explicit solvent conditions. More generally, the use of WHAM is anticipated to be of general utility for the calculation of PMFs from H-REX simulations in a wide range of macromolecular systems.

CHARMM additive and polarizable force fields for biophysics and computer-aided drug design

BACKGROUND

Molecular Mechanics (MM) is the method of choice for computational studies of biomolecular systems owing to its modest computational cost, which makes it possible to routinely perform molecular dynamics (MD) simulations on chemical systems of biophysical and biomedical relevance.

SCOPE OF REVIEW

As one of the main factors limiting the accuracy of MD results is the empirical force field used, the present paper offers a review of recent developments in the CHARMM additive force field, one of the most popular biomolecular force fields. Additionally, we present a detailed discussion of the CHARMM Drude polarizable force field, anticipating a growth in the importance and utilization of polarizable force fields in the near future. Throughout the discussion emphasis is placed on the force fields' parametrization philosophy and methodology.

MAJOR CONCLUSIONS

Recent improvements in the CHARMM additive force field are mostly related to newly found weaknesses in the previous generation of additive force fields. Beyond the additive approximation is the newly available CHARMM Drude polarizable force field, which allows for MD simulations of up to 1μs on proteins, DNA, lipids and carbohydrates.

GENERAL SIGNIFICANCE

Addressing the limitations ensures the reliability of the new CHARMM36 additive force field for the types of calculations that are presently coming into routine computational reach while the availability of the Drude polarizable force fields offers an inherently more accurate model of the underlying physical forces driving macromolecular structures and dynamics. This article is part of a Special Issue entitled "Recent developments of molecular dynamics".

Synthesis and cell-selective antitumor properties of amino acid conjugated tumor-associated carbohydrate antigen-coated gold nanoparticles

The Thomsen Friedenreich antigen (TFag) disaccharide is a tumor-associated carbohydrate antigen (TACA) found primarily on carcinoma cells and rarely expressed in normal tissue. The TFag has been shown to interact with Galectin-3 (Gal-3), one in a family of β-galactoside binding proteins. Galectins have a variety of cellular functions, and Gal-3 has been shown to be the sole galectin with anti-apoptotic activity. We have previously prepared gold nanoparticles (AuNP) coated with the TFag in various presentations as potential anti-adhesive therapeutic tools or antitumor vaccine platforms. Here we describe the synthesis of TFag-glycoamino acid conjugates attached to gold nanoparticles through a combined alkane/PEG linker, where the TFag was attached to either a serine or threonine amino acid. Particles were fully characterized by a host of biophysical techniques, and along with a control particle carrying hydroxyl-terminated linker units, were evaluated in both Gal-3 positive and negative cell lines. We show that the particles bearing the saccharides selectively inhibited tumor cell growth of the Gal-3 positive cells significantly more than the Gal-3 negative cells. In addition, the threonine-attached TF particles were more potent than the serine-attached constructs. These results support the use of AuNP as antitumor therapeutic platforms, targeted against cell lines that express specific lectins that interact with TFag.

Perturbation of long-range water dynamics as the mechanism for the antifreeze activity of antifreeze glycoprotein

Very little is known about the mechanism of antifreeze action of antifreeze glycoproteins (AFGPs) present in Antarctic teleost fish. Recent NMR and CD studies assisted with total synthesis of synthetic AFGP variants have provided insight into the structure of short AFGP glycopeptides, though the observations did not yield information on the antifreeze mechanism of action. In this study, we use Hamiltonian replica exchange (HREX) molecular dynamics simulations to probe the structure and surrounding aqueous environments of both the natural (AFGP8) and synthetic (s-AFGP4) AFGPs. AFGPs can adopt both amphiphilic and pseudoamphiphilic conformations, the preference of which is related to the proline content of the peptide. The arrangement of carbohydrates allows the hydroxyl groups on terminal galactose units to form stable water bridges which in turn influence the hydrogen-bond network, structure, and dynamics of the surrounding solvent. Interestingly, these local effects lead to the perturbation of the tetrahedral environment for water molecules in hydration layers far (10.0-12.0 Å) from the AFGPs. This structure-induced alteration of long-range hydration dynamics is proposed to be the major contributor to antifreeze activity, a conclusion that is in line with terahertz spectroscopy experiments. The detailed structure-mechanism correlation provided in this study could lead to the design of better synthetic AFGP variants.

Conformational properties of α- or β-(1→6)-linked oligosaccharides: Hamiltonian replica exchange MD simulations and NMR experiments

Conformational sampling for a set of 10 α- or β-(1→6)-linked oligosaccharides has been studied using explicit solvent Hamiltonian replica exchange (HREX) simulations and NMR spectroscopy techniques. Validation of the force field and simulation methodology is done by comparing calculated transglycosidic J coupling constants and proton-proton distances with the corresponding NMR data. Initial calculations showed poor agreement, for example, with >3 Hz deviation of the calculated (3)J(H5,H6R) values from the experimental data, prompting optimization of the ω torsion angle parameters associated with (1→6)-linkages. The resulting force field is in overall good agreement (i.e., within ∼0.5 Hz deviation) from experimental (3)J(H5,H6R) values, although some small limitations are evident. Detailed hydrogen bonding analysis indicates that most of the compounds lack direct intramolecular H-bonds between the two monosaccharides; however, minor sampling of the O6···HO2' hydrogen bond is present in three compounds. The results verify the role of the gauche effect between O5 and O6 atoms in gluco- and manno-configured pyranosides causing the ω torsion angle to sample an equilibrium between the gt and gg rotamers. Conversely, galacto-configured pyranosides sample a population distribution in equilibrium between gt and tg rotamers, while the gg rotamer populations are minor. Water radial distribution functions suggest decreased accessibility to the O6 atom in the (1→6)-linkage as compared to the O6' atom in the nonreducing sugar. The role of bridging water molecules between two sugar moieties on the distributions of ω torsion angles in oligosaccharides is also explored.