Immune recognition of infected or malignantly transformed cells relies on antigenic peptides exposed at the cell surface by major histocompatibility complex class I (MHC I) molecules. Selection and loading of peptides onto MHC I is orchestrated by the peptide-loading complex (PLC), a multiprotein assembly whose structure has not yet been resolved. Tapasin, a central component of the PLC, stabilises MHC I and catalyses the exchange of low-affinity against high-affinity, immunodominant peptides. Up to now, the molecular basis of this peptide editing mechanism remained elusive. Here, using all-atom molecular dynamics (MD) simulations, we unravel the atomic details of how tapasin and antigen peptides act on the MHC I binding groove. Force distribution analysis reveals an intriguing molecular tug-of-war mechanism: only high-affinity peptides can exert sufficiently large forces to close the binding groove, thus overcoming the opposite forces exerted by tapasin to open it. Tapasin therefore accelerates the release of low-affinity peptides until a high-affinity antigen binds, promoting subsequent PLC break-down. Fluctuation and entropy analyses show how tapasin chaperones MHC I by stabilising it in a peptide-receptive conformation. Our results explain previous experiments and mark a key step towards a better understanding of adaptive immunity.

Molecular simulations of intrinsically disordered proteins (IDPs) are challenging because they require sampling a very large number of relevant conformations, corresponding to a multitude of shallow minima in a flat free energy landscape. However, in the presence of a binding partner, the free energy landscape of an IDP can be dominated by few deep minima. This characteristic imposes high demands on the accuracy of the force field used to describe the molecular interactions. Here, as a model system for an IDP that is unstructured in solution but folds upon binding to a structured interaction partner, the transactivation domain of c-Myb was studied both in the unbound (free) form and when bound to the KIX domain. Six modern biomolecular force fields were systematically tested and compared in terms of their ability to describe the structural ensemble of the IDP. The protein force field/water model combinations included in this study are AMBER ff99SB-disp with its corresponding water model that was derived from TIP4P-D, CHARMM36m with TIP3P, ff15ipq with SPC/Eb, ff99SB*-ILDNP with TIP3P and TIP4P-D, and FB15 with TIP3P-FB water. Comparing the results from REST2-enhanced sampling simulations with experimental CD spectra and secondary chemical shifts reveals that the ff99SB-disp force field can realistically capture the broad and mildly helical structural ensemble of free c-Myb. The structural ensembles yielded by CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 are also mildly helical; however, each of these force fields can be assigned a specific subset of c-Myb residues for which the simulations could not reproduce the experimental secondary chemical shifts. In addition, microsecond-timescale MD simulations of the KIX/c-Myb complex show that most force fields used preserve a stable helix fold of c-Myb in the complex. Still, all force fields predict a KIX/c-Myb complex interface that differs slightly from the structures provided by NMR because several NOE-derived distances between KIX and c-Myb were exceeded in the simulations. Taken together, the ff99SB-disp force field in the first place but also CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 can be suitable choices for future simulation studies of the coupled folding and binding mechanism of the KIX/c-Myb complex and potentially also other IDPs.

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.

Enhanced sampling techniques are a promising approach to obtain reliable binding free-energy profiles for flexible protein−ligand complexes from molecular dynamics (MD) simulations. To put four popular enhanced sampling techniques to a biologically relevant and challenging test, we studied the partial dissociation of an antigenic peptide from the Major Histocompatibility Complex I (MHC I) HLA-B*35:01 to systematically investigate the performance of umbrella sampling (US), replica exchange with solute tempering 2 (REST2), bias exchange umbrella sampling (BEUS, or replica-exchange umbrella sampling), and well-tempered metadynamics (MTD). With regard to the speed of sampling and convergence, the peptide-MHC I complex (pMHC I) under study showcases intrinsic strengths and weaknesses of the four enhanced sampling techniques used. We found that BEUS can best handle the sampling challenges that arise from the coexistence of an enthalpically and an entropically stabilized free-energy minimum in the pMHC I under study. These findings might also be relevant for other flexible biomolecular systems with competing enthalpically and entropically stabilized minima.

Photoswitchable cages that confine small guest molecules inside their cavities offer a way to control the binding/unbinding process through irradiation with light of different wavelengths. However, a detailed characterization of...