Molecular dynamics simulations to provide insights into epitopes coupled to the soluble and membrane-bound MHC-II complexes

PLoS One. 2013 Aug 19;8(8):e72575. doi: 10.1371/journal.pone.0072575. eCollection 2013.

Abstract

Epitope recognition by major histocompatibility complex II (MHC-II) is essential for the activation of immunological responses to infectious diseases. Several studies have demonstrated that this molecular event takes place in the MHC-II peptide-binding groove constituted by the α and β light chains of the heterodimer. This MHC-II peptide-binding groove has several pockets (P1-P11) involved in peptide recognition and complex stabilization that have been probed through crystallographic experiments and in silico calculations. However, most of these theoretical calculations have been performed without taking into consideration the heavy chains, which could generate misleading information about conformational mobility both in water and in the membrane environment. Therefore, in absence of structural information about the difference in the conformational changes between the peptide-free and peptide-bound states (pMHC-II) when the system is soluble in an aqueous environment or non-covalently bound to a cell membrane, as the physiological environment for MHC-II is. In this study, we explored the mechanistic basis of these MHC-II components using molecular dynamics (MD) simulations in which MHC-II was previously co-crystallized with a small epitope (P7) or coupled by docking procedures to a large (P22) epitope. These MD simulations were performed at 310 K over 100 ns for the water-soluble (MHC-IIw, MHC-II-P(7w), and MHC-II-P(22w)) and 150 ns for the membrane-bound species (MHC-IIm, MHC-II-P(7m), and MHC-II-P(22m)). Our results reveal that despite the different epitope sizes and MD simulation environments, both peptides are stabilized primarily by residues lining P1, P4, and P6-7, and similar noncovalent intermolecular energies were observed for the soluble and membrane-bound complexes. However, there were remarkably differences in the conformational mobility and intramolecular energies upon complex formation, causing some differences with respect to how the two peptides are stabilized in the peptide-binding groove.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Cell Membrane / metabolism*
  • Epitopes / chemistry
  • Epitopes / immunology*
  • Histocompatibility Antigens Class II / chemistry*
  • Histocompatibility Antigens Class II / immunology*
  • Molecular Dynamics Simulation*
  • Peptides / chemistry
  • Peptides / metabolism
  • Protein Binding
  • Protein Conformation
  • Solubility
  • Thermodynamics
  • Water / chemistry

Substances

  • Epitopes
  • Histocompatibility Antigens Class II
  • Peptides
  • Water

Grants and funding

The authors thank Instituto de ciencia y tecnología del DF (PIRIVE09-9), Consejo Nacional de Ciencia y Tecnología (CONACYT), and programa integral de fortalecimiento institucional-secretaria de investif¿gación y posgrado-Comisión de Operación y Fomento de Actividades Académicas/Instituto Politecnico Nacional for financial support. MB thanks CONACYT for scholarship. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.