A.-L.
Noresson
a,
O.
Aurelius
b,
C. T.
Öberg
a,
O.
Engström
a,
A. P.
Sundin
a,
M.
Håkansson
c,
O.
Stenström
d,
M.
Akke
d,
D. T.
Logan
bc,
H.
Leffler
e and
U. J.
Nilsson
*a
aCentre for Analysis and Synthesis, Department of Chemistry, Lund University, Box 124, SE-221 00 Lund, Sweden. E-mail: ulf.nilsson@chem.lu.se
bSection for Biochemistry and Structural Biology, Center for Molecular Protein Science, Department of Chemistry, Lund University, Box 124, SE-221 00 Lund, Sweden
cSARomics Biostructures AB, Medicon Village, SE-223 81 Lund, Sweden
dBiophysical Chemistry, Center for Molecular Protein Science, Department of Chemistry, Lund University, Box 124, SE-221 00 Lund, Sweden
eDepartment of Laboratory Medicine, Section MIG, Lund University, Sölvegatan 23, SE-223 62, Lund, Sweden
First published on 4th December 2017
We investigated galectin-3 binding to 3-benzamido-2-O-sulfo-galactoside and -thiodigalactoside ligands using a combination of site-specific mutagenesis, X-ray crystallography, computational approaches, and binding thermodynamics measurements. The results reveal a conformational variability in a surface-exposed arginine (R144) side chain in response to different aromatic C3-substituents of bound galactoside-based ligands. Fluorinated C3-benzamido substituents induced a shift in the side-chain conformation of R144 to allow for an entropically favored electrostatic interaction between its guanidine group and the 2-O-sulfate of the ligand. By contrast, binding of ligands with non-fluorinated substituents did not trigger a conformational change of R144. Hence, a sulfate–arginine electrostatic interaction can be tuned by the choice of ligand C3-benzamido structures to favor specific interaction modes and geometries. These results have important general implications for ligand design, as the proper choice of arginine–aromatic interacting partners opens up for ligand-controlled protein conformation that in turn may be systematically exploited in ligand design.
The β-D-galactopyranoside-binding protein galectin-3 belongs to a family of lectins that participates in a number of biologically important processes, of which several are pathological under certain circumstances.6–10 An advantageous feature of galectin-3 is the availability of ample structural and thermodynamic data for galectin-3-ligand complexes, which makes galectin-3 a suitable model protein for investigations of design strategies based on protein–ligand conformational adaptability. Structural analysis of galectin-3 in complex with a LacNAc derivative carrying a 4-methoxy-2,3,5,6-tetrafluoro benzamide at C3 revealed that this ligand benzamido moiety induces a large conformational change in the galectin-3C residue R144, moving it away from a water-mediated salt bridge with N148 on the surface of galectin-3 and positioning the benzamide between the protein backbone and R144 (ref. 5) to form an arginine–arene interaction (Fig. 1A). Later, a structural analysis of galectin-3 complex with a LacNAc derivative carrying an unsubstituted benzamide at C3 revealed that R144 maintains its interaction with N148, while the benzamide moiety is positioned on top of R144,4 again forming an arginine–arene interaction but with an alternative complex geometry (Fig. 1B). Furthermore, galectin-3 conformational dynamics and entropies were affected differently by 4-methoxy-2,3,5,6-tetrafluorobenzamide and unsubstituted benzamide, respectively.4 Hence, the arene properties of the ligand determine the interaction mode and the protein conformation and dynamics. The conformational variability of galectin-3 R144 subsequently led to the hypothesis that compounds with electron-rich groups at the galactose 2-O position could form beneficial interactions with the R144 side-chain conformation properly controlled by the ligand C3 aromatic amide moiety, as shown in Fig. 1A. Indeed, electron-rich substituents at galactose 2-O were evidenced to be affinity-enhancing and this was suggested to be due to interactions with R144 (ref. 11 and 12) (Fig. 1C).
Fig. 1 X-ray structures4,5 of galectin-3 C-terminal domain in complex with ligands derived from LacNAc carrying (A) a 4-methoxy-2,3,5,6-tetrafluoro-benzamide 2 or (B) a benzamide 1 at C3′,5 revealing the conformational variability in R144. Galactose O2 is indicated with *. (C) Galectin-3 ligands with a sulfate at the 2-O-position were beneficial for the affinity. This was hypothesized to be due to interactions with the positively charged amino acids.11,12 (D) 3,3-Dibenzamido-thiodigalactosides high-affinity ligands of galectin-3 proposed to be due to benzamide interactions with both R144 and R186.13,14 |
3,3′-Dibenzamido-thiodigalactosides have been shown to be galectin-3 ligands with much higher affinity than the simpler C3-benzamido LacNAc derivatives,13,14 which was explained by two benzamide–arginine interactions involving not only R144 but also R186 (Fig. 1D). The high affinity of the 3,3-dibenzamido-thiodigalactosides for galectin-3 provided an impetus to investigate electron-rich sulfate substituents at one of the galactose 2-O positions of these structures. Herein, we report synthesis, structural, and thermodynamic analysis of 2-sulfo-galactoside and -thiodigalactoside derivatives carrying C3-amides that either induce or do not induce a conformational change in R144 on the surface of galectin-3, as well as the exploitation of such structures to further control R144 conformation and interactions.
2-O- | K d | ΔΔG | |
---|---|---|---|
7 | –H | 1300 ± 390 | |
9 | –SO3– | 65 ± 16 | −7.4 |
10 | –H | 620 ± 160 | |
11 | –SO3– | 370 ± 90 | −1.3 |
Fig. 2 Proposed R144 interaction modes for compounds 7, 9, 10, and 11 in complex with galectin-3 (molecular modelling was performed with the MMFFs force field with water implemented in MacroModel version 9.1, Schrödinger, LLC, New York). The galectin-3 structure and ligand and R144 starting conformations were taken from PDB id 1KJL for 10 and 11 and 1KJR for 7 and 9. |
The introduction of one 2-O-sulfate group on the bis-3,3′-(3-methoxy-benzamido)-thiodigalactoside to give compound 19 did not result in a large affinity improvement over 15 (Table 2). The introduction of a 2-O-sulfate on the bis-3,3′-(4-methoxy-2,3,5,6-tetrafluoro-benzamido)-thiodigalactoside scaffold 20 resulted in a modest affinity improvement over the non-sulfated analog 16, which was only slightly greater in relative terms than that of the bis-3,3′-(3-methoxy-benzamido)-thiodigalactosides 15 and 19 (Table 2). Furthermore, the bis-3,3′-(3-methoxy-benzamido)-thiodigalactosides (15 and 19) displayed a higher affinity for galectin-3 than the fluorinated ones (16 and 20), for both the sulfated and the non-sulfated analogs. Thus, the possible ion pair formed between the 2-O-sulfated bis-3,3′-(4-methoxy-2,3,5,6-tetrafluoro-benzamido)-thiodigalactoside 20 and R144 unexpectedly contributes to only a marginal gain in affinity as compared to the corresponding monosaccharide 9. Furthermore, a higher galectin-3 affinity was observed for the 4-methoxy-2,3,5,6-tetrafluoro-benzamido-LacNAc derivative 2 (Fig. 1, Kd 0.88 μM) than the corresponding 3-methoxy-benzamido-LacNAc derivative (Kd 2.5 μM), which suggests that the added second benzamido moiety in 16 and 20, that is proposed to interact with R186, may negatively influence the affinity for galectin-3. Furthermore, mutant studies support the importance of R144 interacting with both types of benzamide substituents, as a clear decrease in affinity of compounds 15 and 16 is seen for the R144S18 and R144K mutants (Table 2). The sulfated compounds 19 and 20 show a relatively larger loss in affinity for the R144S mutant, but less so for the R144K mutant. This result can be explained by the similar length and charge of lysine and arginine side chains, which both allow for attractive interactions with the sulfate groups of 19 and 20 that cannot occur in the R144S mutant. In case of the R186S18 and R186K mutants, all compounds 15, 16, 19, and 20 have dramatically lower affinity (Table 2), which can be explained by the loss of the favourable interaction between the R186 guanidine group and the benzamide moiety of the ligands.19 By contrast, the limited reduction in affinity observed for the corresponding mutations of R144 (R144S and R144K) might suggest that guanidine–aryl interactions with this residue contributes to free energy of binding to a lesser extent than guanidine–aryl interactions of R186 do. In support of this hypothesis, we have previously observed that R144 undergoes exchange between alternative conformations in all ligand-bound states studied by NMR spectroscopy, whereas R186 is well ordered.4 Altogether, the mutant studies suggest an R144-benzamide interaction in 15, 16, 19, and 20, and the presence of a 2-O-sulfate does affect mutant affinities similarly for both compound pairs (comparing 15 and 19, or 16 and 20).
wt | R144S | R144K | R186S | R186K | |
---|---|---|---|---|---|
15 | 110 ± 20 | 215 ± 28 | 220 ± 73 | 42000 ± 7800 | 26000 ± 7900 |
19 | 57 ± 8 | 550 ± 76 | 76 ± 1 | 33000 ± 5300 | 15000 ± 2200 |
16 | 500 ± 51 | 1800 ± 280 | 1300 ± 210 | 79000 ± 1500 | 53000 ± 1300 |
20 | 210 ± 40 | 2900 ± 320 | 490 ± 46 | 29000 ± 2900 | 15000 ± 2100 |
Fig. 3 (A) X-ray crystal structures of the galectin-3 C-terminal domain in complex with lactose,20 (B) overlaid X-ray crystal structures of the complexes with 15 (pale cyan), 16 (grey), 19 (light blue), and 20 (light orange) depicting hydrogen bond patterns identical to that of bound lactose. Ligand hydrogen bonding oxygens/hydroxyls are red. |
The R186 side chain forms in-plane ion-pairs with E165 and E184, and does not change conformation upon binding of 15, 16, 19, or 20. Face-to-face stacking of one of the benzamides of 15 or 19 onto the R186 side chain is identical for both compounds and strongly supports the hypothesis that such face-to-face stacking of R186 against a ligand phenyl ring is beneficial and contributes to the large affinity gain for bis-amido-thiodigalactosides (e.g.15) over 3′-benzamido-LacNAc derivatives.13,14,19 Alignment of the four complexes of 15, 16, 19, and 20 with galectin-3 reveals that fluorination does not only affect the orientation of R144, but also alters the stacking modes of the benzamido moieties of 16 and 20 onto R186 (Fig. 5A), as compared to the corresponding ones in the complexes of the non-fluorinated 15 and 19 (Fig. 5B). The fluorinated 16 and 20 appears have less optimal π-system overlap with R186 as 16 has the benzamide tilted out of the plane above R186 (Fig. 5B) and 20 is shifted to the edge of the R186 plane (Fig. 5C), which suggests possible explanations for the unexpected relative lower affinity of 16 and 20 for galectin-3 compared to 15 and 19, respectively.
ΔG (kJ mol−1) | ΔH (kJ mol−1) | −TΔS (kJ mol−1) | n | K d (nM) | |
---|---|---|---|---|---|
a ITC data were analyzed by global optimization of multiple data sets. Each compound was fitted using two approaches: the first row shows the results of fitting all parameters to the ITC data, while the second row shows the results obtained when keeping Kd fixed at the value obtained from the competitive FP assay (Table 2). b The relative error in n is less than 3% in all cases. | |||||
15 | −39.40 ± 0.05 | −70.0 ± 0.2 | 30.7 ± 0.2 | 0.80 | 126 ± 2 |
15 | −39.7 ± 0.2 | −69.3 ± 0.5 | 29.6 ± 0.5 | 0.80 | |
19 | −39.9 ± 0.1 | −55.4 ± 0.3 | 15.6 ± 0.4 | 1.01 | 103 ± 56 |
19 | −41.3 ± 0.3 | −53.2 ± 0.4 | 11.9 ± 0.5 | 1.01 | |
16 | −35.8 ± 0.5 | −54.0 ± 3.0 | 18.0 ± 3.0 | 0.89 | 604 ± 129 |
16 | −36.3 ± 0.3 | −52.0 ± 1.0 | 16.0 ± 1.0 | 0.89 | |
20 | −37.9 ± 0.2 | −46.8 ± 0.6 | 8.8 ± 0.6 | 1.06 | 261 ± 19 |
20 | −38.5 ± 0.5 | −46.0 ± 1.0 | 7.2 ± 1.2 | 1.06 |
Fig. 6 Thermodynamic parameters of binding 15, 16, 19, and 20 to galectin-3, as determined from ITC by fitting all parameters. |
The fluorinated compounds (16 and 20) bind less strongly and with less favourable enthalpy, but also less entropic penalty, than their non-fluorinated counterparts (15 and 19): ΔΔH = 16 kJ mol−1 (16vs.15) and 8 kJ mol−1 (20vs.19) and −TΔΔS = −13 kJ mol−1 (15 and 19) and −7 kJ mol−1 (20vs.19). These results can be rationalized in part by the decreased π-electron density of phenyl rings upon fluorine substitution, which decreases the strength of the cation–π interactions between the guanidine groups of R144 and R186 and the phenyl rings, and by the poor solvation of fluorinated aromatic rings that make binding to the protein relatively more favourable in terms of entropy compared to the non-fluorinated analogues. The former explanation is in line with the structural differences shown in Fig. 5C.
The somewhat higher affinities of the 2-O-sulfated 19 and 20, as compared to the non-sulfated 15 and 16, were also confirmed by ITC, with the effect being slightly greater for the fluorinated ligand, as expected. The affinity increase upon sulfation (15vs.19 and 16vs.20, respectively) is driven by entropy alone, whereas the change in binding enthalpy is unfavourable. These changes in binding enthalpy and entropy induced by the sulfate group, viz. ΔΔH > 0 and −TΔΔS < 0, are expected because the attractive Coulomb interaction between opposite charges (i.e., the guanidine and sulfate groups) is entirely due to favourable entropy and involves unfavourable enthalpy.27,28 We note that the sulfate group of 19 had a much greater influence on both enthalpy and entropy than did the sulfate group of 20. However, analyzing the change in binding free energy (ΔΔG) upon introduction of a 2-O-sulfate to give structures 16 and 20, it is clear that the 2-O-sulfate group of 20 achieves a higher gain in free energy than the 2-O-sulfate group of 19. A useful point of reference is the ratio ΔΔS/ΔΔH,27,28 which should be close to 1.35ΔΔG/(−0.35ΔΔG) = −3.85 in the hypothetical case that ΔΔG is caused entirely by electrostatic interactions. We find that ΔΔS/ΔΔH = −1.3 for the comparison of 20vs.16, while ΔΔS/ΔΔH = −1.0 for the comparison of 19vs.15, thus suggesting that electrostatic effects play a greater role in 20 than in 19, as expected from the crystal structures. However, it should be kept in mind that this simplified analysis neglects all other contributions to binding, such as differential desolvation of the ligands, different interactions in the region surrounding R186, and conformational entropy, all of which surely contribute to ΔΔG.
Taken together, this study demonstrates that proper choice of aromatic arginine side-chain interacting partners opens up for ligand-controlled arginine conformations that in turn may be systematically exploited in design of synthetic ligands that induce arginine side chain conformations that form beneficial interactions with electron-rich (anionic) ligand functionalities.
Footnotes |
† Electronic supplementary information (ESI) available: Synthesis experimental procedures and physical data for compounds. Data processing and refinement statistics for the X-ray crystal structures, 2Fo–Fc electron density maps for 15, 16, 19, and 20 when in complex with galectin-3C, electron densities for the alternate positions of 20 compared to 16, 2Fo–Fc electron density maps for galectin-3 in complex and R144 and R186 interactions with 15, 16, 19, and 20, and copies of 1H NMR spectra for compounds 4, 5, 7, 9, 10, 11, 13, 14, 16, 17, 19–22. See DOI: 10.1039/c7sc04749e |
‡ The structures have been deposited in the Protein Data Bank with accession numbers 4BLJ (15), 4BLI (19), 5IUQ (16), and 4BM8 (20). |
This journal is © The Royal Society of Chemistry 2018 |