ORIGINAL RESEARCH article
Partial Dissociation of Truncated Peptides Influences the Structural Dynamics of the MHCI Binding Groove
- Center for Theoretical Chemistry, Faculty of Chemistry and Biochemistry, Ruhr-University, Bochum, Germany
Antigen processing on MHCI involves the exchange of low-affinity peptides by high-affinity, immunodominant ones. This peptide editing process is mediated by tapasin and ERAAP at the peptide C- and N-terminus, respectively. Since tapasin does not contact the peptide directly, a sensing mechanism involving conformational changes likely allows tapasin to distinguish antigen-loaded MHCI molecules from those occupied by weakly bound, non-specific peptides. To understand this mechanism at the atomic level, we performed molecular dynamics simulations of MHCI allele B*44:02 loaded with peptides truncated or modified at the C- or N-terminus. We show that the deletion of peptide anchor residues leads to reversible, partial dissociation of the peptide from MHCI on the microsecond timescale. Fluctuations in the MHCI α2−1 helix segment, bordering the binding groove and cradled by tapasin in the PLC, are influenced by the peptide C-terminus occupying the nearby F-pocket. Simulations of tapasin complexed with MHCI bound to a low-affinity peptide show that tapasin widens the MHCI binding groove near the peptide C-terminus and weakens the attractive forces between MHCI and the peptide. Our simulations thus provide a detailed, spatially resolved picture of MHCI plasticity, revealing how peptide loading status can affect key structural regions contacting tapasin.
To perform their signaling function at the cell surface, major histocompatibility complex class I (MHCI, Figure 1B) molecules (1–4) are first loaded with antigen peptides in the lumen of the endoplasmic reticulum (ER) by the peptide-loading complex (PLC, Figure 1A) (5), a large multi-protein assembly whose structural organization remains to be resolved at the atomic level. A central component of the PLC, tapasin (6, 7), edits the antigen repertoire exposed at the cell surface (8–10) by selecting peptides according to their C-terminus (11). To do so, tapasin accelerates the off-rate of MHCI-bound low-affinity cargo (12–14). Antigen candidates are thus rapidly exchanged until a high-affinity, immunodominant one binds (14, 15). Concomitantly, the aminopeptidase associated with antigen processing in the ER (ERAAP) (16–18) also edits the antigen repertoire, by cleaving peptides at their N-terminus (11) and thus adjusting their length so that they fit optimally into the MHCI binding groove (19).
Figure 1. The peptide-loading complex (PLC) and selected components. (A) The PLC is assembled in the ER. Tapasin, a chaperone and peptide exchange catalyst, bridges the TAP transporter to MHCI. Accessory proteins ERp57 and calreticulin (Crt) stabilize the complex. ERAAP refines peptides through its peptidase activity. (B) MHCI consists of a heavy, variable α chain and the invariant, light β2m chain. The α1 and α2 domains form the peptide-binding groove, while α3 bears the CD8 recognition loop. Both the α2−1 helical segment and the CD8 loop contact tapasin in the PLC. (C) The antigen is bound to the MHCI groove in an extended conformation; the two groove regions interacting with the peptide N- and C-terminus are termed the A- and F-pockets, respectively. Alternating colors distinguish sequential residues.
According to the predicted structure of the tapasin–MHCI complex (14, 15), direct contact between tapasin and the MHCI-bound peptide is very unlikely. Indirect sensing and loading mechanisms have, therefore, been proposed to explain how tapasin influences peptide binding to MHCI. In a previous work (15), we have shown that the tapasin N-terminal domain (TN) acts on MHCI by pulling on its α2−1 region near the peptide C-terminus, thereby widening the binding groove and facilitating the release of low-affinity peptides. We showed that peptides compete with tapasin by pulling the same α2−1 region in the opposite direction, thus tightening the binding groove. Forces exerted by high-affinity peptides overcome those of tapasin, closing the groove and priming the PLC for dissociation (15). This is consistent with the observation that an engineered MHCI variant (KbY84C) whose α2−1 region is linked to α1 by a disulfide bond, which can exert strong closing forces, is able to breach cellular quality control (20). However, other investigations have highlighted the importance of the MHCI α3 domain for peptide editing (21, 22), in line with α3 contacting tapasin at its C-terminal domain (TC). Allosteric coupling between α3- and the α1α2 binding groove domain has been suggested for chicken allele BF2*15:01 (23), and recently also been proposed for the human allele B*44:02 (24), prompting a mechanism according to which tapasin modulates binding groove plasticity through the TC–α3 interface. The direct and allosteric mechanisms are not mutually exclusive, though. They could be involved to different degrees, e.g., in different alleles, in the peptide editing function itself, or in the sensing of properly loaded MHCI.
Several molecular dynamics (MD) simulation studies have probed MHCI dynamics (25–35), how they correlate with peptide cargo (36–38), and their implications for peptide editing by tapasin (14, 15). These studies often compare peptide-loaded (MHCIPL) and -deficient (MHCIPD) molecules (14, 15, 25, 33), because crystal structures of MHCI complexed with low-affinity peptides are usually not available. MHCIPD computational models can be easily prepared from crystal structures of high-affinity antigen-bound MHCI by peptide removal. This approach, however, might suffer from limitations. The PD state may not be particularly relevant since, in vivo, MHCI molecules are loaded with a peptide cargo before they are even recruited by tapasin (39). In the context of the PLC, tapasin stabilizes MHCI in a peptide-receptive state, hence accelerating peptide binding (9). Furthermore, peptide release is known to be the rate-limiting step in peptide editing (14). It follows that MHCIPD is only a transient state, which rapidly leads to the binding of a new peptide. Hence, the essential role of tapasin is to distinguish between MHCI-bound peptides of different affinities, not between the PD and PL forms. MHCIPD may, therefore, not be an optimal proxy for MHCI bound to low-affinity cargo. In that respect, a recent study involving truncated, low-affinity peptides is of particular interest (37). In that work, the authors used both MD simulations and experimental techniques to study the importance of termini residues in two peptides that bind to murine MHCI H-2Kb. From 50 ns MD simulations, they reported increased binding groove RMSD and RMSF for the peptides truncated at either the N- or C-terminus. Thermal denaturation experiments carried out with the same truncated peptides yielded consistent results: the melting point (Tm) of MHCI loaded with truncated peptides was significantly lowered. The peptide C-terminus was also shown to be involved in retaining MHCI at the cell surface.
In the present work, to better understand MHCI structure and dynamics in the context of peptide editing, we have performed microsecond-timescale MD simulations of MHCI allele B*44:02 bound to N- or C-truncated DPA (high-affinity antigen peptide derived from HLA-DPα; Figure 1C). This allows us to correlate MHCI plasticity, including local features of the binding groove, such as the structural stability of the α2−1 region, to specific changes in peptide cargo. We show that the peptide termini anchor it to the groove; removal of a single residue at the C-terminus suffices to trigger partial dissociation of the peptide, while removing two residues at the N-terminus has a similar effect. At the N-terminus, both the peptide main chain and the residue side chains affect stability. Conversely, the side chain of the terminal residue is critical for peptide stability at the C-terminus. Even in the full-length antigen, central residues are somewhat flexible and undergo thermal fluctuations, thus deviating from their position in the X-ray crystal structure over the course of the simulation. Any C-terminal deletion results in a widening of the binding groove near the F-pocket, similar to that observed in MHCIPD. Conversely, N-terminal deletions and peptide removal narrow the binding groove near the A-pocket. The overall structural stability of the groove is heavily influenced by the peptide cargo, with α2−1 being the most affected region, confirming the importance of tapasin TN as an MHCI chaperone. Long-range effects on MHCI plasticity could only be observed in one of the systems under study: B*44:02 loaded with an antigen truncated by a single residue at the N-terminus. In this specific case, the increased fluctuations in the α1α2 domain lead to increased fluctuations, and hence increased configurational entropy, in the α3 domain that contacts tapasin TC via the CD8 recognition loop. Using a C-truncated peptide as proxy for a low-affinity cargo, we performed additional simulations of MHCI in complex with tapasin. We show that tapasin modulates MHCI through an opening of the binding groove that reduces the attractive forces between MHCI and the peptide.
2. Materials and Methods
2.1. System Generation
Coordinates of B*44:02 loaded with the HLA DPA*0201 9-mer peptide were taken from a 1.6-Å resolution X-ray crystal structure (PDB ID 1M6O) (40). The α subunit contains 276 residues and excludes the membrane-spanning helix and cytosolic tail, while the invariant β2m subunit contains 99 residues. Water molecules resolved in the crystal structure were retained. PROPKA 3.1 (41) was used to determine protonation states at pH 7. D156 in the α2−2 helix segment was the only residue with a non-standard protonation state, sharing a proton with D114 in β6 at the bottom of the binding groove.
The structure was first subjected to 60 steps of steepest-descent (SD) energy minimization. Peptide removal, truncation, and side chain replacement were then used to generate from MHCIPL the starting structures for MHCIPD, MHCI loaded with six truncated DPA peptides (Δ1N, Δ2N, Δ3N, Δ1C, Δ2C, Δ3C) and MHCI loaded with three peptides in which one or two residues were replaced by glycines (δ1N, δ2N, δ1C); all peptides were capped with standard (charged) amine and carboxyl termini. All systems were solvated in periodic rhombic dodecahedron cells, with a 10-Å minimal distance between the solute and the box edge. Random water molecules were replaced by Na+ and Cl− ions at a 0.15 M concentration, yielding a neutral net charge. The systems contained ≈69,000 atoms in cells of ≈280,000 Å3. Finally, all systems were again subjected to 500 steps of SD energy minimization.
The tapasin–MHCIΔ1C system was built from the results of a previous simulation of tapasin–MHCI (14, 15) in a solvated periodic rhombic dodecahedron cell with a 10-Å minimal distance between the solute and the box edge, and a 0.15 M concentration of Na+ and Cl− ions. The system contained ≈140,000 atoms in a ≈1,480,000 Å3 cell. The DPA peptide in this system was truncated and capped with a carboxyl terminus prior to 500 steps of SD energy minimization.
2.2. MD Simulations
Simulations were carried out with GROMACS 2016.1 (42, 43). The Amber99SB*-ILDNP protein forcefield (44–49) and TIP3P water model (50) were used. The SETTLE (51) and LINCS (52) algorithms were applied to constrain the internal degrees of freedom of water molecules and the bonds in other molecules, respectively. In combination with virtual site hydrogens (53), this allowed for a 4-fs integration time step. Short-range non-bonded Coulomb and Lennard-Jones 6–12 interactions were treated with a Verlet-buffered pair list (54) with potentials smoothly shifted to zero at a 10-Å cutoff. Long-range Coulomb interactions were treated with the PME method (55) with a grid spacing of 1.2 Å and cubic spline interpolation. Analytical dispersion corrections were applied for energy and pressure to compensate for the truncation of the Lennard-Jones interactions. The thermodynamic ensemble was nPT. Temperature was kept constant at 300 K by a velocity-rescaling thermostat with a stochastic term (56), with coupling time constant 0.1 ps. For constant 1.0-bar pressure, an isotropic Berendsen barostat (57) was used with coupling time constant 0.5 ps and 4.5 × 10−5/bar compressibility.
For each system, five independent 1.0-µs trajectories were acquired (by generating random initial velocities from a Boltzmann distribution), with coordinates recorded every 10 ps. All of these simulations were preceded by 10 ns of equilibration, during which all heavy (non-hydrogen) protein atoms were position-restrained by harmonic potential energy functions with force constants of 1,000 kJ/(mol nm2).
The first 10 ns of simulation without position restraints were considered further equilibration time and were not used for analysis of averaged RMSD, binding groove width, RMSF, configurational entropy, and pairwise forces between residues. In addition, one trajectory involving the Δ1C peptide was excluded from all analyses since the peptide dissociated entirely from the binding groove (see Results for further details).
Configurational entropies were calculated using the quasi-harmonic approximation (QHA) as formulated by Schlitter (58),
which provides an upper bound to the configurational entropy. In this equation, kB is the Boltzmann constant, T the temperature, e Euler’s number, the reduced Planck constant, and M the 3N-dimensional diagonal mass matrix of the N particles. The matrix C is the covariance matrix of particle positions,
where the 3N-dimensional vector x represents the Cartesian coordinates of the N particles for which the entropy is calculated after removing overall translation and rotation by fitting to a reference structure. The starting structure of our simulations was used as the reference structure for this fit. The coordinates of the Cα atoms were used to construct C. By comparing the entropy accumulated for the full 5.0 µs of sampling for each system to the value obtained when discarding a random 10% of the data, we estimated the error and verified that our sampling enables us to evaluate configurational entropies with sufficient statistical reliability.
Force distribution analysis (FDA) was performed using the GROMACS-PF2 implementation (59). Pairwise atom–atom forces were recomputed from the trajectories and were time- and residue-averaged. For the latter, the norm of the sum of the interatomic force vectors between all atoms i and j of residues v and w were calculated, , with i ∈ v and j ∈ w. Attractive and repulsive forces were distinguished by sign.
Starting from an X-ray crystal structure of MHCI allelle B*44:02 loaded with the specific, high-affinity HLA DPA*0201 peptide (40), we set up 8 simulation systems with different peptides. Our reference system is MHCIPL, containing the full-length 9-mer DPA peptide. To simulate lower-affinity peptides, we built 6 truncated systems; these contain the DPA peptide with one, two or three residues deleted at either the N- or the C-terminus, denoted by MHCIΔ1N … MHCIΔ3C. An MHCIPD system was also built by complete peptide removal. Trajectories totaling 5 µs of MD simulation (5 × 1 µs) were acquired for each system.
3.1. Partial Peptide Dissociation on the Microsecond Timescale
Figure 2 shows the Cα RMSD to the X-ray crystal structure for all studied peptides. Deletion of a single C-terminal residue (Δ1C) is sufficient to trigger partial peptide dissociation from the binding groove on the simulation timescale. The C-terminal F9 residue, whose side chain is buried in the F-pocket, is thus crucial for binding to B*44:02. In its absence, central residues cannot stabilize the peptide, and a large segment (G4–S8) dissociates from the MHCI binding groove (Figure 2, Δ1C). Further deletions (Δ2C, Δ3C) have no additional effect; the same G4–S8 region unbinds from the groove. These results are consistent with the observation (36) that removing or modifying the side chain of the peptide C-terminal residue accelerates dipeptide-catalyzed peptide exchange.
Figure 2. Partial dissociation of truncated peptides from B*44:02. Cα RMSD from the initial position in the X-ray crystal structure were averaged for each peptide residue over the five independent 1-µs trajectories acquired for each system. The dotted line represents a chosen dissociation threshold at 3 Å.
In contrast to the C-terminus, a single N-terminal deletion (Δ1N) has a less pronounced effect and leads only to the partial dissociation of residues E2–F3 (Figure 2). Removal of E2, which is buried in the A-pocket, is necessary for dissociation of the F3–R5 region, comparable to the loss of interactions observed for C-terminal deletions. Thus, the peptide is held in place by its termini at E2 and F9, while the residues in between contribute less to MHCI binding. The importance of peptide termini was shown previously from the thermal stability of MHCI molecules loaded with modified peptides (60). Our results are also consistent with the conclusions of a previous MD study of MHCI allele H2-Kb loaded with truncated peptides (37), although the shorter simulation times of that study precluded the observation of partial peptide dissociation.
Due to the thermal fluctuations, even the full-length DPA peptide deviates from its position in the X-ray crystal structure. While the N- and C-terminal residues remain firmly anchored to the A- and F-pocket, respectively, central residues G4, R5, and A6 exhibit an average Cα RMSD of about 1.5 Å, consistent with the idea that residues E2 and F9 are the main determinants of peptide–MHCI stability. In line with this observation, even though the dehydrated environment of the crystal might naturally favor a tight binding groove, crystallographic Debye–Waller factors (B-factors) are slightly higher for residues 3–8 than they are for E1, E2, and F9. This peptide flexibility, observed in our simulations of MHCI in solution, could have implications for antigen recognition. MHC recognition by TCR, reviewed recently (61), involves three TCR loops (CDR1, 2, and 3); while CDR1 and 2 contact MHC, CDR3 is located directly atop the antigen peptide and thus interacts mostly with residues at the center of the peptide. For instance, interaction hotspots at positions 4, 5, and 8 in a 9-mer peptide have been identified in the structure of a TCR complexed with an MHCI molecule loaded with a cytomegalovirus antigen (62). Conformational changes observed in CDR3 upon binding to MHCIPL (61) were interpreted in terms of an induced fit mechanism that could be facilitated by small-scale rearrangements and fluctuations of the peptide central residues.
In all but one of the simulations performed in the present study, the peptide remained partially bound to B*44:02 and did not fully dissociate from the binding groove. However, in one of the Δ1C peptide simulations, the peptide C-terminus dissociated from the groove (after 5 ns), before the N-terminus also dissociated (after 255 ns). This trajectory was excluded from all further analyses. Investigating the possibility of complete peptide dissociation would require much more extensive simulations. To determine the thermodynamic equilibrium, one would need to fully sample both unbinding and rebinding events, which is clearly beyond the scope of the present work.
Partial dissociation is reversible, with truncated peptides mostly sampling conformations away from the binding groove, but also inside the groove to a lesser extent. Figure 3 shows RMSD timeseries from simulations of Δ1C and Δ2N, the two smallest deletions leading to partial dissociation of the central residues. For the C-terminally truncated peptide (Figure 3A), we observe three dissociations and two rebinding events; the peptide spends around 5% of the trajectory time with both termini bound. In the case of Δ2N (Figure 3B), the peptide spends roughly 35% of the trajectory time in the vicinity of the binding groove, although it never goes back to its initial position. Even though a similar behavior was observed in each of the five independent simulations, significantly longer trajectories would be required for statistically more accurate results. Larger peptide truncations (Δ2C, Δ3C, Δ3N) show a behavior similar to that of Δ2N, with the peptide sampling two sets of conformations, either away from or close to the groove.
Figure 3. Representative examples of partial dissociation events. (A) Cα RMSD timeseries of MHCI-bound Δ1C S8 (from its initial position in the X-ray crystal structure). The MHCIΔ1C structure shown is from a trajectory frame at t ≈ 0.35 µs; the corresponding point is labeled (*). (B) Cα RMSD timeseries of MHCIΔ2N DPA F3. The structure shown is from a trajectory frame at t ≈ 0.43 µs. Black lines show the reference MHCIPL system. 1-ns moving averages are overlaid on the timeseries (gray).
Taken together, our results suggest that MHCI loaded with low-affinity peptides can exist in two states: one where both peptide termini are near or inside the binding groove, and another where one terminus is displaced from the groove (see structures in Figure 3). These results are consistent with NMR experiments of MHCI molecules loaded with variant peptides (63), which showed that MHCI at the cell surface exist as a mixture of a major state and a minor state, the former being more hydrated. Since the same study showed that the melting temperature of MHCI increases with the prevalence of the less hydrated, minor state, we speculate that partial peptide dissociation might be the first step in the internalization and recycling of surface-exposed MHCI molecules.
3.2. Influence of Peptide Truncation on B*44:02 Structure and Dynamics
Simulation studies of MHCI often compare dynamics in the PL and PD forms (14, 15, 25, 33). However, in vivo, MHCI molecules are bound to (usually low-affinity) peptides after their synthesis and assembly with β2m and calnexin (39), before peptide editing takes place. We compared the dynamics of the PD and PL forms to those of MHCI loaded with various truncated peptides to assess if MHCIPD is an appropriate proxy for low-affinity peptides.
Figure 4 shows the width of the B*44:02 binding groove, measured at three positions, as a function of peptide cargo. Peptide release widens the binding groove by 2.5 Å at d1, measured near the F-pocket, and by 0.5 Å at d2, measured at the center of the groove. However, no significant change is observed at d3, near the A-pocket. In addition to this opening, the amplitude of the fluctuations of the groove dimensions is larger in the absence of peptide. Interestingly, the effects of N-terminal deletions are very different from those of full peptide removal. The groove remains mostly unchanged at d1 (F-pocket) and d2 (center), but narrows by 1 Å at d3 (A-pocket). This partial collapse suggests that MHCIPD is not an appropriate surrogate for studying the effects of N-terminal peptide editing. The effects of C-terminal deletions are more similar to those of peptide removal, albeit less pronounced. As more residues are removed, the binding groove widens at d1 (F-pocket), as expected, and fluctuations around the average increase at d2 (center); d3 (A-pocket) remains mostly unchanged. In conclusion, these results post a cautionary note when using MHCIPD to study peptide editing; while it might be warranted for C-terminal editing by tapasin (since the binding groove widens in the PD form as it does for the ΔC peptides), groove dynamics near the N-terminus depend not only on the loading state of the A-pocket, but also on that of the opposite F-pocket.
Figure 4. Influence of peptide cargo on B*44:02 binding groove width. Distances were averaged over all simulations for each system; standard deviation is shown as bars. Standard errors of the mean, obtained from a block averaging procedure (64), are below 0.5 Å in all cases. (A) Distance d1 is measured between the Cα atoms of Y85 and T138. (B) d2 between Y74 and V152. (C) d3 between Y59 and Y171.
To understand which domains of B*44:02 are influenced by peptide truncation and to what extent, we analyzed protein Cα RMSF for each system under study. These are shown along B*44:02 sequence in Figure 5 for the α1α2 binding groove domain, and in Figure 6 for the α3 domain. As expected, the binding groove is the region whose fluctuations are most impacted by peptide cargo modification. The α2−1 helix segment known to interact with tapasin (7, 14, 15) and the N-terminal portion of α1 located in the vicinity of the A-pocket are its most affected regions. Stretch 51–54, in α1 N-terminus, was proposed as a hinge that pivots in the transition from peptide-receptive to peptide-loaded MHCI (30, 31) in murine allele H2-Ld. In previous work (14, 15), we proposed a peptide-editing mechanism according to which peptides and tapasin compete for the α2−1 region. The results presented here yield a refined view, in the sense that they suggest there could be two sensing mechanisms at work to distinguish properly loaded MHCI at the A- and F-pockets. The fluctuations in the two affected regions are independent: N-terminal peptide truncation recovers the increased RMSF in α1 (Figure 5B), while C-terminal truncation yields increased fluctuations in α2−1 (Figure 5C).
Figure 5. Influence of peptide cargo on B*44:02 α1α2 fluctuations. Cα RMSF along sequence and the absolute difference compared to MHCIPL are shown for each system. ΔRMSF are also mapped onto the structure of α1α2. Peptide not shown for clarity. (A) MHCIPD. (B) MHCIΔ1N–Δ3N. (C) MHCIΔ1C–Δ3C.
Figure 6. Influence of peptide cargo on B*44:02 α3 fluctuations. Cα RMSF along sequence and the absolute difference compared to MHCIPL are shown for each system. ΔRMSF are also mapped onto the structure of α3. Peptide not shown for clarity. (A) MHCIPD. (B) MHCIΔ1N–Δ3N. (C) MHCIΔ1C–Δ3C.
Possible long-range effects of cargo modification on α3 are important to consider given that this domain contacts tapasin in the PLC and that allosteric effects involving α3 have been proposed (24) as a mechanism for peptide sensing by tapasin. However, apart from the Δ1N peptide (see below), we could not observe any significant effect of peptide truncation on α3 fluctuations (Figure 6). It is possible, however, that signal transduction events are slow compared to the length of our trajectories (1.0 µs) and that increased sampling could yield a different picture. Another possibility is that subtle conformational changes in α3 are not reflected in the magnitude of its fluctuations.
3.3. Configurational Entropy As a Function of Peptide Cargo
To complement our analysis of fluctuations, we computed configurational entropies for the α1α2 and α3 domains of B*44:02. Results are shown in Figure 7, where ΔS is the difference between the entropy of each system and the entropy of the reference system MHCIPL. The important increase of entropy (100–300 J/[K mol]) observed for α1α2 (Figure 7A) is consistent with the increased fluctuations of the individual residues, as discussed above. The small statistical uncertainties (less than ±10 J/[K mol] for most systems) indicate convergence, i.e., that 5.0 µs of MD sampling is sufficient to assess B*44:02 α1α2 dynamics.
Figure 7. Influence of peptide cargo on B*44:02 configurational entropy. ΔS for each system is the difference to MHCIPL computed for the given domain. (A) α1α2. (B) α3. Statistical error estimates are the difference between ΔS computed from the full trajectories to the value obtained when discarding a random 10% of the data, and are shown only when larger than 10 J/(K mol).
Configurational entropy differences in the α3 domain (Figure 7B) are also consistent with residue RMSF: ΔS is very small for most systems. A modest but significant increase (63 J/[K mol]) is observed, however, for Δ1N. While this could be an indication of a signal being transmitted between the binding groove and α3, it also raises the question of why this difference is not observed for any other N- or C-terminal deletion. One possibility is that, since allostery-related protein dynamics often take place on the microsecond–millisecond timescale, our simulations lack the sampling necessary to describe such slow events. Another possibility is that a certain stiffness of the binding groove is required for signal transduction, and that the increased α1α2 plasticity caused by deletions larger than a single residue may disrupt this sensing mechanism; however, this would not explain the apparent absence of allosteric signal in the Δ1C system. Finally, a third possibility is that the presence of tapasin, the tapasin–ERp57 conjugate, or other accessory proteins forming the PLC, is necessary to induce conformational changes in MHCI that mediate the allosteric mechanism.
3.4. Contributions of the Peptide Side Chains and Backbone
Antigenic peptide specificity is often driven by their N- and C-terminal residues (65), which also determine their stability in the MHCI binding groove. In B*44:02DPA, residues E2 and F9 extend their side chains toward the floor of the binding groove, occupying the A- and F-pockets, respectively. Since their truncation leads to partial peptide dissociation in our simulations (see Figure 2), we expect their side chains to be key determinants of peptide stability. To test this hypothesis, we have carried out additional simulations of B*44:02-bound modified DPA peptides. These have had one or two residues replaced by glycines at either the N- or the C-terminus. We denote these systems MHCIδ1N, MHCIδ2N, and MHCIδ1C. As for the other modified peptides, trajectories totaling 5 µs of MD simulations (5 × 1 µs) were acquired for each system.
Partial peptide dissociation was observed for two of the three systems under study (Figure 8). Interestingly, side chain removal at the N-terminus leads to a partial dissociation pattern that is significantly different from the one observed for the corresponding truncated peptides. δ1N does not dissociate at all, showing a stable binding almost identical to that of unmodified DPA. δ2N exhibits partial dissociation, but to a much lesser extent than Δ2N: only E1 and E2 vacate the binding groove. Therefore, at the N-terminus, both the peptide main chain and the E1 and E2 side chains contribute to MHCI binding at the A-pocket. Conversely, at the C-terminus, the δ1C dissociation pattern is similar to that of Δ1C, with residues 4–9 dissociating from the groove. The F9 side chain is thus a critical determinant of peptide binding at the F-pocket.
Figure 8. Partial dissociation of peptides with glycine substitutions from B*44:02. Cα RMSD from the initial position in the X-ray crystal structure were averaged for each peptide residue over the five independent 1-µs trajectories acquired for each system. The dotted line represents a chosen dissociation threshold at 3 Å.
3.5. Tapasin Widens the Binding Groove of MHCI and Weakens Peptide Binding Forces
We performed 5 additional, independent, 1-µs MD simulations of MHCIΔ1C in complex with tapasin, using our previously determined structure of tapasin–MHCI (14, 15) and truncating the full-length DPA peptide present in that structure. This allows us to compare the dynamics of B*44:02 loaded with a low-affinity peptide in the free and tapasin-complexed forms. Results (Figure 9) show that the binding groove width near the F-pocket (d1, as shown previously in Figure 4), is wider by about 1.0 Å in the presence of tapasin. This is consistent with the “tug-of-war” mechanism we previously described (15), according to which tapasin selects peptides depending on their C-terminus by pulling MHCI β7 underneath α2−1 to open the groove and thus accelerate peptide release.
Figure 9. Effect of tapasin on MHCI-binding groove width near the F-pocket. Distance d1 is measured between the Cα atoms of Y85 and T138 as in Figure 4. Distances were averaged over all simulations for each system; standard deviation is shown as bars. Standard errors of the mean, obtained from a block averaging procedure (64), are below 0.5 Å in all cases. The tapasin–MHCI controls for the PL and PD forms are taken from previous simulations (15). Tapasin TC and MHCI α3 not shown for clarity.
Next, to obtain a more detailed and spatially resolved picture, we used force distribution analysis (FDA) to compute the residue-residue pairwise forces between B*44:02 and the Δ1C peptide in tapasin-complexed and -free MHCI (Figure 10). Overall, forces between MHCI and the peptide are on average less attractive (by 150 pN) in the tapasin–MHCI complex, showing that the opening of the groove induced by tapasin decreases the affinity of MHCI for its low-affinity cargo. Interestingly, considering individual pairwise forces at the residue level reveals that tapasin affects the peptide not just in the vicinity of the α2−1 region (C-terminus), but along its whole sequence. Indeed, the largest force difference is observed between the peptide at position R5 and B*44:02 R97 at the bottom of the binding groove; in the tapasin–MHCI complex, this force is more repulsive (by 148 pN) than observed for B*44:02 in the free form. A less attractive binding force is also observed at the peptide N-terminus (E1). Thus, tapasin-induced widening of the binding groove near the F-pocket (see Figure 9) has long-ranging effects that destabilize peptide binding, increasing the latter’s off-rate. Full peptide release, however, was not observed in the present simulations. Much longer timescales would likely be required to investigate the tapasin-catalyzed unbinding event itself.
Figure 10. Effects of tapasin on the intermolecular forces in the binding groove of MHCIΔ1C. Differences (ΔF) in the pairwise residue-residue forces between tapasin-complexed and -free MHCI are shown as cylinders scaled according to force difference magnitude, and show the effects of transitioning MHCIΔ1C from the tapasin-free form to the tapasin-complexed one. Force differences range from 30 to 150 pN; smaller differences are statistically insignificant, and therefore not shown. MHCI α3 not shown for clarity.
Binding groove widening has been suggested before (9) in a biochemical study of tapasin–MHCI. In their mechanism of complex assembly/disassembly, the authors proposed two possible conformational transitions taking place in the MHCI binding groove. The first is an equilibrium between “closed” and “open” states, which exist both in presence and absence of tapasin and could correspond to a tightly packed and loose binding groove, respectively (9, 24, 66). The second is a switch between a canonical and a “tapasin-disrupted” binding groove conformation, which could correspond to the widening we observe in our simulations. Partial peptide dissociation of the Δ1C peptide, in turn, would prevent the transition from the “open” to the “closed” state. Interestingly, the same study (9) also found that mutating certain high-affinity peptides at their N-terminus caused them to exhibit higher tapasin sensitivity. This indicates that tapasin impacts MHCI across the full length of the peptide, which is consistent with our above FDA results that reveal destabilizing effects at E1 and R5.
Using comparative MD simulations of B*44:02 loaded with a variety of truncated DPA peptides, we have shown that positions 2 and 9 are crucial for peptide stability in the binding groove; the corresponding Δ2N and Δ1C are the smallest deletions leading to partial dissociation of the peptide central residues on the microsecond timescale. Even in the full-length antigen, however, central peptide residues enjoy substantial flexibility, potentially facilitating recognition by the TCR CDR3 loop. Transient rebinding to the groove indicates that at least some MHCI molecules exist in an equilibrium of two populations: one with a tightly bound peptide and the other with a partially dissociated peptide. These results are consistent with NMR experiments on allele B*35:01 in complex with different variant peptides (63), and could indicate the first step leading to the internalization and recycling of MHCI. Furthermore, partial dissociation from MHCI at the peptide N-terminus could explain how ERAAP can trim MHCI-bound peptides to their correct final length (19). Finally, partial binding at the N-terminus is likely the first step in dipeptide-catalyzed peptide exchange (36), where the MHCI F-pocket is occupied by short dipeptides that are rapidly replaced upon the addition of high-affinity peptides.
Peptide truncation impacts B*44:02 dynamics in a different manner than complete peptide removal does, especially for N-terminal deletions. This suggests that MHCIPD is a sub-optimal surrogate for MHCI bound to low-affinity peptides, something that needs to be considered for simulation studies of peptide editing, especially N-terminal editing by ERAAP. MHCIΔC systems, on the other hand, are more similar to the PD form as far as binding groove dynamics in the vicinity of the F-pocket are concerned. Fluctuation and entropy analyses confirm the importance of the α2−1 helix segment bordering the F-pocket of the binding groove. Increased flexibility (as observed for the ΔC systems) is necessary for C-terminal peptide editing by the tapasin TN domain as described previously (15): tapasin and the antigen compete for α2−1 to open and close the groove, respectively, leading to accelerated peptide release until the binding of a high-affinity antigen.
MHCI B*44:05, which differs from B*44:02 by only one residue (Y116 instead of D), can efficiently load antigens without tapasin. Since Y116 is located on the floor on the binding groove (β6) underneath the α2−1 helix segment, it is reasonable to expect it to influence the dynamics of that region. Indeed, comparative MD simulation studies (27) of these two alleles have found that B*44:05 exhibits less fluctuations in the vicinity of the F-pocket. By enforcing partial dissociation of the N- and C-terminal residues through umbrella sampling, the authors found that the energy barrier required for peptide dissociation at the C-terminus is higher in B*44:05 (tapasin-independent) than in B*44:02 (tapasin-dependent). A similar conclusion was drawn for alleles B*27:05 and B*27:09 (32), that also differ by a single amino acid in the F-pocket: tapasin-dependent B*27:05 exhibits higher conformational flexibility in the vicinity of the α2−1 region. Therefore, dependence on tapasin could stem, at least for these two alleles, from a requirement for its chaperone activity, which confines MHCI to a peptide-receptive conformation. The widening of the binding groove, conversely, relates apparently only to the catalysis of peptide dissociation. Future work involving the simulation of tapasin-independent allele B*44:05 with modified or truncated peptides is required to determine if partial peptide dissociation at the C-terminus follows a similar pattern as in B*44:02.
We found no strong indication of allosteric effects of peptide truncation on the CD8 loop contacting tapasin. We did, however, observe a possible modulation of α3 dynamics in the Δ1N peptide system. More extensive sampling might be required to elucidate α3 dynamics and reconcile MD simulations with the computational systems models of tapasin function that were recently proposed (24). If C-terminal peptide editing proceeds by a dual mechanism, the next pertinent question would be what is the purpose of each mechanism. One possible hypothesis is that the TN–α1α2 interface provides the mechanical force to open the groove, and therefore accelerates peptide release; the TC–α3 one is a signal for tapasin to recognize MHCI loaded with low-affinity cargo and disengage once a high-affinity antigen is bound. Extended MD simulations of the tapasin–MHCI complex, comparing the effects of low- and high-affinity peptides, could be used to test this hypothesis at the atomic level.
OF and LS designed the study. OF, SW, and LS analyzed and interpreted the results and revised the manuscript. OF performed the simulations and drafted the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We wish to thank Sebastian Springer (Jacobs University Bremen) for stimulating discussions. The authors are grateful to the Gauss Centre for Supercomputing e.V. (www.gauss-centre.eu) for providing computing time on the GCS supercomputer SuperMUC at the Leibniz Supercomputing Centre (www.lrz.de).
This work was funded by the Deutsche Forschungsgemeinschaft (DFG) through Cluster of Excellence RESOLV (EXC 1069) and Emmy-Noether grant SCHA-1574 to LS.
6. Momburg F, Tan P. Tapasin: the keystone of the loading complex optimizing peptide binding by MHC class I molecules in the endoplasmic reticulum. Mol Immunol (2002) 39(3):217–33. doi:10.1016/S0161-5890(02)00103-7
7. Dong G, Wearsch PA, Peaper DR, Cresswell P, Reinisch KM. Insights into MHC class I peptide loading from the structure of the tapasin-ERp57 thiol oxidoreductase heterodimer. Immunity (2009) 30(1):21–32. doi:10.1016/j.immuni.2008.10.018
10. Schneeweiss C, Garstka M, Smith J, Hütt M-T, Springer S. The mechanism of action of tapasin in the peptide exchange on MHC class I molecules determined from kinetics simulation studies. Mol Immunol (2009) 46(10):2054–63. doi:10.1016/j.molimm.2009.02.032
11. Kanaseki T, Lind KC, Escobar H, Nagarajan N, Reyes-Vargas E, Rudd B, et al. ERAAP and tapasin independently edit the amino and carboxyl termini of MHC class I peptides. J Immunol (2013) 191(4):1547–55. doi:10.4049/jimmunol.1301043
13. Dalchau N, Phillips A, Goldstein LD, Howarth M, Cardelli L, Emmott S, et al. A peptide filtering relation quantifies MHC class I peptide optimization. PLoS Comput Biol (2011) 7(10):e1002144. doi:10.1371/journal.pcbi.1002144
14. Fleischmann G, Fisette O, Thomas C, Wieneke R, Tumulka F, Schneeweiss C, et al. Mechanistic basis for epitope proofreading in the peptide-loading complex. J Immunol (2015) 195(9):4503–13. doi:10.4049/jimmunol.1501515
18. Mpakali A, Giastas P, Mathioudakis N, Mavridis IM, Saridakis E, Stratikos E. Structural basis for antigenic peptide recognition and processing by endoplasmic reticulum (ER) aminopeptidase 2. J Biol Chem (2015) 290(43):26021–32. doi:10.1074/jbc.M115.685909
19. Chen H, Li L, Weimershaus M, Evnouchidou I, van Endert P, Bouvier M. ERAP1-ERAP2 dimers trim MHC I-bound precursor peptides; implications for understanding peptide editing. Sci Rep (2016) 6:28902. doi:10.1038/srep28902
20. Hein Z, Uchtenhagen H, Abualrous ET, Saini SK, Janßen L, Van Hateren A, et al. Peptide-independent stabilization of MHC class I molecules breaches cellular quality control. J Cell Sci (2014) 127(Pt 13):2885–97. doi:10.1242/jcs.145334
21. Suh WK, Derby MA, Cohen-Doyle MF, Schoenhals GJ, Früh K, Berzofsky JA, et al. Interaction of murine MHC class I molecules with tapasin and TAP enhances peptide loading and involves the heavy chain alpha3 domain. J Immunol (1999) 162(3):1530–40.
22. Turnquist HR, Vargas SE, Schenk EL, McIlhaney MM, Reber AJ, Solheim JC. The interface between tapasin and MHC class I: identification of amino acid residues in both proteins that influence their interaction. Immunol Res (2002) 25(3):261–9. doi:10.1385/IR:25:3:261
23. Bailey A, van Hateren A, Elliott T, Werner JM. Two polymorphisms facilitate differences in plasticity between two chicken major histocompatibility complex class I proteins. PLoS One (2014) 9(2):e89657. doi:10.1371/journal.pone.0089657
25. Zacharias M, Springer S. Conformational flexibility of the MHC class I alpha1-alpha2 domain in peptide bound and free states: a molecular dynamics simulation study. Biophys J (2004) 87(4):2203–14. doi:10.1529/biophysj.104.044743
27. Sieker F, Straatsma TP, Springer S, Zacharias M. Differential tapasin dependence of MHC class I molecules correlates with conformational changes upon peptide dissociation: a molecular dynamics simulation study. Mol Immunol (2008) 45(14):3714–22. doi:10.1016/j.molimm.2008.06.009
28. Garstka MA, Fritzsche S, Lenart I, Hein Z, Jankevicius G, Boyle LH, et al. Tapasin dependence of major histocompatibility complex class I molecules correlates with their conformational flexibility. FASEB J (2011) 25(11):3989–98. doi:10.1096/fj.11-190249
29. Narzi D, Becker CM, Fiorillo MT, Uchanska-Ziegler B, Ziegler A, Böckmann RA. Dynamical characterization of two differentially disease associated MHC class I proteins in complex with viral and self-peptides. J Mol Biol (2012) 415(2):429–42. doi:10.1016/j.jmb.2011.11.021
30. Mage MG, Dolan MA, Wang R, Boyd LF, Revilleza MJ, Robinson H, et al. The peptide-receptive transition state of MHC class I molecules: insight from structure and molecular dynamics. J Immunol (2012) 189(3):1391–9. doi:10.4049/jimmunol.1200831
31. Mage MG, Dolan MA, Wang R, Boyd LF, Revilleza MJ, Robinson H, et al. A structural and molecular dynamics approach to understanding the peptide-receptive transition state of MHC-I molecules. Mol Immunol (2013) 55(2):123–5. doi:10.1016/j.molimm.2012.10.021
32. Abualrous ET, Fritzsche S, Hein Z, Al-Balushi MS, Reinink P, Boyle LH, et al. F pocket flexibility influences the tapasin dependence of two differentially disease-associated MHC class I proteins. Eur J Immunol (2015) 45(4):1248–57. doi:10.1002/eji.201445307
33. Ostermeir K, Springer S, Zacharias M. Coupling between side chain interactions and binding pocket flexibility in HLA-B*44:02 molecules investigated by molecular dynamics simulations. Mol Immunol (2015) 63(2):312–9. doi:10.1016/j.molimm.2014.07.021
34. Kongkaew S, Yotmanee P, Rungrotmongkol T, Kaiyawet N, Meeprasert A, Kaburaki T, et al. Molecular dynamics simulation reveals the selective binding of human leukocyte antigen alleles associated with Behçet’s disease. PLoS One (2015) 10(9):e0135575. doi:10.1371/journal.pone.0135575
36. Saini SK, Schuster H, Ramnarayan VR, Rammensee H-G, Stevanović S, Springer S. Dipeptides catalyze rapid peptide exchange on MHC class I molecules. Proc Natl Acad Sci U S A (2015) 112(1):202–7. doi:10.1073/pnas.1418690112
37. Abualrous ET, Saini SK, Ramnarayan VR, Ilca FT, Zacharias M, Springer S. The carboxy terminus of the ligand peptide determines the stability of the MHC class I molecule H-2Kb: a combined molecular dynamics and experimental study. PLoS One (2015) 10(8):e0135421. doi:10.1371/journal.pone.0135421
40. Macdonald WA, Purcell AW, Mifsud NA, Ely LK, Williams DS, Chang L, et al. A naturally selected dimorphism within the HLA-B44 supertype alters class I structure, peptide repertoire, and T cell recognition. J Exp Med (2003) 198(5):679–91. doi:10.1084/jem.20030066
41. Olsson MHM, Søndergaard CR, Rostkowski M, Jensen JH. PROPKA3: consistent treatment of internal and surface residues in empirical pKa predictions. J Chem Theory Comput (2011) 7(2):525–37. doi:10.1021/ct100578z
42. Pronk S, Páll S, Schulz R, Larsson P, Bjelkmar P, Apostolov R, et al. GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics (2013) 29(7):845–54. doi:10.1093/bioinformatics/btt055
43. Páll S, Abraham M, Kutzner C, Hess B, Lindahl E. Tackling exascale software challenges in molecular dynamics simulations with GROMACS. In: Markidis S, Laure E, editors. Solving Software Challenges for Exascale. (Vol. 8759), New York, NY: Springer International Publishing (2015). p. 3–27.
44. Cornell WD, Cieplak P, Bayly CI, Gould IR, Merz KM, Ferguson DM, et al. A second generation force field for the simulation of proteins, nucleic acids, and organic molecules. J Am Chem Soc (1995) 117(19):5179–97. doi:10.1021/ja00124a002
45. Wang J, Cieplak P, Kollman PA. How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? J Comput Chem (2000) 21(12):1049–74. doi:10.1002/1096-987X(200009)21:12<1049:AID-JCC3>3.0.CO;2-F
46. Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C. Comparison of multiple amber force fields and development of improved protein backbone parameters. Proteins (2006) 65(3):712–25. doi:10.1002/prot.21123
48. Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, Dror RO, et al. Improved side-chain torsion potentials for the amber ff99SB protein force field. Proteins (2010) 78(8):1950–8. doi:10.1002/prot.22711
49. Aliev AE, Kulke M, Khaneja HS, Chudasama V, Sheppard TD, Lanigan RM. Motional timescale predictions by molecular dynamics simulations: case study using proline and hydroxyproline sidechain dynamics. Proteins (2014) 82(2):195–215. doi:10.1002/prot.24350
53. Feenstra KA, Hess B, Berendsen HJC. Improving efficiency of large time-scale molecular dynamics simulations of hydrogen-rich systems. J Comput Chem (1999) 20(8):786–98. doi:10.1002/(SICI)1096-987X(199906)20:8<786:AID-JCC5>3.0.CO;2-B
61. Rossjohn J, Gras S, Miles JJ, Turner SJ, Godfrey DI, McCluskey J. T cell antigen receptor recognition of antigen-presenting molecules. Annu Rev Immunol (2015) 33:169–200. doi:10.1146/annurev-immunol-032414-112334
62. Gras S, Saulquin X, Reiser J-B, Debeaupuis E, Echasserieau K, Kissenpfennig A, et al. Structural bases for the affinity-driven selection of a public TCR against a dominant human cytomegalovirus epitope. J Immunol (2009) 183(1):430–7. doi:10.4049/jimmunol.0900556
63. Yanaka S, Ueno T, Shi Y, Qi J, Gao GF, Tsumoto K, et al. Peptide-dependent conformational fluctuation determines the stability of the human leukocyte antigen class I complex. J Biol Chem (2014) 289(35):24680–90. doi:10.1074/jbc.M114.566174
Keywords: major histocompatibility complex class I (MHCI), peptide loading complex (PLC), antigen, peptide editing, protein dynamics, molecular dynamics (MD) simulations
Citation: Fisette O, Wingbermühle S and Schäfer LV (2017) Partial Dissociation of Truncated Peptides Influences the Structural Dynamics of the MHCI Binding Groove. Front. Immunol. 8:408. doi: 10.3389/fimmu.2017.00408
Received: 30 November 2016; Accepted: 22 March 2017;
Published: 18 April 2017
Edited by:Tim Elliott, University of Southampton, UK
Reviewed by:Marlene Bouvier, University of Illinois at Chicago, USA
Andras Szilagyi, Hungarian Academy of Sciences, Hungary
Copyright: © 2017 Fisette, Wingbermühle and Schäfer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Olivier Fisette, firstname.lastname@example.org