Thermodynamics of oligosaccharide binding to a monoclonal antibody specific for a Salmonella O-antigen point to hydrophobic interactions in the binding site

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Ground- and Excited-State Interactions of a Psoralen Derivative with Human Telomeric G-Quadruplex DNA

G-quadruplex DNA has been a recent target for anticancer agents, and its binding interactions with small molecules, often used as anticancer drugs, have become an important area of research. Considering that psoralens have long been studied in the context of duplex DNA but that very little is known about their potential as G-quadruplex binders and their excited-state interaction with the latter has not been explored, we have studied herein the binding of a planar water-soluble psoralen derivative, 4'-aminomethyl-4,5',8-trimethylpsoralen (AMT), with the 22-mer human telomeric G-quadruplex-forming sequence, AGGG(TTAGGG)₃, labeled here as (hTel22), and investigated the consequences of photoexcitation of AMT by calorimetric and spectroscopic techniques. The results show an enthalpy-driven 1:1 binding of AMT with hTel22 via end-stacking mode. Fluorescence quenching experiments on 6-fluorescein amidite-labeled oligomers indicate that the binding site is nearer to the 3' end of hTel22 in the diagonal loop region. Femtosecond time-resolved transient absorption measurements indicate electron transfer from the guanine moiety of hTel22 to photoexcited AMT, leading to the formation of a radical pair species (AMT^•-G^•+), which survives for 30 ps and is favored by a parallel/quasi-parallel orientation between the two. The findings reveal psoralens as a prospective class of compounds for the development of anticancer therapeutics by targeting the G-quadruplex DNA.

Conformational Preference of Serogroup B Salmonella O Polysaccharide in Presence and Absence of the Monoclonal Antibody Se155-4

Salmonella infection is a major public health problem worldwide. Antibodies directed toward the O polysaccharide (OPS) of S. Typhimurium, a serogroup B nontyphoidal Salmonella serovar, have protected against fatal infection in animal models. The OPS is known to undergo O-acetylation, though the impact of these modifications on antibody binding is poorly understood. Using molecular simulations, we assessed the conformational properties and antigen-antibody interactions of deacetylated and O-acetylated S. Typhimurium OPS when bound by monoclonal anti-OPS IgG Se155-4. Our findings indicate that (i) the α-d-abequose (8) monosaccharide makes important interactions with Se155-4, (ii) the deacetylated form binds to the antibody in two conformations, (iii) the acetyl group at α-l-rhamnose (5) traps the acetylated O-antigenic saccharide in one of those two conformations when bound to the antibody; (iv) the dominant conformation sampled by both unbound saccharides only occurs in the deacetylated-antibody complex; and (v) both unbound saccharides sample the second bound conformation to a small extent (2-4%). These observations provide insights into the conformational preference of an antigenic saccharide when bound to a well-characterized specific monoclonal antibody, and suggest possible important properties of vaccine induced antibodies following immunization with live attenuated and OPS-based subunit vaccines.

Multivalent display of minimal Clostridium difficile glycan epitopes mimics antigenic properties of larger glycans

Synthetic cell-surface glycans are promising vaccine candidates against Clostridium difficile. The complexity of large, highly antigenic and immunogenic glycans is a synthetic challenge. Less complex antigens providing similar immune responses are desirable for vaccine development. Based on molecular-level glycan–antibody interaction analyses, we here demonstrate that the C. difficile surface polysaccharide-I (PS-I) can be resembled by multivalent display of minimal disaccharide epitopes on a synthetic scaffold that does not participate in binding. We show that antibody avidity as a measure of antigenicity increases by about five orders of magnitude when disaccharides are compared with constructs containing five disaccharides. The synthetic, pentavalent vaccine candidate containing a peptide T-cell epitope elicits weak but highly specific antibody responses to larger PS-I glycans in mice. This study highlights the potential of multivalently displaying small oligosaccharides to achieve antigenicity characteristic of larger glycans. The approach may result in more cost-efficient carbohydrate vaccines with reduced synthetic effort.

Immunologically-active glycans are promising vaccine candidates but can be difficult to synthesize. Here, the authors show that pentavalent display of a minimal disaccharde epitope on a chemical scaffold can mimic a native C. difficile glycan antigen, representing a simple approach to synthetic vaccine production.

Specific Electrostatic Molecular Recognition in Water

The identification of pairs of small peptides that recognize each other in water exclusively through electrostatic interactions is reported. The target peptide and a structure-biased combinatorial ligand library consisting of ≈78 125 compounds were synthesized on different sized beads. Peptide-peptide interactions could conveniently be observed by clustering of the small, fluorescently labeled target beads on the surface of larger ligand-containing beads. Sequences of isolated hits were determined by MS/MS. The interactions of the complex showing the highest affinity were investigated by a novel single-bead binding assay and by 2D NMR spectroscopy. Molecular dynamics (MD) studies revealed a putative mode of interaction for this unusual electrostatic binding event. High binding specificity occurred through a combination of topological matching and electrostatic and hydrogen-bond complementarities. From MD simulations binding also seemed to involve three tightly bound water molecules in the interface between the binding partners. Binding constants in the submicromolar range, useful for biomolecular adhesion and in nanostructure design, were measured.

Antibody recognition of carbohydrate epitopes†

Carbohydrate antigens are valuable as components of vaccines for bacterial infectious agents and human immunodeficiency virus (HIV), and for generating immunotherapeutics against cancer. The crystal structures of anti-carbohydrate antibodies in complex with antigen reveal the key features of antigen recognition and provide information that can guide the design of vaccines, particularly synthetic ones. This review summarizes structural features of anti-carbohydrate antibodies to over 20 antigens, based on six categories of glyco-antigen: (i) the glycan shield of HIV glycoproteins; (ii) tumor epitopes; (iii) glycolipids and blood group A antigen; (iv) internal epitopes of bacterial lipopolysaccharides; (v) terminal epitopes on polysaccharides and oligosaccharides, including a group of antibodies to Kdo-containing Chlamydia epitopes; and (vi) linear homopolysaccharides.

Specificity of furanoside-protein recognition through antibody engineering and molecular modeling

Recognition of furanosides (five-membered ring sugars) by proteins plays important roles in host-pathogen interactions. In comparison to their six-membered ring counterparts (pyranosides), detailed studies of the molecular motifs involved in the recognition of furanosides by proteins are scarce. Here the first in-depth molecular characterization of a furanoside-protein interaction system, between an antibody (CS-35) and cell wall polysaccharides of mycobacteria, including the organism responsible for tuberculosis is reported. The approach was centered on the generation of the single chain variable fragment of CS-35 and a rational library of its mutants. Investigating the interaction from various aspects revealed the structural motifs that govern the interaction, as well as the relative contribution of molecular forces involved in the recognition. The specificity of the recognition was shown to originate mainly from multiple CH-π interactions and, to a lesser degree, hydrogen bonds formed in critical distances and geometries.

Groove-type recognition of chlamydiaceae-specific lipopolysaccharide antigen by a family of antibodies possessing an unusual variable heavy chain N-linked glycan

The structure of the antigen binding fragment of mAb S25-26, determined to 1.95 Å resolution in complex with the Chlamydiaceae family-specific trisaccharide antigen Kdo(2→8)Kdo(2→4)Kdo (Kdo = 3-deoxy-α-d-manno-oct-2-ulopyranosonic acid), displays a germ-line-coded paratope that differs significantly from previously characterized Chlamydiaceae-specific mAbs despite being raised against the identical immunogen. Unlike the terminal Kdo recognition pocket that promotes cross-reactivity in S25-2-type antibodies, S25-26 and the closely related S25-23 utilize a groove composed of germ-line residues to recognize the entire trisaccharide antigen and so confer strict specificity. Interest in S25-23 was sparked by its rare high μm affinity and strict specificity for the family-specific trisaccharide antigen; however, only the related antibody S25-26 proved amenable to crystallization. The structures of three unliganded forms of S25-26 have a labile complementary-determining region H3 adjacent to significant glycosylation of the variable heavy chain on asparagine 85 in Framework Region 3. Analysis of the glycan reveals a heterogeneous mixture with a common root structure that contains an unusually high number of terminal αGal-Gal moieties. One of the few reported structures of glycosylated mAbs containing these epitopes is the therapeutic antibody Cetuximab; however, unlike Cetuximab, one of the unliganded structures in S25-26 shows significant order in the glycan with appropriate electron density for nine residues. The elucidation of the three-dimensional structure of an αGal-containing N-linked glycan on a mAb variable heavy chain has potential clinical interest, as it has been implicated in allergic response in patients receiving therapeutic antibodies.

Recognition of the Thomsen-Friedenreich pancarcinoma carbohydrate antigen by a lamprey variable lymphocyte receptor

Variable lymphocyte receptors (VLRs) are leucine-rich repeat proteins that mediate adaptive immunity in jawless vertebrates. VLRs were recently shown to recognize glycans, such as the tumor-associated Thomsen-Friedenreich antigen (TFα; Galβ1-3GalNAcα), with a selectivity rivaling or exceeding that of lectins and antibodies. To understand the basis for TFα recognition by one such VLR (VLRB.aGPA.23), we measured thermodynamic parameters for the binding interaction and determined the structure of the VLRB.aGPA.23-TFα complex to 2.2 Å resolution. In the structure, four tryptophan residues form a tight hydrophobic cage encasing the TFα disaccharide that completely excludes buried water molecules. This cage together with hydrogen bonding of sugar hydroxyls to polar side chains explains the exquisite selectivity of VLRB.aGPA.23. The topology of the glycan-binding site of VLRB.aGPA.23 differs markedly from those of lectins or antibodies, which typically consist of long, convex grooves for accommodating the oligosaccharide. Instead, the TFα disaccharide is sandwiched between a variable loop and the concave surface of the VLR formed by the β-strands of the leucine-rich repeat modules. Longer oligosaccharides are predicted to extend perpendicularly across the β-strands, requiring them to bend to match the concavity of the VLR solenoid.

Ellipticine binds to a human telomere sequence: an additional mode of action as a putative anticancer agent?

Polyguanine sequences fold into G-quadruplex structures in the presence of monovalent cations. It is accepted that the telomeric DNA region consists of G-quadruplex structure. There are reports that potential G-quadruplex forming motifs are also present in the promoter region of some proto-oncogenes such as c-myc, c-kit, KRAS, etc. Small molecules with the potential to stabilize the telomeric DNA quadruplex have emerged as potential anticancer agents. We have studied the interaction of ellipticine, a putative anticancer agent from a plant source, with a human telomeric DNA sequence (H24). Spectroscopic and calorimetric studies indicate that the association of ellipticine with H24 is an entropically driven phenomenon with a 2:3 (H24:ellipticine) stoichiometry. Though ellipticine binding does not induce any major structural perturbation in H24, the association leads to formation of a complex with enhanced thermal stability. An assay with the telomerase repeat amplification protocol shows that ellipticine inhibits telomerase activity in MDAMB-231 breast cancer cell line extracts. This is the first report of the quadruplex binding ability of ellipticine. Using the results, we propose that along with DNA intercalation and/or topoisomerase II inhibition, interaction with the telomeric DNA region and the resultant inhibition of telomerase activity might be an additional mode of action for its anticancer property.

Molecular analysis of carbohydrate-antibody interactions: case study using a Bacillus anthracis tetrasaccharide

The process for selecting potent and effective carbohydrate antigens is not well-established. A combination of synthetic glycan microarray screening, surface plasmon resonance analysis, and saturation transfer difference NMR spectroscopy was used to dissect the antibody-binding surface of a carbohydrate antigen, revealing crucial binding elements with atomic-level detail. This analysis takes the first step toward uncovering the rules for structure-based design of carbohydrate antigens.

Molecular simulations of carbohydrates and protein-carbohydrate interactions: motivation, issues and prospects

The characterization of the 3D structure of oligosaccharides, their conjugates and analogs is particularly challenging for traditional experimental methods. Molecular simulation methods provide a basis for interpreting sparse experimental data and for independently predicting conformational and dynamic properties of glycans. Here, we summarize and analyze the issues associated with modeling carbohydrates, with a detailed discussion of four of the most recently developed carbohydrate force fields, reviewed in terms of applicability to natural glycans, carbohydrate-protein complexes and the emerging area of glycomimetic drugs. In addition, we discuss prospectives and new applications of carbohydrate modeling in drug discovery.

Structural basis for substrate selectivity in human maltase-glucoamylase and sucrase-isomaltase N-terminal domains

Human maltase-glucoamylase (MGAM) and sucrase-isomaltase (SI) are small intestinal enzymes that work concurrently to hydrolyze the mixture of linear alpha-1,4- and branched alpha-1,6-oligosaccharide substrates that typically make up terminal starch digestion products. MGAM and SI are each composed of duplicated catalytic domains, N- and C-terminal, which display overlapping substrate specificities. The N-terminal catalytic domain of human MGAM (ntMGAM) has a preference for short linear alpha-1,4-oligosaccharides, whereas N-terminal SI (ntSI) has a broader specificity for both alpha-1,4- and alpha-1,6-oligosaccharides. Here we present the crystal structure of the human ntSI, in apo form to 3.2 A and in complex with the inhibitor kotalanol to 2.15 A resolution. Structural comparison with the previously solved structure of ntMGAM reveals key active site differences in ntSI, including a narrow hydrophobic +1 subsite, which may account for its additional substrate specificity for alpha-1,6 substrates.

Sensitive and specific detection of the non-human sialic Acid N-glycolylneuraminic acid in human tissues and biotherapeutic products

BACKGROUND

Humans are genetically defective in synthesizing the common mammalian sialic acid N-glycolylneuraminic acid (Neu5Gc), but can metabolically incorporate it from dietary sources (particularly red meat and milk) into glycoproteins and glycolipids of human tumors, fetuses and some normal tissues. Metabolic incorporation of Neu5Gc from animal-derived cells and medium components also results in variable contamination of molecules and cells intended for human therapies. These Neu5Gc-incorporation phenomena are practically significant, because normal humans can have high levels of circulating anti-Neu5Gc antibodies. Thus, there is need for the sensitive and specific detection of Neu5Gc in human tissues and biotherapeutic products. Unlike monoclonal antibodies that recognize Neu5Gc only in the context of underlying structures, chicken immunoglobulin Y (IgY) polyclonal antibodies can recognize Neu5Gc in broader contexts. However, prior preparations of such antibodies (including our own) suffered from some non-specificity, as well as some cross-reactivity with the human sialic acid N-acetylneuraminic acid (Neu5Ac).

METHODOLOGY/PRINCIPAL FINDINGS

We have developed a novel affinity method utilizing sequential columns of immobilized human and chimpanzee serum sialoglycoproteins, followed by specific elution from the latter column by free Neu5Gc. The resulting mono-specific antibody shows no staining in tissues or cells from mice with a human-like defect in Neu5Gc production. It allows sensitive and specific detection of Neu5Gc in all underlying glycan structural contexts studied, and is applicable to immunohistochemical, enzyme-linked immunosorbent assay (ELISA), Western blot and flow cytometry analyses. Non-immune chicken IgY is used as a reliable negative control. We show that these approaches allow sensitive detection of Neu5Gc in human tissue samples and in some biotherapeutic products, and finally show an example of how Neu5Gc might be eliminated from such products, by using a human cell line grown under defined conditions.

CONCLUSIONS

We report a reliable antibody-based method for highly sensitive and specific detection of the non-human sialic acid Neu5Gc in human tissues and biotherapeutic products that has not been previously described.

Optimization of immobilized bacterial disaccharides for surface plasmon resonance imaging measurements of antibody binding

The interactions between proteins and immobilized carbohydrates are crucial to biological events such as cell signaling and immune response. The modification of surfaces with carbohydrates to create sensing platforms provides a pathway to study these interactions in a laboratory setting. In this work, a family of structurally related Salmonella disaccharide epitopes is immobilized on thin gold films in an array format to probe antibody binding with surface plasmon resonance (SPR) imaging. The disaccharides are modified with an alkyl thiol linker for facile immobilization to gold. Small differences in the stereochemistry of the immobilized, modified disaccharides are shown to greatly influence the binding of a monoclonal antibody. Specifically, binding is only observed to an immobilized abequose dideoxyhexose relative to a tyvelose or a paratose analogue. However, both the amount and relative strength of bound antibody depends on the distribution of disaccharide moieties in a mixed monolayer of the epitope and a nonbinding diluent molecule. We thoroughly characterize the mixed monolayers with a variety of techniques to understand the optimal density and distribution of the disaccharide for antibody capture. This work reinforces the importance of controlling the density of ligands at the interface for optimized surface based bioassays.

Thermodynamics and density of binding of a panel of antibodies to high-molecular-weight capsular polysaccharides

The interaction between antipolysaccharide (anti-PS) antibodies and their antigens was investigated by the use of isothermal titration calorimetry to determine the thermodynamic binding constant (K), the change in the enthalpy of binding (DeltaH), and the binding density (N) to high-molecular-weight PSs. From these values, the change in the entropy of binding (DeltaS) was calculated. The thermodynamic parameters of binding to high-molecular-weight capsular PSs are reported for two monoclonal antibodies (MAbs) with different specificities for meningococcal serogroup C PS, five MAbs specific for different pneumococcal serotypes, and the Fab fragments of two antipneumococcal MAbs. The K values were in the range of 10(6) to 10(7) M(-1), and these values were 1 to 2 orders of magnitude greater than the previously reported K values derived from antibody-oligosaccharide interactions. The DeltaH associated with binding was favorable for each MAb and Fab fragment. The DeltaS associated with binding was also generally favorable for both the MAbs and the Fab fragments, with the exception of the anti-serotype 14 MAb and its Fab fragment. N provides information regarding how densely MAbs or Fabs can bind along PS chains and, as expressed in terms of monosaccharides, was very similar for the seven MAbs, with an average of 12 monosaccharides per bound MAb. The value of N for each Fab was smaller, with five or seven monosaccharides per bound Fab. These results suggest that steric interactions between antibody molecules are a major influence on the values of N of high-affinity MAbs to capsular PSs.

Isothermal titration calorimetry reveals differential binding thermodynamics of variable region-identical antibodies differing in constant region for a univalent ligand

The classical view of immunoglobulin molecules posits two functional domains defined by the variable (V) and constant (C) regions, which are responsible for antigen binding and antibody effector functions, respectively. These two domains are thought to function independently. However, several lines of evidence strongly suggest that C region domains can affect the specificity and affinity of an antibody for its antigen (Ag), independent of avidity-type effects. In this study, we used isothermal titration calorimetry to investigate the thermodynamic properties of the interactions of four V region-identical monoclonal antibodies with a univalent peptide antigen. Comparison of the binding of IgG1, IgG2a, IgG2b, and IgG3 with a 12-mer peptide mimetic of Cryptococcus neoformans polysaccharide revealed a stoichiometry of 1.9-2.0 with significant differences in thermodynamic binding parameters. Binding of this peptide to the antibodies was dominated by favorable entropy. The interaction of these antibodies with biotinylated peptides manifested greater enthalpy than for native peptides indicating that biotin labeling affected the types of Ag-Ab complexes formed. Our results provide unambiguous thermodynamic evidence for the notion that the C region can affect the interaction of the V region with an Ag.

Diversity in specificity, abundance, and composition of anti-Neu5Gc antibodies in normal humans: potential implications for disease

Human heterophile antibodies that agglutinate animal erythrocytes are known to detect the nonhuman sialic acid N-glycolylneuraminic acid (Neu5Gc). This monosaccharide cannot by itself fill the binding site (paratope) of an antibody and can also be modified and presented in various linkages, on diverse underlying glycans. Thus, we hypothesized that the human anti-Neu5Gc antibody response is diverse and polyclonal. Here, we use a novel set of natural and chemoenzymatically synthesized glycans to show that normal humans have an abundant and diverse spectrum of such anti-Neu5Gc antibodies, directed against a variety of Neu5Gc-containing epitopes. High sensitivity and specificity assays were achieved by using N-acetylneuraminic acid (Neu5Ac)-containing probes (differing from Neu5Gc by one less oxygen atom) as optimal background controls. The commonest anti-Neu5Gc antibodies are of the IgG class. Moreover, the range of reactivity and Ig classes of antibodies vary greatly amongst normal humans, with some individuals having remarkably large amounts, even surpassing levels of some well-known natural blood group and xenoreactive antibodies. We purified these anti-Neu5Gc antibodies from individual human sera using a newly developed affinity method and showed that they bind to wild-type but not Neu5Gc-deficient mouse tissues. Moreover, they bind back to human carcinomas that have accumulated Neu5Gc in vivo. As dietary Neu5Gc is primarily found in red meat and milk products, we suggest that this ongoing antigen-antibody reaction may generate chronic inflammation, possibly contributing to the high frequency of diet-related carcinomas and other diseases in humans.

Detection of specific solvent rearrangement regions of an enzyme: NMR and ITC studies with aminoglycoside phosphotransferase(3')-IIIa

This work describes differential effects of solvent in complexes of the aminoglycoside phosphotransferase(3')-IIIa (APH) with different aminoglycosides and the detection of change in solvent structure at specific sites away from substrates. Binding of kanamycins to APH occurs with a larger negative DeltaH in H2O relative to D2O (DeltaDeltaH(H2O-D2O) < 0), while the reverse is true for neomycins. Unusually large negative DeltaCp values were observed for binding of aminoglycosides to APH. DeltaCp for the APH-neomycin complex was -1.6 kcal x mol(-1) x deg(-1). A break at 30 degrees C was observed in the APH-kanamycin complex yielding DeltaCp values of -0.7 kcal x mol(-1) x deg(-1) and -3.8 kcal x mol(-1) x deg(-1) below and above 30 degrees C, respectively. Neither the change in accessible surface area (DeltaASA) nor contributions from heats of ionization were sufficient to explain the large negative DeltaCp values. Most significantly, 15N-1H HSQC experiments showed that temperature-dependent shifts of the backbone amide protons of Leu 88, Ser 91, Cys 98, and Leu143 revealed a break at 30 degrees C only in the APH-kanamycin complex in spectra collected between 21 degrees C and 38 degrees C. These amino acids represent solvent reorganization sites that experience a change in solvent structure in their immediate environment as structurally different ligands bind to the enzyme. These residues were away from the substrate binding site and distributed in three hydrophobic patches in APH. Overall, our results show that a large number of factors affect DeltaCp and binding of structurally different ligand groups cause different solvent structure in the active site as well as differentially affecting specific sites away from the ligand binding site.

The contribution of conformational adjustments and long-range electrostatic forces to the CD2/CD58 interaction

CD2 is a T cell surface molecule that enhances T and natural killer cell function by binding its ligands CD58 (humans) and CD48 (rodents) on antigen-presenting or target cells. Here we show that the CD2/CD58 interaction is enthalpically driven and accompanied by unfavorable entropic changes. Taken together with structural studies, this indicates that binding is accompanied by energetically significant conformational adjustments. Despite having a highly charged binding interface, neither the affinity nor the rate constants of the CD2/CD58 interaction were affected by changes in ionic strength, indicating that long-range electrostatic forces make no net contribution to binding.

Oxalone and Lactone Moieties of Podophyllotoxin Exhibit Properties of Both the B and C Rings of Colchicine in Its Binding with Tubulin

Thermodynamics of podophyllotoxin binding to tubulin and its multiple points of attachment with tubulin has been studied in detail using isothermal titration calorimetry. The calorimetric enthalpy of the association of podophyllotoxin with tubulin is negative and occurs with a negative heat capacity change (DeltaC(p) = -2.47 kJ mol(-)(1) K(-)(1)). The binding is unique with a simultaneous participation of both hydrophobic and hydrogen-bonding forces with unfavorable negative entropic contribution at higher temperature, favored with an enthalpy-entropy compensation. Interestingly, the binding of 2-methoxy-5-(2',3',4'-trimethoxyphenyl)tropone (AC, a colchicine analogue without the B ring) with tubulin is enthalpy-favored. However, the podophyllotoxin-tubulin association depending upon the temperature of the reaction has a favorable entropic and enthalpic component, which resembles both B- and C-ring properties of colchicine. On the basis of the crystal structure of the podophyllotoxin-tubulin complex, distance calculations have indicated a possible interaction between threonine 179 of alpha-tubulin and the hydroxy group on the D ring of podophyllotoxin. To confirm the involvement of the oxalone moiety as well as the lactone ring of podophyllotoxin in tubulin binding, analogues of podophyllotoxin are synthesized with methoxy substitution at the 4' position of ring D along with its isomer and another analogue epimerized at ring E. From these results, involvement of oxalone as well as the lactone ring of the drug in a specific orientation inclusive of ring A is indicated for podophyllotoxin-tubulin binding. Therefore, podophyllotoxin, like colchicine, behaves as a bifunctional ligand having properties of both the B and C rings of colchicine by making more than one point of attachment with the protein tubulin.

Thermodynamic and kinetic characterization of ligand binding to the purine riboswitch aptamer domain

Riboswitches are cis-acting genetic regulatory elements found commonly in bacterial mRNAs that consist of a metabolite-responsive aptamer domain coupled to a regulatory switch. Purine riboswitches respond to intracellular concentrations of either adenine or guanine/hypoxanthine to control gene expression. The aptamer domain of the purine riboswitch contains a pyrimidine residue (Y74) that forms a Watson-Crick base-pairing interaction with the bound purine nucleobase ligand that discriminates between adenine and guanine. We sought to understand the structural basis of this specificity and the mechanism of ligand recognition by the purine riboswitch. Here, we present the 2,6-diaminopurine-bound structure of a C74U mutant of the xpt-pbuX guanine riboswitch, along with a detailed thermodynamic and kinetic analysis of nucleobase recognition by both the native and mutant riboswitches. These studies demonstrate clearly that the pyrimidine at position 74 is the sole determinant of purine riboswitch specificity. In addition, the mutant riboswitch binds adenine and adenine derivatives well compared with the guanine-responsive riboswitch. Under our experimental conditions, 2,6-diaminopurine binds the RNA with DeltaH=-40.3 kcal mol(-1), DeltaS=-97.6 cal mol(-1)K(-1), and DeltaG=-10.73 kcal mol(-1). A kinetic determination of the slow rate (0.15 x 10(5)M(-1)s(-1) and 2.1 x 10(5)mM(-1)s(-1) for 2-aminopurine binding the adenine-responsive mutant riboswitch and 7-deazaguanine-binding guanine riboswitch, respectively) of association under varying experimental conditions allowed us to propose a mechanism for ligand recognition by the purine riboswitch. A conformationally dynamic unliganded state for the binding pocket is stabilized first by the Watson-Crick base pairing between the ligand and Y74, and by the subsequent ordering of the J2/3 loop, enclosing the ligand within the three-way junction.

Modeling protein recognition of carbohydrates

We have a limited understanding of the details of molecular recognition of carbohydrates by proteins, which is critical to a multitude of biological processes. Furthermore, carbohydrate-modifying proteins such as glycosyl hydrolases and phosphorylases are of growing importance as potential drug targets. Interactions between proteins and carbohydrates have complex thermodynamics, and in general the specific positioning of only a few hydroxyl groups determines their binding affinities. A thorough understanding of both carbohydrate and protein structures is thus essential to predict these interactions. An atomic-level view of carbohydrate recognition through structures of carbohydrate-active enzymes complexed with transition-state inhibitors reveals some of the distinctive molecular features unique to protein-carbohydrate complexes. However, the inherent flexibility of carbohydrates and their often water-mediated hydrogen bonding to proteins makes simulation of their complexes difficult. Nonetheless, recent developments such as the parameterization of specific force fields and docking scoring functions have greatly improved our ability to predict protein-carbohydrate interactions. We review protein-carbohydrate complexes having defined molecular requirements for specific carbohydrate recognition by proteins, providing an overview of the different computational techniques available to model them.

Measurements of binding thermodynamics in drug discovery

Thermodynamics governs the process of biomolecular recognition. The steps of characterizing, understanding and exploiting binding thermodynamics have the potential to contribute to an improved rational drug design process that is more robust and reliable. It is only relatively recently that instrumentation capable of direct and full thermodynamic characterization has been improved, giving impetus to the application of thermodynamic measurements in drug discovery. This review highlights current instruments and methods that can be employed to measure binding thermodynamics and their use in studies of biomolecular recognition and drug discovery.

The interaction of carbohydrate-binding modules with insoluble non-crystalline cellulose is enthalpically driven

Natural cellulose exists as a composite of cellulose forms, which can be broadly characterized as crystalline or non-crystalline. The recognition of both of these forms of cellulose by the CBMs (carbohydrate-binding modules) of microbial glycoside hydrolases is important for the efficient natural and biotechnological conversion of cellulosic biomass. The category of CBM that binds insoluble non-crystalline cellulose does so with an affinity approx. 10-20-fold greater than their affinity for cello-oligosaccharides and/or soluble polysaccharides. This phenomenon has been assumed to originate from the effects of changes in configurational entropy upon binding. The loss of configurational entropy is thought to be less profound upon binding to conformationally restrained insoluble non-crystalline cellulose, resulting in larger free energies of binding. However, using isothermal titration calorimetry, it is shown that this is not the case for the high-affinity interactions of CcCBM17 (the family 17 CBM from EngF of Clostridium cellulovorans) and BspCBM28 (the family 28 CBM from Cel5A of Bacillus species 1139) with regenerated cellulose, an insoluble preparation of primarily non-crystalline cellulose. The enhanced free energy of binding of non-crystalline cellulose relative to cello-oligosaccharides is by virtue of improved enthalpy, not entropy.

Energetics of carbohydrate binding to Momordica charantia (bitter gourd) lectin: an isothermal titration calorimetric study

Physico-chemical and carbohydrate binding studies have been carried out on the Momordica charantia (bitter gourd) seed lectin (MCL). The lectin activity is maximal in the pH range 7.4-11.0, but decreases steeply below pH 7.0. The lectin activity is mostly unaffected in the temperature range 4-50 degrees C, but a sharp decrease is seen between 50 and 60 degrees C, which could be correlated to changes in the structure of the protein as seen by circular dichroism and fluorescence spectroscopy. Isothermal titration calorimetric studies show that the tetrameric MCL binds two sugar molecules and the binding constants (Kb), determined at 288.15 K, for various saccharides were found to vary between 7.3 x 10(3) and 1.52 x 10(4)M(-1). The binding reactions for all the saccharides investigated were essentially enthalpy driven, with the binding enthalpies (DeltaHb) at 288.15 K being in the range of -50.99 and -43.39 kJ mol(-1), whereas the contribution to the binding reaction from the entropy of binding was negative, with values of binding entropy (DeltaSb) ranging between -99.2 and -72.0 J mol(-1)K(-1) at 288.15 K. Changes in heat capacity (DeltaCp) for the binding of disaccharides, lactose and lactulose, were significantly larger in magnitude than those obtained for the monosaccharides, methyl-beta-D-galactopyranoside, and methyl-alpha-D-galactopyranoside, and could be correlated reasonably well with the surface areas of these ligands. Enthalpy-entropy compensation was observed for all the sugars studied, suggesting that water structure plays an important role in the overall binding reaction. CD spectroscopy indicates that carbohydrate binding does not lead to significant changes in the secondary and tertiary structures of MCL, suggesting that the carbohydrate binding sites on this lectin are mostly preformed.

Thermodynamics of binding of D-galactose and deoxy derivatives thereof to the L-arabinose-binding protein

We report the thermodynamics of binding of d-galactose and deoxy derivatives thereof to the arabinose binding protein (ABP). The "intrinsic" (solute-solute) free energy of binding DeltaG degrees (int) at 308 K for the 1-, 2-, 3-, and 6-hydroxyl groups of galactose is remarkably constant (approximately -30 kJ/mol), despite the fact that each hydroxyl group subtends different numbers of hydrogen bonds in the complex. The substantially unfavorable enthalpy of binding (approximately 30 kJ/mol) of 1-deoxygalactose, 2-deoxygalactose, and 3-deoxygalactose in comparison with galactose, cannot be readily accounted for by differences in solvation, suggesting that solute-solute hydrogen bonds are enthalpically significantly more favorable than solute-solvent hydrogen bonds. In contrast, the substantially higher affinity for 2-deoxygalactose in comparison with either 1-deoxygalactose or 3-deoxygalactose derives from differences in the solvation free energies of the free ligands.

Thermodynamic origin of the chiral recognition of tryptophan on teicoplanin and teicoplanin aglycone stationary phases

The D-, L-tryptophan binding and the chiral recognition properties of the teicoplanin and teicoplanin aglycone (TAG) chiral stationary phase (CSPs) were compared at various column temperatures. The solute adsorption isotherms (bi-Langmuir model) were determined for both the two CSPs using the perturbation method. It was demonstrated that the sugar units were involved in the reduction of the apparent enantioselectivity through two phenomena: (i) the inhibition of some enantioselective contacts with low-affinity binding regions of the aglycone and (ii) a decrease in the stereoselective properties of the aglycone high-affinity binding pocket. The phenomenon (ii) was governed by both a decrease in the ratio of the enantiomer adsorption constant and a strong reduction of the site accessibility for D- and L-tryptophan. In addition, a temperature effect study was performed to investigate the chiral recognition mechanism at the aglycone high-affinity pocket. An enthalpy-entropy compensation analysis derived from the Grunwald model as well as the comparison with the literature data demonstrated that the enantioselective binding mode was dependent on an interface dehydration process. The change in the enantioselective process observed between the TAG and teicoplanin CSP was characterized by a difference of ca. 2-3 ordered water molecules released from the species interface.

Evaluation of fluorescence-based thermal shift assays for hit identification in drug discovery

The fluorescence-based thermal shift assay is a general method for identification of inhibitors of target proteins from compound libraries. Using an environmentally sensitive fluorescent dye to monitor protein thermal unfolding, the ligand-binding affinity can be assessed from the shift of the unfolding temperature (Delta Tm) obtained in the presence of ligands relative to that obtained in the absence of ligands. In this article, we report that the thermal shift assay can be conducted in an inexpensive, commercially available device for temperature control and fluorescence detection. The binding affinities obtained from thermal shift assays are compared with the binding affinities measured by isothermal titration calorimetry and with the IC(50) values from enzymatic assays. The potential pitfalls in the data analysis of thermal shift assays are also discussed.

Iminosugar Glycosidase Inhibitors: Structural and Thermodynamic Dissection of the Binding of Isofagomine and 1-Deoxynojirimycin to β-Glucosidases

The design and synthesis of transition-state mimics reflects the growing need both to understand enzymatic catalysis and to influence strategies for therapeutic intervention. Iminosugars are among the most potent inhibitors of glycosidases. Here, the binding of 1-deoxynojirimycin and (+)-isofagomine to the "family GH-1" beta-glucosidase of Thermotoga maritima is investigated by kinetic analysis, isothermal titration calorimetry, and X-ray crystallography. The binding of both of these iminosugar inhibitors is driven by a large and favorable enthalpy. The greater inhibitory power of isofagomine, relative to 1-deoxynojirimycin, however, resides in its significantly more favorable entropy; indeed the differing thermodynamic signatures of these inhibitors are further highlighted by the markedly different heat capacity values for binding. The pH dependence of catalysis and of inhibition suggests that the inhibitory species are protonated inhibitors bound to enzymes whose acid/base and nucleophile are ionized, while calorimetry indicates that one proton is released from the enzyme upon binding at the pH optimum of catalysis (pH 5.8). Given that these results contradict earlier proposals that the binding of racemic isofagomine to sweet almond beta-glucosidase was entropically driven (Bülow, A. et al. J. Am. Chem. Soc. 2000, 122, 8567-8568), we reinvestigated the binding of 1-deoxynojirimycin and isofagomine to the sweet almond enzyme. Calorimetry confirms that the binding of isofagomine to sweet almond beta-glucosidases is, as observed for the T. maritima enzyme, driven by a large favorable enthalpy. The crystallographic structures of the native T. maritima beta-glucosidase, and its complexes with isofagomine and 1-deoxynojirimycin, all at approximately 2.1 A resolution, reveal that additional ordering of bound solvent may present an entropic penalty to 1-deoxynojirimycin binding that does not penalize isofagomine.

The influence of competing factors on the endothermic nature of bile salt binding by cationic adsorbent

The binding of sodium chenodeoxycholate, a hydrophobic bile salt, by a polyacrylamide resin with N,N,N-trimethylammonium dodecyl chloride (QPDA12) pendant groups was studied in the presence of elevated concentrations of competing anions. The equilibrium bile salt concentration was determined from HPLC data for a range of initial bile salt concentrations. Binding constants extracted from the fit of the isotherms to the Langmuir equation were obtained for data for temperatures from 24 to 60°C. A reduction in chenodeoxycholate binding affinity was observed in comparison with that reported previously at lower overall anion concentrations. A van't Hoff analysis of the data revealed that a decrease in the favourable entropy of binding was responsible for this reduced affinity. Decomposition of the binding constants into specific contributions for the competing buffer anions and the chenodeoxycholate anion returned positive values of ΔH° and ΔS° for the binding of both the competing ions and bile salt. These findings reveal that the weakening of the bile salt binding that occurs with added salt is not due to Debye screening of an attraction between ions of opposite charge. A binding mechanism consistent with the thermodynamics observed was proposed in which the dominant role of hydration changes that occur on binding provides the principal driving force for the process.Key words: bile salts, endothermic binding, polymer sorbents, solvent randomization.

Synthesis of threeSalmonellaepitopes for biosensor studies of carbohydrate–antibody interactions

Disaccharides 1-3 corresponding to the antigenic determinants of Salmonella serotypes A, B, and D1were synthesized in a form suited for use in biosensors. The disaccharide determinants each contain a unique 3,6-dideoxyhexose, namely abequose (3,6-dideoxy-D-xylo-hexose), paratose (3,6-dideoxy-D-ribohexose), and tyvelose (3,6-dideoxy-D-arabino-hexose), are α-linked to the 3-position of D-mannopyranose. The disaccharides were further derivatized with a linear aglycon that has a terminal amino group, and can be readily coupled to pertinent chains carrying a terminal thiol for the construction of self-assembled monolayers (SAMs). Efficient routes that employed a single 3,6-dideoxygenation step were developed for the synthesis of paratoside 15 and tyveloside 22.Key words: Salmonella O-antigens, lipopolysaccharide, abequose, paratose, tyvelose, 3,6-dideoxyhexose, deoxygenation, glycoside tethers, immobilization via pentenyl glycosides.

Oligosaccharide recognition by antibodies: Synthesis and evaluation of talose oligosaccharide analogues

A series of monosaccharide (4–6), disaccharide (3,7–12), and trisaccharide (13–15) analogs of the native ligand 2, which fills the binding site of monoclonal antibody Se 155.4, have been synthesized and their bioactivity measured by solid- and solution-phase assays. The syntheses of disaccharide analogs sought to replace galactose by various alkyl groups at the O-2 position of mannose. The activity of one of these O-2 alkyl analogs was 75% of that observed for the trisaccharide and points to only weak net bonding between the solvent exposed galactose residue and the antibody binding site. The synthesis of talose analogs 13 and 14, where the mannose or galactose residues of 2 were replaced by talose produced ligands with activities from one-third to one-half of that seen for the native ligand 2. These activity changes did not exhibit discernable correlations with the ability of talose to disrupt water of solvation.Key words: abequose, 3,6-dideoxy-D-xylo-hexose, talose disaccharide and trisaccharide, antibody oligosaccharide interactions, molecular recognition of carbohydrates, water in antibody complexes, Salmonella LPS, monoclonal antibody Se 155.4, bacterial O-antigen.

Carbohydrate-antibody interactions by NMR for a13C-labelled disaccharide ligand

Incorporation of a13C label into a carbohydrate ligand, methyl 3-O-(3,6-dideoxy-α-D-xylohexopyranosyl)-2-O-methyl-α-D-mannopyranoside permitted by NMR spectroscopy the study of its binding to the Fab from a monoclonal antibody, Se 155-4. The signals of the free and bound form were observed in the13C spectrum of the carbohydrate-protein complex. The dissociation rate constants were consequently determined by full lineshape analysis of the13C spectrum. Comparison with simplified analyses relying only on the linewidth of the1H and13C signals of the free ligand were made and the justifications of underlying assumptions used in these analyses were discussed. For1H NMR, the protein resonances were purged with a13C filter to observe only the ligand resonances and NOEs between the ligand and the protein.Key words: carbohydrate, binding, NMR, C-13 label, chemical exchange.

Differential oligosaccharide recognition by evolutionarily-related beta-1,4 and beta-1,3 glucan-binding modules

Enzymes active on complex carbohydrate polymers frequently have modular structures in which a catalytic domain is appended to one or more carbohydrate-binding modules (CBMs). Although CBMs have been classified into a number of families based upon sequence, many closely related CBMs are specific for different polysaccharides. In order to provide a structural rationale for the recognition of different polysaccharides by CBMs displaying a conserved fold, we have studied the thermodynamics of binding and three-dimensional structures of the related family 4 CBMs from Cellulomonas fimi Cel9B and Thermotoga maritima Lam16A in complex with their ligands, beta-1,4 and beta-1,3 linked gluco-oligosaccharides, respectively. These two CBMs use a structurally conserved constellation of aromatic and polar amino acid side-chains that interact with sugars in two of the five binding subsites. Differences in the length and conformation of loops in non-conserved regions create binding-site topographies that complement the known solution conformations of their respective ligands. Thermodynamics interpreted in the light of structural information highlights the differential role of water in the interaction of these CBMs with their respective oligosaccharide ligands.

Thermodynamics of aminoglycoside and acyl-coenzyme A binding to the Salmonella enterica AAC(6')-Iy aminoglycoside N-acetyltransferase

Kinetic and mechanistic studies on the chromosomally encoded aminoglycoside 6'-N-acetyltransferase, AAC(6')-Iy, of Salmonella enterica that confers resistance toward aminoglycosides have been previously reported [Magnet et al. (2001) Biochemistry 40, 3700-3709]. In the present study, equilibrium binding and the thermodynamic parameters of binding of aminoglycosides and acyl-coenzyme A derivatives to AAC(6')-Iy and of two mutants, C109A and the C109A/C70A double mutant, have been studied using fluorescence spectroscopy and isothermal titration calorimetry (ITC). Association constants for different aminoglycosides varied greatly (4 x 10(4)-150 x 10(4)) while the association constants of several acyl-coenzyme A derivatives were similar (3.2 x 10(4)-4.5 x 10(4)). The association constants and van't Hoff enthalpy changes derived from intrinsic protein fluorescence changes were in agreement with independently measured values from isothermal titration calorimetry studies. Binding of both aminoglycosides and acyl-coenzyme A derivatives is strongly enthalpically driven and revealed opposing negative entropy changes, resulting in enthalpy-entropy compensation. The acetyltransferase exhibited a temperature-dependent binding of tobramycin with a negative heat capacity value of 410 cal mol(-1) K(-1). Isothermal titration studies of acetyl-coenzyme A and tobramycin binding to mutant forms of the enzyme indicated that completely conserved C109 does not play any direct role in the binding of either of the substrates, while C70 is directly involved in aminoglycoside binding. These results are discussed and compared with previous steady-state kinetic studies of the enzyme.

van't Hoff and calorimetric enthalpies II: effects of linked equilibria

The complexity of binding reactions, including the linkage with other equilibria, is becoming increasingly apparent in biological processes such as signal transduction. Understanding these interactions requires obtaining thermodynamic profiles for each of the equilibria that occur in a binding event. Concern has been raised as to whether linked equilibria contribute differently to thermodynamics, such as DeltaH degrees and DeltaC(p), obtained from calorimetric and van't Hoff methods. We have previously shown that linked equilibria do not contribute differently to the van't Hoff and calorimetrically determined DeltaH degrees for processes such as linked folding or hydration. Here, examples of proton and ion linkage are examined. We show that there is no reason to expect the calorimetric and van't Hoff DeltaH degrees to be different, even without prior knowledge of the presence or absence of linked equilibria, as long as the system is permitted to equilibrate. However, it is possible to create experimental scenarios that result in and discrepancies. Furthermore, it is found that the presence of linked equilibria in all cases can result in "nonconventional" DeltaH degrees and DeltaC(p) profiles, making data analysis nontrivial.

Temperature dependence of thermodynamic properties for DNA/DNA and RNA/DNA duplex formation

A clear difference in the enthalpy changes derived from spectroscopic and calorimetric measurements has recently been shown. The exact interpretation of this deviation varied from study to study, but it was generally attributed to the non-two-state transition and heat capacity change. Although the temperature-dependent thermodynamics of the duplex formation was often implied, systemic and extensive studies have been lacking in universally assigning the appropriate thermodynamic parameter sets. In the present study, the 24 DNA/DNA and 41 RNA/DNA oligonucleotide duplexes, designed to avoid the formation of hairpin or slipped duplex structures and to limit the base pair length less than 12 bp, were selected to evaluate the heat capacity changes and temperature-dependent thermodynamic properties of duplex formation. Direct comparison reveals that the temperature-independent thermodynamic parameters could provide a reasonable approximation only when the temperature of interest has a small deviation from the mean melting temperature over the experimental range. The heat capacity changes depend on the base composition and sequences and are generally limited in the range of -160 to approximately -40 cal.mol-1.K-1 per base pair. In contrast to the enthalpy and entropy changes, the free energy change and melting temperature are relatively insensitive to the heat capacity change. Finally, the 16 NN-model free energy parameters and one helix initiation at physiological temperature were extracted from the temperature-dependent thermodynamic data of the 41 RNA/DNA hybrids.

Isothermal titration calorimetry in drug discovery

Isothermal titration calorimetry (ITC) follows the heat change when a test compound binds to a target protein. It allows precise measurement of affinity. The method is direct, making interpretation facile, because there is no requirement for competing molecules. Titration in the presence of other ligands rapidly provides information on the mechanism of action of the test compound, identifying the intermolecular complexes that are relevant for structure-based design. Calorimetry allows measurement of stoichiometry and so evaluation of the proportion of the sample that is functional. ITC can characterize protein fragments and catalytically inactive mutant enzymes. It is the only technique which directly measures the enthalpy of binding (delta H degree). Interpretation of delta H degree and its temperature dependence (delta Cp) is usually qualitative, not quantitative. This is because of complicated contributions from linked equilibria and a single change in structure giving modification of several physicochemical properties. Measured delta H degree values allow characterization of proton movement linked to the association of protein and ligand, giving information on the ionization of groups involved in binding. Biochemical systems characteristically exhibit enthalpy-entropy compensation where increased bonding is offset by an entropic penalty, reducing the magnitude of change in affinity. This also causes a lack of correlation between the free energy of binding (delta G degree) and delta H degree. When characterizing structure-activity relationships (SAR), most groups involved in binding can be detected as contributing to delta H degree, but not to affinity. Large enthalpy changes may reflect a modified binding mode, or protein conformation changes. Thus, delta H degree values may highlight a potential discontinuity in SAR, so that experimental structural data are likely to be particularly valuable in molecular design.

Binding specificity and thermodynamics of a family 9 carbohydrate-binding module from Thermotoga maritima xylanase 10A

The C-terminal family 9 carbohydrate-binding module of xylanase 10A from Thermotoga maritima (CBM9-2) binds to amorphous cellulose, crystalline cellulose, and the insoluble fraction of oat spelt xylan. The association constants (K(a)) for adsorption to insoluble polysaccharides are 1 x 10(5) to 3 x 10(5) M(-1). Of the soluble polysaccharides tested, CBM9-2 binds to barley beta-glucan, xyloglucan, and xylan. CBM9-2 binds specifically to the reducing ends of cellulose and soluble polysaccharides, a property that is currently unique to this CBM. CBM9-2 also binds glucose, xylose, galactose, arabinose, cellooligosaccharides, xylooligosaccharides, maltose, and lactose, with affinities ranging from 10(3) M(-1) for monosaccharides to 10(6) M(-1) for disaccharides and oligosaccharides. Cellooligosaccharides longer than two glucose units do not bind with improved affinity, indicating that cellobiose is sufficient to occupy the entire binding site. In general, the binding reaction is dominated by favorable changes in enthalpy, which are partially compensated by unfavorable entropy changes.