|
|
||||||||
J. Biol. Chem., Vol. 281, Issue 4, 2317-2332, January 27, 2006
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



1
From the
Unité de RMN des Biomolécules, URA CNRS 2185, Institut Pasteur, **Unité de Chimie Organique, URA CNRS 2128, Institut Pasteur, ¶Unité de Pathogénie Microbienne Moléculaire, Institut Pasteur, 28 Rue du Dr. Roux, 75724 Paris Cedex 15,
DPM UMR5063 UJF/CNRS, 5 Avenue de Verdun 38240 Meylan, France, and the ||CERMAV-CNRS (affiliated with Université Joseph Fourier), 38041 Grenoble BP53, Cedex 09, France
Received for publication, September 15, 2005 , and in revised form, October 25, 2005.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
However, access to the required carbohydrate haptens is often a roadblock. Therefore, besides the investigation of anti-idiotype antibody (7), expanding the concept of mimicry led in the recent past to extensive exploring of the potential mimicking of polysaccharide and/or complex oligosaccharide antigens by peptides (8-10). These peptide mimotopes, i.e. peptide mimics inducing an anti-carbohydrate antibody response upon immunization, have been proposed as potential surrogate antigens of carbohydrates in vaccine development (10). Indeed, because of their ease of manufacture and their intrinsic immunogenic properties, peptide mimotopes may have greater advantage over complex carbohydrate haptens issued from bacterial cell cultures or low yielding multistep syntheses. However, not all peptide mimics of carbohydrate antigens behave as mimotopes. Despite the large number of known peptide mimics, only few peptide mimotope-based experimental vaccines have been reported so far (11-16).
It is believed that a better understanding of the molecular basis of peptide-carbohydrate mimicry could help the rational design of potent peptide mimotope-based vaccines. In particular, whether mimicry is structural, functional, or both remains an unsolved question (9, 17-19). Available x-ray data of carbohydrate-protein and of the corresponding peptide mimotope-protein complexes along with information on the thermodynamics of peptide mimic-protein binding are somewhat scarce (17, 20, 21). Thus to date, although analysis of the topography of ligand-receptor complementarity may be performed by a variety of methods, available knowledge on the molecular features of peptide-carbohydrate mimicry mostly relies on data obtained from NMR and molecular modeling studies as reviewed recently (22).
|
Besides, based on a combination of NMR and molecular modeling studies, we proposed a conformational model for the S. flexneri 5a O-SP whose biological repeating unit is the branched pentasaccharide I (Structure 1) (31). Study of both the antigenicity and the conformation of the four synthetic frame-shifted pentasaccharides corresponding to pentasaccharide I (32) suggested that the DA(E)BC sequence is the structure that best mimics the native O-SP antigen (33). More recently, the pentasaccharide DA(E)BC was shown to act as a mimotope.4
Here we report the antigenicity and the NMR findings on the preferred conformation of p100c and p115 peptide mimotopes both in their free and mIgA-bound forms. Analysis was also performed using peptide p22 (KRHFLSQRQ, mIgA C5- and mIgA I3-specific), one of the 17 nonimmunogenic peptides selected during the original screening (28). Antibody-bound conformations and epitope mapping were derived from transferred NOE (trNOE) (34, 35) and saturation transfer difference (STD) experiments (36), respectively. The conformational preferences observed for the peptides were tentatively related to those derived from NMR and molecular modeling analysis of the DA(E)BC-mIgA complexes that led to a theoretical model of the recognition of S. flexneri 5a O-SP by mIgA I3. This contribution adds to the few reports investigating molecular mimicry by analyzing both peptide mimic-mAb and carbohydrate-mAb recognition features (37-39).
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Inhibition ELISACharacterization of the oligosaccharide determinant recognized by the mIgA was performed by measuring the mIgA-oligosaccharide interaction as follows. First of all, a standard curve was established for each mIgA tested. Different concentrations of the mAb were incubated overnight at 4 °C on microtiter plates coated with purified S. flexneri 5a LPS at a concentration of 5 µg/ml in carbonate buffer, pH 9.6, and subsequently incubated with 1% PBS/BSA for 30 min at 4 °C. After washing with PBS/Tween 20 (0.05%), alkaline phosphatase-conjugated anti-mouse IgA was added at a dilution of 1:5,000 (Sigma) for 1 h at 37 °C. After washing with PBS/Tween 20 (0.05%), the substrate was added (12 mg of p-nitrophenyl phosphate in 1.2 ml of 1 M Tris-HCl buffer, pH 8.8, and 10.8 ml of 5 M NaCl). Once the color developed, the plate was read at 405 nm (Dynatech MR 4000 microplate reader). A standard curve A = f([Ab]) was fitted to the quadratic equation Y = aX2 + bX + c, where Y is the absorbance and X is the Ab concentration. Correlation factor (r2) of 0.99 was routinely obtained.
Then the amount of oligosaccharides giving 50% inhibition of mIgA binding to LPS (IC50) was determined as follows. Each mIgA at a given concentration (chosen as the minimal concentration of Ab which gives the maximal absorbance on the standard curve) was incubated overnight at 4 °C with various concentrations of each of the oligosaccharides to be tested, in 1% PBS/BSA. Measurement of unbound mIgA was performed as described above using microtiter plates coated with purified LPS from S. flexneri 5a, and the mAb concentration was deduced from the standard curve.
The recognition capacity of anti-LPS mIgA for LPS was determined as described above using various concentrations of LPS that were incubated in solution overnight at 4 °C with the predefined concentration of each mIgA. IC50 was defined as the concentration of oligosaccharides required to inhibit 50% of mIgA binding to LPS.
NMR SpectroscopyAll 1H NMR experiments were recorded at 298 K on a Varian Unity Inova spectrometer operating at 1H frequencies of 500 MHz. 1H chemical shifts were given relative to an external standard of 4,4-dimethyl-4-silapentane sodium sulfonate at 0 ppm.
Free PeptidesThe samples were prepared in 90% H2O and 10% D2O at pH 5 for p115 and p100c and at pH 6.5 for p22. The solution concentrations were about 10, 3, and 8 mM for p115, p100c, and p22, respectively. DQF-COSY (42), TOCSY (43), and ROESY (44) experiments were recorded with 512 increments and 16 scans at 298 K. The TOCSY and ROESY experiments were acquired using a mixing time of 80 and 400 ms, respectively. Water suppression was performed using the WATERGATE pulse sequence (45). All NMR spectra were collected in the phase-sensitive mode using the States-Haberkorn method (46).
Ligand-Antibody InteractionsShigemi tubes were used for all samples. In order to prepare NMR samples of pentasaccharide DA(E)BC-OMe in the presence of the antibodies mIgA C5 and mIgA I3, mAbs were concentrated after repeated cycles of exchange with D2O buffer (50 mM deuterated sodium phosphate, 100 mM NaCl, pH 6.5) in Amicon Centriprep-10 concentrators. trNOE experiments (34, 35, 47) performed on different pentasaccharide:binding site ratios (5:1, 10:1, 15:1, 20:1, and 30:1) showed that the most favorable ratio for trNOE was 20:1. So the final samples were prepared with 3.75 µM antibody and 0.3 mM pentasaccharide in 380 µl of the above mentioned D2O buffer. trNOE, trROE, and STD experiments (36) on pentasaccharide DA(E)BC-OMe in the presence of mIgA C5 and mIgA I3 were recorded at 500 and 600 MHz, respectively. trNOE experiments were performed with mixing times of 100, 150, 250, 300, and 400 ms at 303 K to obtain build-up curves and trROE with a mixing time of 400 ms.
The conformation of the free peptides was studied at pH 5. However, with this pH value being close to the isoelectric points of the mIgAs, a study of peptides in their bound conformation was performed at pH 6.5 to avoid precipitation of the mAb. Similarly to the DA(E)BC-mIgA complexes, a peptide:antibody-binding site ratio of 20:1 was used (0.3 mM:3.75 µM). trNOE experiments were performed with mixing times of 100, 150, 250, 300, and 400 ms. To be sure that the observed negative cross-peaks were real trNOEs, NOESY spectra were recorded under the same pH, temperature, and concentration values with the peptides alone. Furthermore, to discard any impact on NOE effects of viscosity increase as a result of the mAb presence, a NOESY spectrum (
m = 200 ms) of p115 was registered in the presence of BSA at the same concentration ratio as that used with the mIgAs. Because no negative NOE cross-peaks were observed in either case, it was assumed that the negative NOEs observed in the presence of mIgA were trNOEs.
Selective saturation of antibody resonances were performed for all STD-NMR experiments at 0.3 ppm (30 ppm for reference spectra) using a series of 40 gaussian-shaped pulses (50- and 10-ms delay between pulses, excitation width
B1/2
, approximately 50 Hz) for a total saturation time of 2.4 s. The one-dimensional STD spectra were recorded with 4096 scans at 288 and 298 K for the pentasaccharide and the peptides, respectively. Subtraction of saturated spectra from reference spectra was obtained by phase cycling (36). For DA(E)BC-OMe, two STD-TOCSY experiments (48) were recorded with selective saturation at 0.3 and 30 ppm, respectively. Differences between the two spectra were performed using the VNMR software. No attempt here was made to quantify STD-NMR intensities, as it is known that these exhibit a complex dependence on relaxation times, correlation times, exchange rates, and on binding site proton density. Indeed, only when short saturation times are used, i.e. less than 1 s, can intensities reflect ligand proton-protein proton distances (37). Here the saturation time of 2.4 s prevented us from quantitative analysis.
Distance and Angle ConstraintsThe cross-peak volumes from trNOESY and trROESY experiments of the pentasaccharide in the presence of mIgAs were measured with the VNMR software. Distances between neighboring protons were calculated by the usual 1/r6 NOE/distance relationship (49). NOE-derived and trROE distances were obtained from initial NOE build-up rates, which were calculated by NOE volumes fitting during different mixing times. The intra-residue distance of 2.52 Å between the H-1 and H-2 protons of the
-L-rhamnopyranosyl unit B was used as a reference for distance calibration.
Distance constraints of free peptides were obtained from the ROESY spectrum run at 298 K with a 400-ms mixing time. For peptides in the presence of mIgA, distance constraints were obtained from the trNOESY spectra run at 298 K with a 200-ms mixing time. NOE intensities were evaluated from the height of the cross-peaks. For structure calculations, upper limit distances of 2.8, 3.5, and 5 Å were used for strong, medium, and weak NOEs, respectively (50). The 3JNH-H
values were used to restrain
angles as follows: for 3J >9 Hz, -155° <
< -85°; for 8 Hz < 3J <9 Hz, -175° <
< -65°; for 5 Hz < 3J <7 Hz, -105° <
< -55°; for 3J <5 Hz, -90° <
< -40° (51).
Structure CalculationsStructure calculations of free and bound peptides were run on a Silicon Graphics work station using the standard protocol of the DYANA program (52). A total of 100 structures were calculated using the torsion angle dynamics protocol. The structures were sorted according to the final value of the target function, and the 20 best structures were analyzed in terms of distance and angle violations. Of these 20 structures, the 10 best structures were visualized by using MOLMOL (53).
Homology Modeling of the IgA I3 Fab Fragment and DockingThe search for structures with sequences similarities was performed with Blast (54) on sequences of all proteins with known three-dimensional structure in the Protein Data Bank (55). Five structures of interest were downloaded and used as template by the Composer program for the building of VL and VH chains of IgA I3 (56).
The Tripos force field (57) option of the Sybyl program (SYBYL) was used to minimize the energy of the resulting model whose stereochemical features were validated with the PROCHECK program (58).
The Autodock3 program (59) was used for docking oligosaccharides and peptides in the binding site of modeled IgA I3 Fab. Because the goal was to model the behavior of the O-SP, calculations were performed on the largest possible fragment compatible with the limitations of the software, in that case a nonasaccharide. The 9-carbohydrate residue fragment was thus chosen as BCDA(E)BCDA in which the key pentasaccharide DA(E)BC is flanked by two residues on each side. The two conformations that were shown previously to correspond to helical shapes of the O-SP (33) were used as starting models. Hydroxyl and N-acetyl bonds were considered as flexible, whereas glycosidic bonds were considered as rigid to keep the helical conformation, resulting in 28 degrees of freedom. AMBER force field charges were assigned to all protein atoms, and partial charges were assigned to the atoms according to the PIM force field (60). Grids of probe atom interaction energies and electrostatic potential were generated around the whole protein by the AutoGrid program present in Autodock3 with a spacing of 0.5 Å. All probes were placed arbitrarily at a distance of 10 Å from the protein surface, and their exocyclic torsion angles were allowed to rotate freely. For each monosaccharide, one job of 240 docking runs was performed using a population of 100 individuals and an energy evaluation number of 10 x 106. Clustering of solutions was done by root mean square fitting (<1 Å). The best solution of each cluster was used to propagate the helices to 20 residues while keeping the conformations determined previously (33). Twenty different conformers of the p100c peptide were also docked in the mIgA I3 Fab-binding site using the rigid body approach of the Autodock3 program. For each of them, one docking run was performed using a population of 100 individuals and an energy evaluation number of 0.75 x 106.
| RESULTS |
|---|
|
|
|---|
IC50 values for recognition of the mimotopes by mIgAs revealed that p100c was better recognized by mIgA I3 (IC50 = 75 ± 29 µM) than by mIgA C5 (IC50 >1000 µM). In contrast, p115 was better recognized by mIgA C5 (IC50 = 197 ± 39 µM) than by mIgA I3 (IC50 >1000 µM). Most interestingly, p22 exhibited a higher IC50 value for both mIgAs (70 ± 11 µM and 0.03 ± 0.01 µM for mIgA I3 and mIgA C5, respectively) than those measured for p100c and p115.
|
coupling constants, and proton-proton dipolar interactions observed in the ROESY spectra. Dihedral angles and distance constraints deduced from these data were then used to model the averaged solution structure of each peptide analyzed with the DYANA program (52).
Peptide 115 (KVPPWARTA)The NMR spectrum revealed that p115 displayed four different conformers as a result of the cis-trans isomerization of the amide bonds involving the two sequential prolines, Val2-Pro3 and Pro3-Pro4, respectively. Based on signal intensities, it was estimated that the major conformer represented 80% of the different species, whereas the three other forms altogether made up for the remaining 20%. Because no information was available for the conformer recognized upon selection from the phage displayed peptide library, structural analysis was conducted for this major conformer only. Chemical shifts and three bond 3JNH-H
coupling constants are reported in Table 1. Significant deviations from random coil values are only observed for the H-
protons of Pro3 and Pro4, whereas all three bond 3JNH-H
coupling constants are those expected for flexible peptides. Inter-residue dipolar interactions observed in the ROESY experiment between the H-
of the residue preceding a proline and the H-
proton of the proline indicates that both Pro3 and Pro4 adopt a trans-conformation in the peptide major conformer. In addition to standard sequential interactions, several medium range interactions were also observed between side chain protons of Val2, Pro3, and Pro4 and the CH3-
protons of Ala6 (Fig. 1). These ROE connectivities were used as distance constraints to model the conformation of the p115 major conformer using the DYANA program. The 10 best structures, i.e. with the lowest energy function, showed that p115 adopts a rather organized conformation comprising residues Pro3 to Arg7, whereas the N- and C-terminal ends are quite flexible (Fig. 2) as expected for peptides of this size. Based on C-
i-C-
i+3 distances as well as on
and
angle values (Table 2), the conformation of the Pro3-Arg7 fragment can be described as two sequential
-turns, a nonclassified one for PPWA and a type I
-turn for PWAR (61, 62).
|
|
protons of residues Tyr2, Lys3, Pro4, Leu5, and Gly6 suggesting some restricted flexibility along the Tyr2-Gly6 sequence. Furthermore, the three bond 3JNH-H-
coupling constant values for residues Leu5, Gly6, and Ala7 are slightly smaller than those measured for the other residues, 5 Hz versus 7-8 Hz (Table 3), strengthening the hypothesis of a probable structuring of the Lys3-Ala7 segment. Inter-residues dipolar interactions observed in the ROESY experiment between H-
of Lys3 and H-
of Pro4 indicate that Pro4 adopts a trans-conformation in the major conformer of p100c. In addition to standard sequential interactions, several medium range interactions were also observed between side chain protons of residues Pro4 to Ala7 as for example between all protons of Pro4 and the methyl group of Ala7 (Fig. 1). Moreover, four long range ROE connectivities were observed between Tyr2 and His10 protons, confirming the cyclic nature of the peptide (Fig. 1). ROE derived distances and coupling constants were used as constraints to generate a family of structures for p100c using DYANA. The 10 best structures indicate that the Pro4-Pro7 fragment of p100c is conformationally organized into a type I
-turn (Fig. 2 and Table 2) (61, 62). Because both the cyclic form and Pro4 can induce this type of conformational behavior, the structural analysis was extended to reduced p100c (data not shown). The type I
-turn remained, suggesting that Pro4 alone is responsible for its formation, although the cyclic structure might contribute to its stabilization.
|
|
proton, chemical shifts do not significantly deviate from standard values. Meanwhile, none of the coupling constants of internal residues could be measured because of extensive signals overlaps. Nevertheless, medium range ROE connectivities were observed between side chains protons of residues His3 and Leu5 as well as between those of Phe4 and Ser6 (Fig. 1). The 10 structures of lowest energy matching those distance constraints show that fragment His3 to Ser6 of p22 is organized into a nonclassified
-turn (Table 2), although the peptide N- and C-terminal ends remain quite flexible (Fig. 2). That available chemical shift and coupling constant values do not reflect such an organized conformation in solution suggests a weaker stability of the
-turn.
|
Binding of DA(E)BC-OMe to mIgA I3 and mIgA C5Bound pentasaccharide conformation. Independently of the mIgA tested, the best DA(E)BC-OMe:mIgA ratio to observe trNOEs was shown to be 20:1 in binding sites. Because the highest attainable mIgA concentration was 3.75 µM, a 0.3 mM concentration of DA(E)BC-OMe was used to fulfill the 20:1 ratio requirement. Because NOE intensities depend, among other parameters, on the correlation time for reorientation and therefore on temperature, the later parameter was optimized so that 
c equals 1 for the free pentasaccharide, thus allowing us to distinguish NOE from trNOE connectivities. Indeed, the NOESY spectrum of the free pentasaccharide at 30 °C displayed only few positive and weak NOE connectivities characteristic of small molecules, whereas negative NOE connectivities were observed in the trNOESY spectrum of DA(E)BC-OMe when interacting with either mIgA C5 or mIgA I3. Because these effects did not result from the increased solution viscosity as probed by recording NOESY spectra of the pentasaccharide in the presence of BSA at the same w/v concentration (37), it was assumed that they corresponded to trNOE connectivities for the mIgA-bound DA(E)BC-OMe.
trNOESY spectra obtained with several mixing times, ranging from 100 to 400 ms, allowed us to trace the build-up curves (trNOEs intensities versus
m) from which the distance information was extracted. In addition, inter-residue 1H-1H distances were also calculated from a trROESY spectrum obtained with a mixing time of 400 ms to take spin diffusion into account, if any. Comparison of these distances with those measured for unbound DA(E)BC-OMe (33) suggested that the pentasaccharide conformation was not significantly modified upon binding to the mIgAs (Table 5).
|
Binding of the Peptide Mimics to mIgAsBased on available IC50 values, analysis was run on the p115-mIgA C5, p100c-mIgA I3, p22-mIgA C5, and p22-mIgA I3 complexes.
Interaction of Peptide 115 (KVPPWARTA) with mIgA C5trNOE experiments recorded for the p115-mIgA C5 complex (see supplemental Fig. 4S) showed new NOE connectivities when compared with those observed for the free peptide. These additional cross-peaks, such as those observed between residues Trp5 and Ala6, were clearly identified as representative of the p115-bound form. Interestingly, most of the NOE connectivities involving amide protons of the free p115 were no longer observed in bound p115 with the exception of the Trp5, Ala6 amide proton connectivity (see supplemental Fig. 1S). The pH increase from 5 in the free peptide to 6.5 in the peptide:mIgA solution might account for such experimental observations since amide protons exchange faster at higher pH. Nevertheless, the medium range NOE connectivities observed between the side chain proton of residues Val2 and Pro3, and the CH3-
of Ala6 remained (see supplemental Fig. 1S). Distance constraints derived from trNOE intensities were used to establish the conformation of p115 when bound to mIgA C5. Superimposition of the 10 lowest energy backbone conformations of free p115 to those of mAb-bound p115 showed that only the turn involving residues Pro3 to Ala6, observed in the free form, is maintained in the bound form. Based on C-
i-C-
i+3 distances as well as on
and
angle values (Table 2) the type I
-turn observed between residues Pro4 and Arg7 for the free peptide is no longer present. Whereas in the free form fragment Pro3-Ala6 adopts a nonclassified type of
-turn, in the bound form a well defined type II
-turn is clearly present (61, 62) (Fig. 2 and Table 2).
|
-turn is crucial for p115:mIgA C5 recognition. It is worth noting that the major form of p115 was that recognized by mIgA C5.
Interaction of Peptide p100c (CYKPLGALTHC) with mIgA I3As compared with data corresponding to the free form, the trNOESY spectrum of p100c in interaction with mIgA I3 displayed new data specific for the bound form. These include new sequential NOEs such as those observed between Lys3 and Pro4 or between Pro4 and Leu5, whereas sequential NOEs between Ala7 and Leu8, or between Leu8 and Thr9, were no longer visible. However, medium range NOE connectivities observed between Pro4 and Ala7 in the free form remained. In addition, medium range interactions specific to the bound form were observed, such as those involving side chain protons of Tyr2 and Leu5, Lys3 and Leu8, as well as Tyr2 and His10 (see supplemental Fig. 2S). Superimposition of the backbone (Pro4 to Ala7 fragment) of the 10 lowest energy conformations of free and mIgA I3-bound p100c matching the distance constraints showed that, as observed for p115, the turn observed in the free form was maintained in the bound form. Furthermore, data pointed to a switch from a type I
-turn in the free form to a type II
-turn in the bound form (Table 2) (61, 62). Additional significant rearrangements were observed for the rest of the backbone (Fig. 2).
|
|
protons, were in contact with the mIgA I3 (Fig. 7). Within the turn, only Leu5 and Ala7 methyl groups contacted the mAb. Similarly to p115, the methyl groups of the hydrophobic residues and the aromatic residue of p100c were involved in the interaction with the mIgA I3.
|
-turn involving the His3-Ser6 segment in the free form changed to a type II
-turn in the bound form, as observed for p115 and p100c. Additional rearrangements were observed for the rest of the backbone (Fig. 2). trNOE experiments for p22 bound to the mIgA C5 were unsuccessful because sparse negative NOEs were observed (data not shown). The high affinity of mIgA C5 (IC50 of 0.03 µM), most probably associated with an equilibrium constant for dissociation (Kd) below 10-6 M and thus not compatible with trNOE observations (10-3 > Kd > 10-6 M), might be responsible for this effect. Despite this drawback, epitope mapping by STD experiments was successfully undertaken as the lower limit for exchange was less stringent in terms of Kd (10-8 M). Epitope identification for p22 bound to mIgA I3 and mIgA C5, respectively, was obtained from the one-dimensional STD experiments (Fig. 8). Whether mIgA I3 or IgA C5 was concerned, p22 residues in direct contact with the mAbs were identical. They included the His3 imidazole protons, Phe4 aromatic protons, as well as Leu5 methyl protons, all corresponding to amino acids taking part in the turn observed in free p22. Although the conformation of mIgA C5-bound p22 remained undisclosed, available data suggest that the nonclassified turn naturally adopted by p22 contributed to both mIgA I3 and mIgA C5 recognition.
|
|
-D-glucopyranose residue E is deeply buried into the central pocket of the groove. In all of the four solutions, the mAb features involved in carbohydrate recognition are the three loops from the heavy chain and the L1 and L3 loops from the light chain. The variable loop H3 plays the major role, with its two Asp residues (Asp91 and Asp92) involved in recognition for most binding modes. When the nonasaccharide-bound conformations were propagated into polysaccharide structures comprising four repeating units, only two docking modes appeared to be stable with additional contacts created on both sides of the binding site. In each case, this ability corresponded to the parallel arrangement of BCDA(E)BCDA. Therefore, only these two solutions, i.e. docking of E and O conformations in parallel mode, were considered as possible mimics of O-SP binding to mIgA I3. The two possible docking modes of BCDA(E)BCDA and the polysaccharide of DP4 are displayed in Fig. 10 (A-D) and the contacts of interest are listed in Table 6. For both conformations of each ligand, the central trisaccharide A(E)B makes most of the binding contribution, and the additional contacts established by the GlcNAc (D) residue are minor.
|
|
| DISCUSSION |
|---|
|
|
|---|
Here, by aiming at designing new vaccine strategies against Shigella infection, we developed synthetic mimics, carbohydrates and peptides, of S. flexneri 5a O-SP. Previous NMR and molecular modeling studies from our laboratory have shown that, among the four possible frame-shifted pentasaccharides representative of the O-SP, DA(E)BC-OMe best mimics the conformational features of S. flexneri 5a O-SP (33). More importantly, the trNOE data reported here indicate that the conformation of DA(E)BC-OMe when bound to O-SP-specific mIgA I3 and mIg C5 does not differ from its conformation when free in solution. Noteworthy, selection of free solution conformers often predominates in mAb-carbohydrate recognition processes (70). However, it is not always so as exemplified with mAb Se155-4 binding to a trisaccharide antigenic determinant of the Salmonella paratyphi B branched O-SP (71).
Another interesting example of such induced conformational change is the mAb SYA/J6 binding to the pentasaccharide ABCDA' fragment of the linear O-SP defining S. flexneri serotype Y (72). Modeling of the linear heteropolysaccharide has shown that it is structured into a left-handled helical chain of three ABCD repeating units (73, 74). Most interestingly, extension of the modeling study to S. flexneri 5a O-SP suggested that residue E, which is associated with serotype specificity, crucially impacts the overall O-SP conformation. In fact, the branched heteropolysaccharide, whose repeating unit (I) bears the E side chains, behaves as a right-handed 3-fold helix with residue E protruding outwardly (33). Along this line, the interaction of DA(E)BC-OMe with the serotype-specific mIgAs is mainly driven by the branched E residue as evidenced by the large number of E-specific signals enhanced in the STD spectra of the pentasaccharide in complex with mIgA I3. Furthermore, NMR experiments showed that all rhamnose methyl groups are also in close contact with the mAb, with the methyl group of rhamnose B on which E is branched being the major contributor. The N-acetyl group of residue D gives only weak contacts with the mIgA protons, indicating that it probably lies at the surface of the mAbs. These data are supported by the inhibition ELISA results showing that all frame-shifted tri-, tetra-, and pentasaccharides, bearing A(E)B-branched trisaccharide characteristic of S. flexneri 5a serotype, are recognized by a protective serotype 5a-specific mIgG (33). Further insights on the central role played by the branched
-D-glucopyranosyl residue E in mAb recognition derives from docking of the nonasaccharide BCDA(E)BCDA and fragments of the O-SP comprising four repeating units in the mIgA I3 Fab-binding site. The latter appears to have a distinct groove character with a deep central pocket, a type of binding site often encountered for mAbs binding internal polysaccharide sequences (75). Most interestingly, this finding is identical to that observed for Strep 9, a mouse mAb of the IgG3 subclass directed against the cell wall polysaccharide of group A Streptococcus (76) made of repeats comprising a branched
-N-acetyl-D-glucosamine residue linked to a linear di-rhamnopyranosyl backbone. As found earlier, sides of the mIgA I3 groove are flanked by CDR2 of the heavy chain and CDR1 of the light chain. Moreover, aromatic residues such as Tyr45 of the heavy chain and Tyr34 of the light chain define the pocket region, pointing once more to the importance of such amino acids in carbohydrate recognition (77, 78). Indeed, whatever the orientation of nonasaccharide BCDA(E)BCDA, parallel or perpendicular relative to the groove binding site, the glucose residue E was always deeply buried in the central pocket of the groove and was poorly accessible to solvent. Most interestingly, relying on molecular modeling only, a heptasaccharide related to Brucella abortus O-SP exemplifies another O-SP-mAb interaction for which the mAb-binding site identified as a groove bearing a pocket in its center could also accommodate two binding modes of an O-SP fragment (79). As shown here, docking of S. flexneri 5a O-SP large fragments in the mIgA I3-binding site pointed to only one possible binding mode of the O-SP, namely the parallel mode, independently of the length of the helical repeat taken into account. Thus, in addition to the branched glucopyranosyl residue E behaving as an anchor and to the trisaccharide A(E)B providing the critical epitope exposed on the O-SP, chain elongation also takes part in O-SP:mIgA recognition, highlighting the essential contribution of some kind of conformational epitope or presentation in an extended surface. However, we are aware that small changes at the VL:VH interface of the mAb may result in significant alteration of the binding mode, which cannot be predicted at this stage (80). Thus, data provided here are only meant to provide a model of S. flexneri 5a O-SP binding to a homologous protective mIgA, which needs to be further assessed based on crystallographic data.
|
30 nM), our model fits to others, such as that on Cryptococcus neoformans (18, 84) and that on N. meningococcus C (30), suggesting that commonly used parameters for selecting peptide mimicking polysaccharide antigens, such as high-affinity binding to mAb, are not predictive of the ability of the selected peptides to act as mimotopes. Indeed, based on x-ray analysis, several lines of evidence support the idea that peptide binding to an anti-polysaccharide mAb may differ significantly from that of the natural antigen. On one hand, data on the dodecapeptide PA1 mimicking C. neoformans CPS suggest poor steric complementarity between PA1 and the heavy chain of the mAb used for selection, which may explain why PA1 acts only as a partial mimotope (20). On the other hand, an octapeptide functional mimic of S. flexneri serotype Y O-SP was found to complement the shape of the groove-type binding site of mAb SYA/J6, used for selection, much better than the ABCDA' pentasaccharide fragment of the O-SP (21). However, the peptide does not fully complement the deep pocket located in the center of the groove and occupied by rhamnose C upon binding of ABCDA'. This may explain, at least in part, the poor ability of the octapeptide to behave as an immunogenic mimic (9).
Not surprisingly, although the three peptides do not share any consensus sequence, the conformations they adopt in their free form encompass many rapidly interconverting conformers with short internal sequences spending long lifetimes organized in turn like motifs. Although p22 was found much more flexible than p115 or p100c, all three peptides adopted
-turn conformations, either of nonclassified type or of type I. This appears to be a rather common conformational feature for short peptides representative of antigenic regions of proteins (85) or polysaccharide antigens such as group A Streptococcus CP (9) and group B Streptococcus CP (38). Indeed, it has been suggested that a
-turn allows appropriate exposure of side chain residues for optimal fit within the mAb combining site. Most interestingly, in the later example, peptide FDTGAFDPDWPA, a molecular mimotope of the CPS, was earlier thought to adopt a nonrandom coil conformation in aqueous solution assimilated to a nascent helix that could potentially mimic the extended helical form of the natural carbohydrate epitope (29). This discrepancy underscores the high complexity of conformational studies dealing with short peptides. The relative heterogeneity of
-turn types adopted by p115, p100c, and p22 in the free form is completely lost in the bound form, as the three peptides adopt a type II
-turn conformation, which appears to be crucial for binding independently of the involved mAb. Thus, the lack of consensus sequence among the selected peptides seems to be compensated by structural consensus induced upon fitting to the mAb combining sites. Moreover, the type II
-turn structure starts from a proline (Pro3 and Pro4, respectively) and ends with an alanine (Ala6 and Ala7, respectively) for both p115 and p100c, underscoring the partial structural resemblance between the two peptides. Major contributions of the p115-turn to binding involve aromatic Trp5 and cyclic Pro4, whereas hydrophobic Leu5 was the residue most involved in mAb binding to the p100c-turn. p22 differs notably from p115 and p100c because it has no proline. In this case, aromatic His3 and Phe4 together with hydrophobic Leu5 are the major turn components contributing to binding independently of the mAb. Noteworthy, additional residues do not seem to be engaged in mAb recognition, which may explain the ability of p22 to bind the two mIgAs. On the contrary, going from His1 to Ala9, p115 binding to mIgA C5 necessitates that most residues along the peptide chain, especially those at the N terminus, make specific contacts with the combining site. Similarly, residues at the N terminus of p100c appear critical for peptide binding to mIgA I3. In particular, the two docking models obtained for p100c interacting with mIgA I3 reveal a salt bridge formed between the peptide Lys3 residue and the residue Asp92 within the CDR H3 loop, similarly to those observed between rhamnose A or glucose E and Asp92 or Asp91 of the mIgA CDR H3 loop, respectively, upon DA(E)BC binding. Although not probed at this stage, analogous ionic contributions to peptide-mAb interactions may be anticipated because all selected peptides share basic residues. However, available data suggest that independently of the peptide mimic under study, all mIgA-peptide interactions derived mostly from the direct contact of peptide aromatic residues and methyl groups with the mAb-binding site, suggesting that recognition was basically driven by hydrophobic and van der Waals contacts. In that matter, data provided for the S. flexneri 5a system fully support previous observations made for other models implicating peptides mimicking polysaccharide antigens in complex with specific mAbs (11).
In addition, all data reported here strongly emphasize the crucial role played by the type II
-turn topology in governing molecular mimicry of S. flexneri 5a O-antigen. Most interestingly, superimposition of family of conformations obtained for bound p22 with this obtained for bound DA(E)BC-OMe shows that the type II
-turn in the peptide seems to mimic the nascent helicoidal shape of the oligosaccharide main chain with the aromatic ring of Phe4 having the same orientation as the crucial branched glucose E (Fig. 11). This is in agreement with previous observations showing that aromatic amino acids are considered as ideal residues for mimicking glycan side chain structures (69, 86) and that
-turn/extended structures may be accurate conformational mimics of helices (87, 88). Yet, except for the number of residues involved, discriminating between the binding modes of the three peptides remains difficult. Thus, based on binding complementarity only, rules governing the sel